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ABSTRACT 

Consider  the  nonlinear  Volterra  equation 

(V)  u( t )  +  (b*Au) ( t)  a  f(t)  <0  <  t  <  «) 

in  the  general  setting  b  :  [0,“)  -►  R  a  given  kernel,  A  a  nonlinear  m- 
accretive  operator  on  a  real  Banach  space  X,  f  :  [0,“)  *  X  a  given  function, 
and  *  the  convolution.  This  paper,  based  on  lectures  delivered  at  West 
Virginia  University,  discusses  existing  and  recent  results  for  the  following 
problems  concerning  (V):  JL.  the  global  existence  and  uniqueness  of  solutions 
and  their  continuous  dependence  on  the  data,  _2 •  the  boundedness  and 
asymptotic  behaviour  as  t  -*•  “  in  the  special  cases  when  X  =  H  is  a  real 
Hilbert  space  and  A  is  either  a  maximal  monotone  operator  on  H  or  A  is  a 
subdifferential  of  a  proper,  convex,  lower  semicontinuous  function 
<P  ;  H  -*•  (-*,+“],  3_.  the  existence,  boundedness,  and  asymptotic  behaviour  of 

positive  solutions  in  the  general  setting.  The  theory  is  used  to  study  one 
possible  model  problem  for  heat  flow  in  a  material  with  "memory"  which  can  be 
transformed  to  the  equivalent  form  (V)  under  physically  reasonable 
assumptions;  the  latter  provide  a  motivation  for  the  natural  setting  of  much 
of  the  theory  developed  here.  This  and  various  other  models  for  heat  flow  in 
such  materials  are  formulated  from  physical  principals  and  discussed  in  an 
introductory  chapter. 
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SIGNIFICANCE  AND  EXPLANATION 


Thege  notes  on  lectures  delivered  at  West  Virginia  University  preceding 
the  conference  on  Volterra  and  Functional  Differential  Equations  in  June  1979, 
concern  the  existence,  uniqueness,  positivity,  boundedness,  and  asymptotic 
behaviour  as  t  ♦  •  of  solutions  of  the  nonlinear  Volterra  equation 

(V)  u(t)  +  (b*Au)(t)  a  f(t)  (0  <  t  <  «)  . 

The  general  setting  for  V  is  as  follows:  b  :  (0,"»)  ♦  R  is  a  given 
kernel,  A  is  a  m-accretive,  possibly  multivalued,  operator  on  a  real  Banach 
space  X,  f  :  (0,m)  +  X  is  a  given  function,  and  *  denotes  the 
convolution;  the  integral  in  V  is  taken  in  the  sense  of  Bochner.  The 
special  cases  of  A  maximal  monotone  on  a  real  Hilbert  space  H,  and 
A  =  3^,  the  subdifferential  of  a  proper,  convex,  lower  semicontinuous 
function  ip  i  H  *  (-»,<*>]  will  also  play  a  prominent  role,  primarily  in  the 
boundedness  and  asymptotic  theory  for  (V) • 

It  should  be  observed  that  if  b  =  1  and  f  e  W1 ' 1(0,“»;X) ,  where  w 
denotes  the  usual  Sobolev  space,  equation  (V)  is  formally  equivalent  to  the 
evolution  problem 

(E)  $  +  Au  »  f •  (0  <  t  <  •>,  u(0 )  -  uft  -  f(0)  . 

at  0 

Thus  the  theory  for  (V)  is  to  a  considerable  extent  a  generalization  of  the 
theory  of  evolution  equations,  and  uses  most  of  the  techniques  for  the  latter 
combined  with  techniques  for  Volterra  equations  developed  in  recent  years. 

Chapter  1  is  primarily  intended  for  motivation.  Beginning  from  simple 
physical  principles  equation  (V)  is  derived  as  one  possible  mathematical  model 
for  nonlinear  heat  flow  in  a  homogeneous  body  of  material  with  memory 
following  ideas  of  B.  D.  Coleman,  N.  E.  Gurtin,  R.  C.  Mac  Camy,  W.  Noll, 

J.  H.  Nunziato,  and  A.  C.  Pipkin.  While  the  derivation  is  restricted  to  one 
space  dimension,  |he  modification  for  heat  flow  in  a  homogeneous  body 
0  in  R  or  R  of  isotropic  material  is  also  indicated.  One  purpose  of 
Chapter  1  is  to  arrive  at  a  physically  reasonable  set  of  conditions  concerning 
the  kernel  b  in  (V),  under  which  one  can  expect  boundedness  and  the  kind  of 
asymptotic  behaviour  of  solutions  of  (V)  solutions  studied  later  in  Chapter 
3.  The  kernel  b  in  (V)  does  not  arise  directly  from  physical  principles; 
rather,  it  is  the  case  that  b  is  expressed  in  terms  of  two  physically 
measurable  quantities  (at  least  in  principle)  about  which  one  can  make 
appropriate  assumptions  based  on  physical  considerations.  The  types  of 
assumptions  made  concerning  the  operator  A  and  the  function  f  are  also 
motivated.  In  particular,  A  is  a  nonlinear  second-order  elliptic  partial 
differential  operator  in  the  space  variables  which  incorporates  the  boundary 
conditions,  and  f  depends  on  the  external  heat  supply,  the  initial 
temperature  distribution  in  the  body  and  the  history  of  temperature  in  the 
body.  If  b  =  1 ,  then  in  the  application  (V)  is  equivalent  to  the  classical 
nonlinear  heat  equation  in  an  ordinary  body. 

Another  purpose  of  Chapter  1  is  to  point  out  that  heat  flow  in  certain 
materials  with  memory  can  also  be  modelled  by  integrodifferential  equations 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  mrc,  and  not  with  the  author  of  this  report. 


other  than  (V).  Limitations  of  space  do  not  permit  us  to  discuss  the 
qualitative  aspects  of  these  alternative  models  in  any  detail;  however, 
references  to  relevant  existing  literature  are  given. 


Chapter  2  discusses  the  theory  of  existence,  uniqueness,  and  continuous 
dependence  of  global  solutions  of  (V),  both  in  the  general  case  of  a  real 
Banach  space  and  in  the  special  case  of  a  real  Hilbert  space,  and  for  A  m 
where  the  results  are  stronger.  The  development  is  primarily  based  on  recent 
joint  work  with  M.  G.  Crandall  (26) ,  and  partly  on  a  recent  paper  by 
G.  Gripenberg  [34].  References  to  earlier  and  related  literature  are  given. 

Chapter  3  develops  the  theory  of  boundedness  and  asymptotic  behaviour  of 
solutions  of  (V)  as  t  ♦  *,  under  assumptions  partly  motivated  by  the  heat 
flow  problem  formulated  in  Chapter  1;  application  of  the  theory  to  this 
problem  is  given.  The  development  is  based  on  forthcoming  joint  work  with 
P.  Clement  and  R.  C.  Mac  Camy  [20].  References  to  other  pertinent  literature 
are  given. 

Chapter  4,  based  on  recent  and  forthcoming  joint  work  with  P.  Clement, 
[18],  [19],  as  well  as  recent  work  by  Clement  [17],  discusses  the  existence, 
boundedness  and  asymptotic  behaviour  as  t  ®  of  positive  solutions  of  (V) 
under  assumptions  which  are  also  motivated  by  the  heat  flow  problem  in  Chapter 
1.  This  problem  is  then  used  to  illustrate  the  theory.  The  reader  should 
recall  that  it  is  classical  that  solutions  of  the  heat  equation  are  positive, 
if  the  initial  temperature  distribution  and  the  external  heat  supply  are 
positive.  Some  of  the  results  of  this  chapter  were  presented  at  the  West 
Virginia  Conference  by  Clement. 
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NONLINEAR  VOLTE RRA  EQUATIONS  FOR  HEAT  FLOW  IN  MATERIALS  WITH  MEMORY 

John  A.  Nohel 


Chapter  1 

Volterra  Equations  Occurring  In  Heat  Flow  in  Materlala  with  Memory 
1.1.  Introductory  Remarks.  The  purpose  of  this  chapter  is  to  derive  frca  physical 
considerations  several  mathematical  models  for  nonlinear  heat  flow  in  materials  with 
memory.  For  simplicity  we  shall  limit  most  of  our  considerations  to  heat  flow  in  one  space 
dimension  and  remark  about  the  situation  in  the  multidimensional  ease.  The  primary 
objective  is  to  arrive  at  equation  (V)  in  Section  1.2  below  under  physically  reasonable 
assumptions  on  the  kernel  b,  the  operator  A  and  the  forcing  term  f.  This  particular 
mathematical  model  motivates  many  of  the  considerations  in  Chapters  2,  3.  and  4.  In 
Sections  1.3  and  1.4  we  shall  also  derive  two  other  models  for  heat  flow  in  materials  with 
memory  which  have  been  and  are  being  studied,  but  due  to  limitations  of  space  we  will  only 
refer  the  reader  to  the  relevant  mathematical  literature  for  their  analysis. 

The  mathematical  models  are  derived  from  the  following  general  considerations  of 
energy  balance  for  heat  transfer  in  a  body  B  in  if  (n  -  1,2,3).  If  C(t,x)  represents 
the  internal  energy,  q(t,x)  represents  the  heat  flux,  and  h(t,x)  represents  the  heat 
supply,  where  t  is  the  nine  and  x  is  the  position  in  the  body,  then  the  energy  balance 
equation  is 

(1«1)  ■  -divq  ♦  h  (t  >  0,  x  «  B)  . 

The  classical  linear  heat  equation  which  accurately  describes  heat  transfer  by  conduction 
in  uny  materials  is  derived  from  (1.1)  by  assuming  that  the  heat  flux  obeys  Fourier's  laws 

q  -  -cp  grad  u 

where  cfl  >  0  is  the  constant  thermal  conductivity  and  u  represents  the  temperature  in 
the  body  at  time  t  and  position  x.  It  is  also  assumed  that  the  internal  energy  depends 
linearly  on  the  temperature 


where  bg  >  0  is  the  constant  heat  capacity  and  cg  >  0  is  a  constant,  the  energy 
balance  (1.1)  then  yields  the  linear  heat  equation 

Sponsored  by  the  United  States  Army  under  Contract  Nos.  DAAG29-7 5-C-0024  and  DAAG29-B0- 
C-0041. 


which  adequately  describes  the  evolution  of  temperature  in  moat  homogeneous  and  isotropic 
rigid  bodies. 

However,  some  materials  exhibit  memory  effects  (materials  of  fading  memory  type,  see 
Coleman  and  Mizel  (231 )  for  which  the  classical  theory  is  unable  to  account.  Heat  flow  in 
such  material  is  modelled  by  assuming  that  the  internal  energy  c  and  the  heat  flux  q 
are  respectively  functionals  (rather  than  functions)  of  the  temperature  and  of  the  gradient 
of  temperature.  The  considerations  which  follow  are  based  on  extensive  research  by 
Coleman,  Curtin,  Noll,  Pipkin,  Mae  Canty,  Mizel,  and  Nunziato  (see  especially  Coleman  (211, 
Coleman  and  Gurtln  (22] ,  Coleman  and  Mizel  (23] ,  Gurtin  and  Pipkin  (39] ,  Mac  Camy  (57] , 

(581.  (59],  Nunziato  (68]). 

2.  A  Model  for  Nonlinear  Heat  Flow  In  a  Material  with  Memory.  We  consider  nonlinear  heat 
flew  in  a  homogeneous  bar  of  unit  length  of  material  with  memory  with  the  temperature 
u  •  u(t,x)  maintained  at  zero  at  x  »  0  and  x  »  1.  We  shall  assume  that  the  history 
of  u  is  prescribed  for  t  <  0  and  8  <  x  <  1.  The  equation  satisfied  by  u  in  such  a 
material  is  derived  from  the  assumptions  that  the  internal  energy  c  and  the  heat  flux 
q  are  functionals  (rather  than  functions)  of  u  and  of  the  gradient  of  u 
respectively.  According  to  the  theory  developed  by  Coleman,  Curtin,  Noll,  Pipkin,  Mac  Camy 
and  Nunziato  (see  e.g.  Mac  easy  (57],  (59]  and  Nunziato  [68])  for  heat  flow  lit  materials  of 
fading  memory  type  the  functionals  C  and  .  q  are  taken  respectively  as* 

(1*2)  s(t,x)  •  b  u(t,x)  ♦  /  $(t  -  s)u(s,x)ds  (t  >  0,  0  <  x  <  1)  , 

*  0 

t 

(1.1)  q(t,x)  -  -c#«(u|{t,x})  ♦  /  y(t  -  s)e(ux(s,x))ds  (t>0,  0  <  x  <  1)  • 

0 

e 

to  writing  the  functionals  c  and  q  we  have  assumed  for  simplicity  and  without  loss  of 


■F 


queer elity  that  Uw  history  of  the  temperature  u  is  prescribed  es  zero  for  t  <  0  (If 
this  is  sot  the  csss  end  if  the  history  of  u  is  sufficiently  smooth  for  t  <  0  and 
•  <  *  <  1,  this  has  the  effect  of  altering  the  forcing  term  h  in  equation  (1.4)  below  - 
and  consequently  also  C  in  (1.6)  below  -  by  additional  known  forcing  terms).  In  (1.2), 
(1.3)  bg  >  0,  Cg  >  0  are  given  constants,  t,Y  t  (0,**)  *  R  are  given  sufficiently  smooth 
functions  which  we  call  the  internal  energy  and  heat  flux  relaxation  functions 
respectively. 

flu  real  function  c  s  R  ♦  R  in  (1.3)  will  be  assumed  to  satisfy  the  assumptions 
(o)  of  Lemma  1.3  below.  It  should  be  noted  that  the  case  0(r)  £  r  gives  rise  to  the 
linear  model  derived  in  Nunziato  (68],  and  that  (1.3)  is  one  reasonable  generalization  of 
the  heat  flux  functional  for  nonlinear  heat  flow. 

In  the  physical  literature  it  is  customary  to  define 

«(t)  -  b  +  /  0(s)ds  (0  <  t  <  -) 

U  0 

as  the  internal  energy  relaxation  function;  thus  «(0)  -  b#  and  o’ (t) 

t 

K(t)  -  c  -  /  t(s)ds  (0  <  t  <  -) 

y  0 

is  defined  as  the  heat  flux  relaxation  function,  so  that  ie(0)  -  c_  and 

0 

then  (1.2),  (1.3)  are  replaced  respectively  by 

« 

0.3*)  c(t.x)  •  +  a(0)u(t,x)  ♦  /  o’(t  -  s)u(s,x)ds,  t  >  0  , 

0 

t 

0.3*)  q(t,x)  -  -x(0)a(u  )  -  /  *'(t  -  s)a(u  (s,x))ds  » 

*  0  * 

mo  (68,  (5.13),  (5.11)]  where  the  linear  ease  o(r)  I  r  la  considered.  The  quantity 
•(0)  is  called  the  instantaneous  heat  capacity  while  o(«)  is  the  equilibria  heat 


■  8<t).  Similarly, 


x*(t)  -  -T(t). 


-3- 


capacity  and  similar  definitions  for  the  k(0 )  and  ■((•>).  Xt  Is  shown  in  [66]  that 
«(•)  >  a(0)  >  0,  and  that  k|0),  *(«•)  are  nonnegative. 

In  the  physical  literature  6,  Y  are  usually  assumed  to  be  decaying  exponentials  with 
positive  coefficients.  As  we  shall  see  the  theory  developed  in  Chapters  2-4  permits  a  much 
greater  generality,  and  we  shall  merely  have  to  require  that  6(0)  >  0,  Y CO )  >  0,  that  6 
and  y  *  l»* (0,“),  and  that 


t 

b ♦  /  $<T)dT  >  0, 

"  it 


t 

e  -  /  Y(tJ6x  >  6 
8  0 


(0  <  t  <  -)  , 


which  corresponds  to  the  physically  reasonable  assumptions  a(t)  >  0,  r(t)  >0,  0  <  t  <  •• 
MS  shall  also  assume  that  the  conditions 

<pw)  b0>  Re  6(in)  >0  (n  «  R>  . 

m 

ft)  cfl  -  /  Y(t)dT  >  0  , 

0 

*  • 

where  8(iu)  •  /  6(t)exp(-int)dt,  are  satisfied;  assumption  (y)  states  that  <(•)  >  0. 

0 

the  above  assumptions  will  be  motivated  presently.  Remark  4.8  in  Chapter  3  below  shows 

t  t 

that  the  physically  reasonable  assumptions  b p  *  f  6(T)dT  >  0  and  c q  ”  /  Y(*)dT  >  0 

0  0 

(0  <  t  <  “)  are  actually  not  essential  for  the  theory  developed  in  Chapter  3  to  apply. 
Xt  b  •  h(t,x)  c  L?  (0,“;L2(0,1))  represents  the  external  heat  supply  added  to  the 

IOC 

rod  for  t  >  0  and  0  <  x  <  t,  and  if  u(0,x)  -  Ug(x),  0  <  x  <  1,  is  the  given  initial 
temperature  distribution,  the  law  of  balance  of  heat  (1.1)  shows  that  in  one  space 
dimension  the  temperature  u  satisfies  the  initial-boundary  value  problem 


Jb0u  6*u]  “  e06<ttx>x  “  Y*«(«x)x  ♦  b  (0  <  t  <  •,  0  <  x  <  1) 
u(0,x)  -  uQ(x)  (0  <  x  <  1),  u(t,0)  -  u(t,1)  5  0  (t  >  0)  , 


where  subscripts  denote  differentiation  with  respect  to  x  end  where  *  denotes  the 
convolution  on  (0 ,t] .  Note  that  in  an  ordinary  material  6  ”  Y  5  0,  end  (1.4)  becomes 
the  nonlinear  heat  equation  in  one  space  dimension. 


The  next  task  is  to  transform  (1.4)  to  the  equivalent  form  (V)  below  which  will  be 
used  for  the  analysis.  Define 

t 

O.S)  C(t)  »  e  •  /  y(t)dx 

0 

t 

(1.*)  G(tfx)  “  b  u  (x)  +  /  h(x,x)dt 

0 

Noting  that 

(C*o(ux)x)(t,x)  -  e0o(ux(t,x))x  -  (y*o(ux)x)(t,x)  , 
and  integrating  (1.4)  using  the  initial  condition,  and  (1.6)  yields  the  equivalent  Volterra 
equation  (to  (4.3)): 

(1.7)  D0u(t#x)  ♦  (S*u)(t,x)  -  (C*o(ux)x)(t,x)  ♦  C(t,x)  (0  <  t  <  “,  0  <  x  <  1)  , 

where  u  satisfies  the  boundary  conditions  u(t,0)  -  u(t,1)  S  0  (t  >  0). 

Wa  next  define  the  nonlinear  operator  A  formally  by  the  relation 
(*)  An  “  -o(u  J  where  uft.O)  •  u(t,1)  SO. 

Xn  order  not  to  interrupt  this  development  we  postpone  a  precise  definition  of  A  to  Lemma 
1.3  below.  Then  the  Volterra  equation  (1.7)  has  the  abstract  form 
(V|)  bflu  ♦  #«u  ♦  C*Au  »  G  (0  <  t  <  «)  . 

To  transform  (Vt)  to  the  equivalent  form  (V)  below  define  p  :  (0,»)  *1  to  be  the  unique 

solution  of  the  linear  Volterra  equation  (called  the  resolvent  kernel  of  6)* 

<P>  b  p(t>  ♦  (d*p)(t)  -  -  (0  <  t  <  -)  . 

*  *9 

Xt  la  standard  that  if  >  0  and  $  e  equation  (p)  has  a  unique  solution 

9  •  (0.-).  Applying  the  variation  of  constants  formula  for  Volterra  equations  (63) 

<*0*  *  t*r  •  9  <”">  y  «  ♦  p*g) 

finally  yields  that  (V^)  is  equivalent  to  the  abstract  equation 
(»)  U  ♦  b*Au  ■  f  (0  <  t  <  -) 


(0  <  t  <  •)  . 

(0  <  t  <  •»  0  <  x  <  1)  . 


with  the  definitions 


(0  <  t  <  -> 


CM  bCt)  -  ♦  (p*c)(t) 

bo 

<n  m>  -  +  (p*G)(t,«)  to  <«<•*)  . 

bo 

Similar  considerations  show  that  for  heat  flow  in  a  bounded  homogeneous  body  Q  of 
isotropic  material  with  memory  in  R2  or  R2  with  a  smooth  boundary  T,  the 
temperature  u  will  also  satisfy  the  abstract  Volterra  equation  (V)  with  the  kernel  b 
and  forcing  term  f  given  exactly  as  above,  but  with  the  nonlinear  operator  A  defined 
precisely  in  Remark  1.4. 

He  next  comment  on  the  significance  of  the  assumptions  concerning  0,  y  as  well  as 

(PH)  and  (y).  Since  the  relaxation  functions  0  and  Y  are  generally  taken  as  decaying 

exponentials  with  positive  coefficients  in  the  physical  literature,  it  is  certainly 

reasonable  to  assume  that  0,Y  «  L1  {<),*•)  and  that  0(0)  >  0,  Y(0)  >  0.  He  next  motivate 

t 

the  assumption  that  b  ♦  /  0(T)dr  >0  (0  <  t  <  <■).  A  similar  reasoning  motivate* 

t  o 

c ■  -  /  Y(T)dt  >0  (0  <  t  <  •).  Consider  the  internal  energy  e  defined  by  (1.2)  and 

0 

suppose  that  the  temperature  u  is  maintained  at  zero  up  to  time  t„  and  at  a  state 

Uj(x)  >  0  (0  <  x  .<  1)  for  t  >  tg.  (Ate  would  then  expect  the  internal  energy  to  be 

positive  for  t  >  t  .  Zf  the  function  0  is  positive  for  t  >  0  this  is  automatically 

t 

the  case.  However  if  not,  the  assumption  b_  ♦  /  0(T)dr  >0  (0  <  t  <  ")  is  natural  in 

B  0 

view  of  the  fact  that  in  this  situation 

_  t 

«(t,X)  ■  u  (x)(b  +  /  0(T)dT)  (t  <  t  <  •)  . 

tQ 

Since  0  c  I.1  (0,«)  equation  (4.1)  shows  that  c  is  bounded  whenever  u  is 

m 

bounded.  The  assumption  (PH)  implies  that  b_  ♦  /  0(t)dt  >  0  (taka  n  ■  0)t  thus  it 

_  0 

u(x,t)  tends  to  sn  equilibrium  state  u(x)  >0  as  t  ♦  •,  (1.2)  implies  that  the 

corresponding  limiting  internal  energy  7(x)  >0  as  is  to  be  expected.  For  physical 
reasons  it  is  also  to  be  expected  that  if  c  is  bounded  the  temperature  should  be 
bounded.  Applying  the  variations  of  constants  formula  to  (1.2)  yields 


u(t,x)  -  +  (p*e)  (t,x)  (0  <  t  <  -.  0  <  x  <  1)  . 

b0 

where  p  is  the  resolvent  kernel  of  8  defined  by  equation  {p>.  Thus  to  have  u  bounde 
whenever  e  is  bounded  it  is  sufficient  to  require  that  P  e  l'(0,**).  But  by  the  Paley- 
Wiener  theorem  [69]  applied  to  equation  (p),  8  t  l\o,»)  implies  that  p  f  l'(0,®)  if 

and  only  if 

bg  +  B(z)  *  0  for  Re  z  >  0  • 

The  condition  (PW)  now  results  from  taking  the  real  part  of  this  expression,  noting  that 
for  physical  reasons  one  wants  b^  +  8(0)  >  0,  and  arguing  as  in  the  proof  of  Lemma  2.2, 
Chapter  3. 

To  motivate  assumption  (y)  suppose  that  u(t,x)  *  u(x)  as  t  *  •  and  that 

u(x)  >  0,  implying  that  c(^)  >  0  (see  assumptions  (o)  below).  One  then  expects 
ax  ax 

that  the  limiting  heat  flux  q(x)  in  equation  (1.3)  is  strictly  negative,  if  the  process 
being  modelled  represents  "forward”  heat  flow;  condition  (y)  insures  that  this  is  the 
case. 

We  shall  next  see  that  the  physically  reasonable  assumptions  bQ  >  0,  cQ  >  0,  (PW) 
and  (y),  together  with  some  mild  additional  technical  assumptions,  imply  that  the  kernel 
b  in  the  Volterra  equation  (V)  defined  (b)  satisfies  the  assumptions  - 
(Hb)  b(t;  -  b^  +  B(t),  b(0)  >  0,  b„  >  0,  B,B'  C  lVo,-)  ; 

these  will  play  an  important  role  in  the  application  of  the  boundedness  and  asymptotic 
theory  developed  in  Chapter  3,  Section  3,  and  in  the  application  of  that  theory  to  heat 
flow  described  by  the  problem  (1.4)  above  in  Chapter  3,  Section  4.  One  has  the  following 
result  whose  elementary  proof  is  omitted: 

Lemma  1.1.  Let  bQ  >  0  cQ  >  0,  8,  y,  tB,  ty  «  L1  (<>,•),  and  let  assumptions  (PW) 
and  (y)  be  satisfied.  Define 


c0  -  /  T(t)dt 
0  0 _ 

m 

b.  ♦  /  8  ( t )  dt 

"  A 


(1.8) 


(1.9) 


B(t)  -  *  <p* C)(t)  -  , 

t 

where  C(t)  =  c  -  /  y(T)dT,  and  p  is  the  resolvent  of  6  uniquely  defined  by 
0  . 

equation  (p).  Then  b  >  0  and  B,B'  £  L  (0,“),  and  b(t)  -  b^  +  B(t)  satisfies  (H^) 
c0  co 

with  b(0 )  =  —  >  0,  B(0)  =  -  b  >0. 

-  b0  b0 

The  next  elementary  result  gives  physically  reasonable  sufficient  conditions  on  the 
relaxation  function  8#  the  initial  temperature  distribution  Uq,  and  the  external  heat 
supply  h  in  order  that  the  forcing  term  f  in  (V)  defined  by  equation  (f)  will  satisfy 
the  assumption 

(H.)  f(t)  «  f  +  F(t),  F  £  W,1,2([0,“);H),  F*  e  L2(0,«»H)  , 

r  oo  loc 

where  H  is  the  real  Hilbert  space  L2<0,1)  and  W  is  the  usual  Sobolev  space. 
Assumption  (H ^ )  will  play  an  important  role  in  the  boundedness  and  asymptotic  theory  of 
Chapter  3. 

Lemma  1.2.  Let  H  -  L2(l,1)  and  uQ  E  hJ(0,1).  Let  B  e  l’iO,-)  n  l2(0,«)  and  let 

assumption  (PW)  be  satisfied.  Finally,  assume  that 

(h)  h  e  L1(0,«;R)  ^  L2(0, >•,•«)  . 

Then  the  function  f  :  [0,“)  x  (0,1)  ♦  H,  defined  by  equations  (f),  (G),  where  p  is  the 
resolvent  of  B,  satisfies  f  e  W^2(0,»,-H)  and  f(0,x)  =  uQ(x)  e  h’(0,1).  Moreover, 
f(t,x)  -  fjx)  +  F(t,x)  (0  <  t  <  0  <  x  <  1)  , 

where 

m  m 

(1.10)  f„(x)  -  (b  u  (x)  +  /  h(T,x)dr)  (—  ♦  /  p(T)dt)  , 

'  0  0  0 


(1.11) 


and 


JF 

at 


F(t,x)  -  G(-gfX)  +  (p*G)(t,x)  -  f„(x)  =  -  r-  J  h(T,x)dT 
“o  0  t 

t  *>  "  • 

-  /  p(t-s)  /  h(T,x)dTds  -  f  p(T)dr(b  u  (x)  +  /  h(T,x)dt)  , 

0  s  t  0 

E  L2(0,®;H).  If  in  addition  tS  f  I.1  ( 0 ,—)  and  th  e  L^O.-jH),  then  F  t  L2(0,®;H). 


Sketch  of  Proof  of  Lemma  1.;.  The  assumptions  g  e  L  (0,*»)  and  (PW),  together  with 

the  Paley-Wiener  theorem  [69],  applied  to  the  resolvent  equation  (p)  imply  that 

p  f  i/fO,**).  But  then  the  assumption  g  e  L2(0,»)  and  the  fact  that  p  t  L^O,**)  imply 
2 

that  also  p  e  L  (0,“)  from  the  resolvent  equation.  These  facts  combined  with  the 

definition  of  f  in  (f)  and  assumption  (h)  yield  the  formulae  (1.10)  for  f^  and  (1.11) 

for  F  given  in  the  statement,  as  well  as  f  t  W?'2(0,»>H).  From  formula  (1.11)  one 

loc 

easily  proves  that 

3F  1 

0*12)  3^  (t.x)  -  —  h(t,x)  +  bQuQ(x)p(t)  +  (p*h)(t,x)  (0  <  t  <  •>,  0  <  x  <  1 )  » 

°0 

then  t  l2(0,“»H)  follows  from  h  t  L2(0,"iH)  and  p  t  L^O,*)  n  l2(0,«).  Finally, 

(PW)  and  t6  e  l'io,"),  together  with  p  c  L^O,-)  imply  that  tp  *  l’io,*)  from  the 

resolvent  equation.  This,  together  with  the  assumption  the  L^O,**)  and  routine 

2 

estimates  applied  to  the  formula  (1.11)  yield  Fe  L  (0,»jH).  This  completes  the  proof. 

The  next  task  is  to  jive  a  precise  definition  of  the  operator  A  in  the  abstract 
equation  (V)  for  the  heat  flow  problem  under  study.  Let  H  -  L2(0,1)  be  the  real  Hilbert 
space  of  square  integrable  functions  on  (0,1).  Let  c  :  R  ♦  R  satisfy  the  assumptions 
(o)  a  t  c’dO,  o(0)  ■  0,  o’  (?)  >  p0  >  0  (5  *  R)  , 

for  some  pQ  >  0.  Define  W  s  R  ♦  R+  by 

r  p.  , 

W(r)  -  /  o(C)dC  (>  —  r  )  (r  «  R)  , 

0  1 


and  define  •fi  :  H  +  (-«•,+"]  by 


(1.13) 


*(y) 


if 


t  *  h’{0,1) 


otherwise  . 


Lemma  1.3.  Let  the  assumptions  (o)  be  satisfied  and  let  •fi  i  H  ♦  (-»,♦«•]  be  the  function 
defined  by  (1.13).  Then  *'  is  convex,  l.s.c.  and  proper  on  H,  and 


-9- 


Ay  -  3*<y)  -  -  ^  D(3'f)  "  {v  *  h’(0,1)  :  o(|^)  c  L2<0,1)}  , 

where  3^  denotes  the  subdifferential.  Moreover,  <?(y)  >  0,  *(y)  ♦  •  as 

1  2  '  2 
/  |y(x) |Z dx  ♦  ",  and  (y,3^(y))  >  p.n  /  |y(x)|  dx,  where  (*,♦)  denotes  the  scalar 

0  0 
product  In  H. 

Sketch  of  Proof  of  Lemma  1.3.  The  first  result  is  standard  see  Brezis  (114], (15)). 
To  prove  the  last  two  conclusions  let  y  e  H^(0,1)f  then  from  the  definition  of  and 
the  Poincare  inequality  one  has 


•fly )  >  ^  /  l|^  <x)|2dx  >  v2  I  |y(x)|2dx  >  0  , 

0  0 


and  <e  (y)  *  "  as  /  |y(x)|  dx  ♦  Moreover, 

0 


(y 


,3v>(y))  -  -  /  y(x)  { x ) ) dx  , 


and  an  integration  by  parts,  y  t  Hfl (0, 1 )  and  the  Poincare  inequality  give 


(y.**(y ))  -  /  (x)o(^  (x))dx  >  P0  /  (X)|2dx  >  p  X2  /  |y(x)|2dx. 


„  dx  '"''•'dx  - ‘'O  ‘  ’dx  - '0"  . 

0  0  0 


While  our  considerations  of  the  heat  flow  problem  are  primarily  in  one  space 
dimension,  we  indicate  how  to  define  the  nonlinear  operator  A  •  3*  in  (V)  for  the  heat 
flow  problem  in  two  or  three  space  dimensions  such  that  the  function  v>  satisfies  the 
conclusions  of  Lemma  1.3  under  the  physically  reasonable  assumption  (1)  below. 

Remark  1.4.  Let  0  be  a  bounded  domain  in  Rn  (for  heat  flow  n  •  2  or  3)  with  smooth 
boundary  T.  Let  X  :  R+  ♦  R  be  a  given  smooth  function  satisfying  the  assumption 


(A) 


{A ( 0 )  >  0,  there  exists  pQ  >  0  such  that  1(5)  *  pQ  and 
5A’(C>  +  A<$)  >  pQ  l(  (  «|  . 


Define  A  :  R 


R 


by 


Y.' 
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r  p.  , 

A(r)  -  /  (X(()d(  (>  r*  r  )  (r  «  *)  . 

•  * 


Let  H  ■  l2{fl)  end  define 


1/  A( | Vu| )dx 
0 


if  u  e  Hg (0) 
otherwise  . 


Then  it  is  readily  verified  that  i  H  ♦  (-■,+■]  is  convex,  1 .s.c. ,  end  proper  on  H  end 

An  -  9«u)  -  -V*  (X(  I  Vo|  )Vu) 

with 

D(J^)  -  {u  e  H^(0)  t  V*(XC|7u|)7u)  *  l*(Q )}  . 

Cleerly  *(u)  >0  (u  <  H)  end  by  the  Poincare  inequality  *»(u)  ♦  "  as  |u[  ♦  •.  Osinp 
Integration  by  parts  end  the  Poincare  inequality  one  also  has  the  coercivity  condition 

(ta,n)  >  kp0Iu|*  , 

where  k  >  0  is  the  constant  in  the  Poincare  inequality. 

The  results  of  Lesraas  1.1-1. 3  and  of  Remark  1.4  will  be  used  in  Chapter  3,  Section  4, 
to  discuss  the  global  existence,  uniqueness,  boundedness  and  asymptotic  behavior  of  the 
solution  of  the  initial-boundary  value  problem  (1.4)  using  the  theory  which  will  be 
developed  for  the  equivalent  abstract  Vol terra  equation  (V). 


3.  h  General  (Parabolic)  Heat  Flow  Model.  He  consider  the  same  heat  flow  problem  as  in 
Section  2.  In  this  model  we  assume  that  the  internal  energy  is  given  by  the  functional 
(1.2),  but  we  assume  the  following  more  general  form  of  the  heat  flux  functional 

(1.14)  «(t,x)  -  -f(u  )  -  /  a(t  -  s)o(u  (s,x))ds  (t  >  0)  , 

*  j  * 

where  we  again  assume  that  the  history  of  temperature  is  prescribed  as  sero  for  t  <  0 
and  0  <  x  <  1  •  The  real  function  7  i  It  ♦  X,  1(0 )  »  0  satisfies  the  same  assumptions  as 


0  la  Lei 


1.3,  Section  2.  Xf  ♦(•)  «  c.o(«),  and  a  •  -r  (l.a.  in  the  notation  of 


(t.l')  fc'(t)  «  a(t))  (1.14)-  reduces  to  (1.3).  We  ittwe  that  $,  a,  and  the  external 
heat  supply  h  satisfy  the  same  type  of  smoothness  assumptions  respectively  as  8  and  y 
in  Section  2.  We  assume  again  that 


b-  ♦  /  8(T)dt  >0  <0  <  t  <  •>  . 


and  that  a(0)  >  0.  For  physical  reasons  it  is  also  reasonable  to  assume  that 


♦«)  ♦  (/  a(T)dt)o(5) 

0 


0  if  C  >  0#  t  >  0 
0  if  C  <  0,  t  >  0 


and  that 


♦(C)  ♦  (/  a(t)dt)oCC) 

0 


0  If  l  >  0 
0  if  C  <  0 


applying  the  energy  balance  equation  (1.1).  and  using  (1.2),  (1.14)  and  h  shows  that 
the  temperature  u  satisfies  the  initial-boundary  value  problem 

b0«t  +  8(0)u  ♦  8**u  -  ♦  a*o(«  )x  ♦  h  (0  <  t  <  •»,  0  <  x  <  1) 

(1.15)  u(t,0)  -  u(t,1)  SO  (t  >  0) 

u(0,x)  -  »0<*)  (0  <  x  <  1)  . 

Bote  that  if  a  »  8  *  0,  (1.15)  reduces  to  the  nonlinear  heat  equation.  Defining  the 
epertors  A,  B,  L  by  the  relations: 

An  »  where  n(t,0)  ■  u(t,1)  s  0 

Bu  ■  where  u(t,0)  ■  u(t,1)  s  0 

In  “  8(0)u  ♦  8**u  , 

and  taking  (without  loss  of  generality)  bg  »  1,  one  sees  that  the  problem  (1.15)  has  the 
Abstract  form 


jj*  ♦  Bn  ♦  a*Au  ♦  tu  -  h  (0  <  t  <  •) 


•(•»  "  *o 


Xf  the  functions  o  and  f  satisfy  tha  assumptions  (o).  Lemma  1.3  shows  that  A  and  B 
are  subdifferentials  of  proper  convex,  l.s.c.  functions  defined  on  the  real  Hilbert  space 
■  “  L2(0,1).  If  the  relaxation  function  S  has  ft*  <  the  linear  operator  L  is 

troll  defined. 

The  abstract  problem  (V* )  has  been  investigated  by  a  number  of  authors  combining 
techniques  of  Volterra  equations  and  the  theory  of  monotone  operators.  If  0  =  0  v.  Barbu 
(5] ,  (7],  Barbu  and  Malik  (11]  studied  the  problem  of  global  existence;  a  store  complete 
existence  theory  for  considerably  more  general  kernels  a,  as  veil  as  a  discussion  of 
boundedness  and  asymptotic  behaviour,  was  developed  by  Crandall,  London,  and  HOhel  [28]. 

The  latter,  see  (28,  p.  717],  also  permits  in  the  existence  theory  Lipschits  type 
perturbations  of  the  operator  Aj  this  essentially  covers  the  case  of  the  specific 
operator  L  in  the  present  application;  the  details  of  this  generalization  of  the 
existence  theory  have  recently  been  vorked  out  by  M.  J.  Luo  as  a  part  of  a  forthcoming 
Ph.D.  thesis  at  the  University  of  Wisconsin-Madison. 

The  general  assumptions  for  the  existence  theory  in  [28]  concerning  the  operators 
A  and  B  are  roughly  speaking  as  follows.  The  operators  A  and  B  are  subdifferentials 
of  proper,  convex,  l.s.c.  functions  defined  on  a  real  Hilbert  space  H;  B  dominates  A  in 
a  certain  precise  sense  (see  (28],  Inequality  (1.7);  the  case  B  ■  kA,  k  >  0  is  not 
excluded),  and  B  satisfies  a  compactness  assumption  (see  (1.8)  in  [28];  this  compactness 
condition  excludes  the  possibility  B  5  0.  The  kernel  a  satisfies  an  abstract  condition 
(see  (28;  (1.10)])  which  is  shown  to  hold  for  two  physically  important  classes  of  kernels; 

|a(0)  >0,  a  locally  absolutely  continuous  on  (0,*>)  , 
a*  locally  of  bounded  variation  on  (0,“)  ; 


(na> 


a(0)  >  0,  a  c  c[0 »•)  n  c2(0,«»),  and 
a  nonnegative,  nonincreasing,  convex  on  (0,*) 


St  should  be  noted  that  no  satisfactory  uniqueness  theorem  has  to  date  been  discovered 
for  the  general  problem  (V£)  with  B  *  kA,  k  >  0  a  constant,  even  when  L  S  0* 

Other  approaches  to  the  study  of  (V£)  include  Aizicovici  (11,  (2),  (3]  and  Gripenberg 
(32],  [36].  While  these  studies  are  in  a  more  general  setting  than  those  referred  to 
above,  they  do  not  seem  to  shed  new  light  on  the  physical  problem  (1.15),  and  indeed  the 
existence  results  of  Aizicovici  do  not  include  those  of  (28) .  An  extension  of  the  latter 
to  nonconvolution  kernels  has  been  studied  by  Rennolet  [71].  A  semigroup  approach  to  a 
special  case  of  (V*)  has  been  investigated  by  Vrabie  (83).  Another  interesting  "parabolic" 

Is 

heat  flow  problem  involving  an  analysis  of  evolution  inequalities  has  recently  been  studied 
by  J.  Haumann  (64). 

4.  A  Hyperbolic  nonlinear  Volterra  Equation  for  Heat  Conduction  with  Finite  Wave  Speeds. 
The  parabolic  models  for  heat  flow  in  materials  with  memory  formulated  in  Sections  2  and  3 
both  predict  that  a  thermal  disturbance  at  any  point  of  the  body  is  Instantly  felt 
everywhere  in  the  body  (though  not  with  equal  strength).  This  implies  that  finite 
discontinuities  propagate  with  Infinite  speed.  This  situation  is  unrealisitie  for  some 
materials,  particularly  at  low  temperatures. 

Curtin  and  pipkin  [39],  see  also  thiaziato  [68],  have  proposed  a  model  for  heat  flow 
which  exhibits  a  finite  speed  of  propagation)  they  study  of  the  linear  model  0(r)  3  r. 

We  present  briefly  «uch  a  nonlinear  model  which  has  been  investigated  by  Mac  Camy  [59]  by 
the  method  of  characteristics,  and  by  Dafemos  and  Nohel  [29]  by  an  energy  method.  Another 
interesting  variant  of  the  energy  method  was  recently  developed  by  Stef fans  [80]. 

Consider  the  heat  flow  problem  of  Section  2.  Define  the  internal  energy  by  the 
functional  (1.2),  with  8  satisfying  the  conditions  of  Section  2.  In  place  of  (1.3) 
assume  that  the  heat  flux  is  given  by  the  functional 

« 

(1*16)  q{t,x)  -  -  /  c(t  -  s]e(u  (s,x))ds  (0  <  t  <  •,  0  <  x  <  1)  , 

9  * 


1.3 


where  the  relaxation  function  c  «  L  (0,<*)#  c(0)  >  0,  f  c(s)ds  >  0  (0  <  t  <  •), 

•  0 

/  c(a)ds  >  0,  and  where  the  real  function  a  satisfiee  assumptions  (o)  of  Lemma 

• 

(actually  o,  smooth,  o(0)  ■  0,  o'(0)  >  0  la  sufficient  for  the  development  in  (29]). 

Rote  that  from  (1.16)  the  heat  flux  depends  only,  on  the  history  of  the  gradient  of  u,  and 
is  independent  of  the  present  value  of  the  gradient  of  u.  Evidently,  this  model  of  the 
heat  flux  results  from  (1.3)  by  talcing  c„  «  0  and  by  replacing  Y  by  -c»  or  from 
(1.3')  by  taking  k(0)  -  0,  x*(t)  •  c(t).  The  model  (1.16)  for  the  heat  flux  also  results 
by  taking  p  3  0  in  (1.14),  a  case  which  is  excluded  in  the  theory  and  the  referenced 
accompanying  mathematical  literature  in  Section  1.3. 

Applying  the  energy  balance  (1.1)  and  using  (1.2),  (1.16)  and  the  external  heat 
supply  h  leads  to  the  equation 

b0ut  ♦  6*u  -  0*0(0^  ♦  h  (0  <  t  <  •,  0  <  x  <1)  . 

noting  that  (3*u  -  u*S.  carrying  out  the  differentiation,  and  imposing  the  boundary  and 
initial  conditions  leads  to  the  following  initial-boundary  value  problem  for  the  heat  flow 
problem  under  the  assumptions  of  this  section 

bflut  ♦  e*ufc  *  0*0(0^  ♦  h(t,x)  -  0(t)ug(x)  (0  <  t  <  ■»,  0  <  x  <  1)  , 

(1.17)  u(t,0)  -  u(t.l)  SO  (t  >  0) 

,  u(0,x)  -  ufl(x)  . 

The  problem  (1.17)  is  transformed  to  the  more  standard  abstract  form  (V')  below  as 

A 

follows.  As  in  Section  2  define  the  resolvent  kernel  p  of  6  by  the  linear  Volte rra 

equation  (p)t  by  assumption  (PH)  p  <  1/(0, •).  Definu  the  function  a  1  (0,**)  ♦  R  by 

a(t)  -  c(t)  ♦  (p*c)(t)  (0  6  t  <  •)  1 

b0 

dafine  the  function  9  i  (0,«)  »  (0,1)  ♦  I  by 

g(t,x)  -  ^  (h(t.xj  -  S(t)u^(x))  ♦  p*(h(t,x)  -  0(t)«^(x))  (0  <  t  <  -,0  <  x  <  1)  . 

finally  dafine  the  operator  A  as  in  Section  2.  Applying  the  variation  of  constants 
formula  for  Volterra  equations,  and  these  definitions  to  (1.17)  shows  that  (1.17)  la 
equivalent  to  the  abstract  nonlinear  Volterra  equation 


to  <  t  <  -) 


*  a*Au  -  g(t,») 

U(O)  ■  B((»)  . 


Xfc  is  easy  to  in  (compare  Learns  1.1)  thst  under  the  present  assumptions  on  0  and  c, 
a  *  l\o,-)  and  that 


/  a(t)dt 


/  c(t)dt 
0 _ 


>  0  i 


b ^  /  B(t)dt 


If  also  0(0)  >  0.  e(0)  >  0,  e'(0)  <  0,  then  a(0)  >  0  and  a*(0)  <  0.  It  is  also 


evident  that  the  forcing  tern  g  «  L2(0,«»La(0,1))  if  he  L2(0,«>|L3(0,1) )  and  also 


M  La<0,-). 


Notice  that  if  0  3  0  and  f  3  y(0)  >0  (or  equivalently  a  3  ■■)  the  problem 

b0 


(1.1?)  reduces  to  the  nonlinear  wave  equation  problem 


00 


Vtt  *  rt0>  "Vx  *  \ 


(0  «  t  <  •,  0  <  *  <  1) 


e(t,0)  -  u(t,1)  30  (t  >  0) 


u(0,x)  -  uQ{x),  u^(0,x)  -  u^(x)  »  g(0,x)  . 


If  the  real  function  o  is  "genuinely  nonlinear"  (o"(C)  |  0,  (  •  K).  Lax  (SO]  has  shown 
that  (W)  falls  to  have  global  smooth  (C2)  solutions  in  time,  even  if  h  3  0,  no  matter  how 
smooth  and  "small"  one  takes  the  initial  functions  uQ,  u,.  Zf  the  function  o  is  convex 
the  derivatives  of  the  solution  u  of  (W)  develop  singularities  due  to  the  crossing  of 
characteristics  in  finite  time  ("shocks"). 

The  objective  of  the  analysis  by  Kac  Camy  (59]  (which  uses  Niemann  invariants  and  is 
therefore  restricted  to  one  space  dimension),  and  a  different  analysis  by  Dafermos  and 
•abel  (29,  Theorem  4.1],  Staffans  [80],  which  are  applicable  to  several  space  dimensions 
and  both  of  which  use  energy  methods,  is  to  show  that  under  the  present  assumptions  on 
0  and  T  (together  with  some  technical  ones  and  some  other  physically 


reasonable  ones  which  imply  that  the  kernel  a  Is  strongly  positive  on  [0  ,**))«  the 
Integral  in  (V*)  has  the  effect  of  a  frictional  damping  mechanism  which  prevents  the 
formation  of  shocks,  provided  the  data  Ug  and  the  forcing  term  g  are  sufficiently 
smooth  and  small  in  certain  H  norms.  This  analysis  leads  to  global  existence, 
uniqueness ,  and  decay  of  smooth  solutions  of  the  problem  (1.17)  for  sufficiently  smooth  and 
small  data  Ug  and  h  (see  especially  [29,  Theorems  4.t,  6.1);  for  a  physically 
reasonable  two-dimensional  heat  flow  problem  with  the  same  kernel  a  see  (29,  Theorem 
7*1].  It  is  also  evident  from  the  analysis  in  [29,  Section  3]  dealing  with  the  local 
existence  and  uniqueness  of  the  problem  (1.17),  resp.  (V^) ,  that  solutions  of  (1.17), 
reap.  (V'),  possess  the  property  of  finite  speed  of  propagation.  For  an  analysis  of  the 

a 

existence  and  uniqueness  of  classical  solutions  for  "small"  data  in  the  simpler  case  of  a 
nonlinear  wave  equation  with  frictional  damping  see  Nishida  (65) ,  and  Nohel  [66] . 

We  remark  also  that  the  abstract  problem  (V*)  has  recently  been  studied  by  S.  O. 

H 

Londen  [55],  [56]  for  a  class  of  kernels  a  which  are  positive,  decreasing,  and  convex  on 
(0,*)  and  which  satisfy  the  crucial  (for  his  method)  condition  a'(0+)  «  His  method 

is  a  significant  generalization  of  that  of  Crandall,  Londen,  and  Nohel  [28]  for  the 
parabolic  problems  discussed  briefly  in  Section  3.  However,  the  assumption  a'  (0+)  - 
is  inappropriate  in  the  present  physical  context;  moreover,  the  type  of  solution  obtained 
by  Londen  in  [55]  and  [56]  need  not  be  regular  in  the  sense  of  smooth  solutions,  and  no 
decay  results  of  solutions  comparable  to  [29]  are  obtainable  by  his  methods.  It  should 
also  be  remarked  that  for  the  linear  problem  (V^)  (e.g.  Au  •  -V2 u) ,  interesting  and 
useful  results  using  deep  results  of  Shea  and  wainger  [73]  have  been  obtained  in  a  series 
of  papers  by  K.  S.  Hannsgen  [41-47]  and  by  Carr  and  Hannsgen  [16]. 

Finally,  we  note  the  model  problem  (▼')  is  similar  to  a  particular  model  problem  for 
nonlinear  viscoelastic  motion  in  one  space  dimension  in  which,  however,  the  kernel  a  has 
the  form  a(t)  m  *m  +  Mt),  as  >  I,  A  positive,  decreasing  convex  on  (0,*),  and  for 
which  the  analysis  is  considerably  more  complicated.  This  problem  has  been  extensively 
studied  by  Hae  Cany  [60] ,  Defence  and  Nohel  [29] ,  and  Staffans  [SO]  t  see  also  an 
interesting  asymptotic  result  by  Staffans  [81]  motivated  by  this  problem. 


Existence  an<i  Uniqueness  of  Solutions  of  Abstract  Volterra  Equation* 


2.1.  Introduction.  In  this  chapter  we  study  the  abstract  nonlinear  Volterra  equation 
(V)  u(t)  ♦  b*Au(  t)  a  f(t)  (0  <  t  <  T)  , 

where  T  >  0  is  arbitrary,  in  the  setting:  A  is  an  m-accretive  (possibly  multivalued) 

operator  in  a  real  Banach  space  X,  the  given  kernel  b  is  a  real  absolutely  continuous 

t 

function  on  (0,T] ,  b*g(t)  “  /  b(t  -  s)g(s)ds  with  the  integral  in  (V)  interpreted  as  the 

0 

usual  Bochner  integral,  and  the  given  function  f  e  W1#1(0,*|X)  where  W1'1  is  the  usual 
Sobolev  space. 

We  treat  the  problem  of  existence,  uniqueness,  dependence  on  data,  and  regularity  of 
solutions  of  (V)  on  [0,1]  by  means  of  a  simple  method  developed  jointly  with  M.  6. 
Crandall  (26]  to  which  the  reader  is  referred  for  more  details*  the  results  obtained  for 
(V)  generalize  and  simplify  considerably  earlier  work  on  existence  and  uniqueness  obtained 
by  Barbu  [6],  (8],  London  [54],  Gripenberg  (30),  for  the  case  X  -  B  a  real  Hilbert 
Space.  A  different  approach  to  the  study  of  (V)  in  the  same  general  setting  was  developed 
Independently  by  Gripenberg  (31],  (32].  The  general  theory  will  be  used  in  Chapters  3  and 
4.  He  will  also  comment  on  the  special  cases:  (i)  A  maximal  monotone  on  H,  and 
(11)  A  -  3 <f,  where  c  i  d  ♦  (•*,+*]  is  a  proper,  convex,  l.s.c.  function  and  3 <f  denotes 
the  subdifferential  of  *  (see  Brezis  (14]);  these  special  cases  will  be  important  in 
Qiapter  3.  He  shall  also  consider  briefly  a  recent  generalization  of  (26)  by  Gripenberg 
(34),  which  will  also  be  jsed  in  Chapter  4. 

Our  method  involves  reducing  the  study  of  (V)  to  that  of  an  equivalent  functional 
differential  equation  of  the  form 


(EDS) 


^  ♦  Au  a  G(u)  (0  <  t  <  T) 

,  B(0)  -  X  -  f(«)  , 


where  6  t  C([0,T]>  0(A))  ♦  L^(0,T;X)  is  a  particular  mapping,  and  developing  the  theory 
for  (E0«).  Our  results  are  also  directly  applicable  to  certain  lntagrodlffarentlal 


-Id- 


equations  studied  by  Mac  Camy  (58]  via  a  Galerkin  argument  which  necessitates  further 
restrictions. 


He  observe  that  if  b  5  1,  equation  (V)  is  equivalent  to  the  evolution  problem 

(E)  ♦  Au  a  f*,  u(0)  -  x  -  f(0),  (0  <  t  <  T)  . 

at 

Our  method  of  studying  (FOB)  consists  of  generalizing  known  results  for  (E)  due  primarily 
to  Benilan  (13] »  the  latter  are  reviewed  in  Section  2.  He  recall  also  that  the  initial¬ 
boundary  value  problem  for  a  linear  or  nonlinear  diffusion  problem  is  a  special  case  of 


i 

| 

♦ 

i 


(E). 

2.2.  Preliminaries  on  Evolution  Equations.  For  further  background  and  details  of  this 
section  we  refer  the  reader  to  (7],  [25],  [26].  Let  X  be  a  real  Banach  space  with  norm 
1*1.  A  mapping  A  t  X  ♦  2*  is  called  an  operator  in  Xi  its  domain 
D(A)  «  {x  «  X  i  Ax  *  and  its  ranee  RCA)  -  u  {Ax  >  x  e  DC  A)};  A  is  single-valued 
if  Ax  is  a  singleton.  An  operator  A  in  X  is  accretive  iff  ■  (1  +  XA)  '  is  a 
contraction  in  X  for  X  >  0.  It  follows  immediately!  A  is  accretive  iff 
(2.1)  l(xt  ♦  Xy1)  -  (i^  ♦  Xy2)l  >  lx1  -  x^l  for  y±  e  Ax±  Cl  -  1,2)  . 

An  operator  A  in  X  is  called  m-accretlve  iff  A  is  accretive  and  RCI  ♦  XA)  *  X  for 


X  >  0. 


He  shall  be  concerned  with  applying  some  known  facts  about  the  abstract  evolution 
equation 

(E^)  ♦  Av  a  g,  vCO)  -  x 

to  the  study  of  CFOE).  He  asume  throughout  that  g  *  L^CO.TjX),  T  >  0. 

Definition  2.1.  A  function  v  «  [0,TJ  ♦  X  la  a  strong  solution  of  (Eg)  on  [0,1]  if 
v(0)  -  x,  v  c  CC  [0,T]  iX)  n  h’^CO.TiX),  v(t)  «  OCA)  a.e.  on  [0,T]  and  there  exists 
v  «  Av  such  that  w(t)  t  AvCt)  and  v’Ct)  ♦  w(t>  »  g(t)  a.e.  on  [0,T). 

Definition  2.2.  v  t  [0,T]  ♦  X  is  a  weak  solution  of  (Eg)  _on  (0,T|  if  there  is  a 
sequence  {(v  ,g  ))^  «  C((0,I]|X)  x  l'co.TjX)  such  that  vR  is  a  strong  solution  of 
(EgJ  on  [0,T]  and  (v^.g,,)  ♦  (v.g)  in  C([0,T]>X)  *  l’cO.TiX). 
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ncm 


For  our  considerations  we  require  a  third  concept  of  solution  of  (Eg) ,  namely  the 
notion  of  integral  solution.  First,  let  |  ,  )j  t  X  x  X  ♦  R  be  defined  for  X  >  0  by 

|x,ylx  «  j  (lx  ♦  Xyl  -  Ixl)  , 

which  is  a  nondecreasing  function  of  X.  Define 

[x,y]  “  lim  [x,y] .  -  inf  [x,y]. 

X+0  A  X>0 

lx,y]  -  lim-  [x,y] .  «  sup  [x,y]  . 

XtO  A  X<0  A 

Thus  lx  +  Xyl  >  Ixl  for  X  >  0  iff  [x,y]+  >  0,  so  that  A  is  accretive  iff 

(2.2)  *X1  “  *2'  y1  "  y2*+  *  0  £or  yl  *  **1  * 

Definition  2.3.  v  »  [0,T]  ♦  X  is  an  integral  solution  of  (Eg)  on  (0,T)  if 

V«C([0,T]jX  and 

t 

(2.3)  lv(t)  -  xl  -  lv(s)  -  xl  <  /  [v(a)  -  x,  g(o)  -  y]+da 

for  t  >  s,  (t,s)  *  ’  [0,T] ,  x  «  D(  A)  and  y  e  Ax.  We  note  since  |[x,y]+|  <  lyl,  and 

since  g  c  t1(0,T;X),  the  integral  in  (2.3)  is  well  defined.  A  straightforward 

calculation,  see  [25],  shows  that  the  notion  of  integral  solution  only  makes  sense  \<hen 

A  is  accretive.  We  shall  apply  the  following  result  on  existence,  uniqueness,  dependence 

on  data,  and  regularity  about  integral  solutions  of  (Eg)  due  to  Benilan  [13]. 

Theorem  A.  If  A  is  m-accretive,  x  t  D(A) ,  and  g  e  L1(0,T;X)  then  (Eg)  has  a  unique 

_______  * 

integral  solution  v  «  C([0,T]|  0(A))  on'  [0,T],  and  if  v,v  are  integral  solutions  of 
(Eg),  (E“)  on  [0,T]  corresponding  to  initial  values  x,x  respectively  then 

*  .  t  . 

(2.4)  lv(t)  -  v(t)l  <  lx  -  xl  +  /  lg(o)  -  g(o)tdo,  0  <  t  <  T  . 

0 

Moreover,  if  g  *  BV([0,T]jX)  and  x  «  D( A) ,  then 

(2.5)  hr(C)  -  V(n)l  <  I?  -  nlOg(0+)  -  yl  ♦  var(g  «  [0,t] )) 

for  y  *  Ax,  and  0  <  C,  n  <  t,  t  e  (0,T).  In  particular,  the  Integral  solution  v  is 
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Llpschitz  continuous.  If,  In  addition,  X  is  reflexive,  then  v  is  a  strong  solution  of 
(Eg)  on  [0 ,T] . 


3.  Discussion  of  Existence  and  Uniqueness  Results.  Ke  shall  reduce  the  study  of  existence 
and  uniqueness  of  solutions  of  the  nonlinear  Volterra  equation  (V)  on  [0,T]  to  studying 
the  abstract  functional  differential  equation 


(FDE) 


Aa  9  G(u)  (0  <  t  <  T) 

u(0)  «  X  , 


where  A  is  a  given  m-aecretive  operator  on  X,  and  where  G  is  a  given  mapping 

G  :  C([0,T1)  D(  A) )  +  l’iO.T.-X)  . 

Let  v  »  H(g)  denote  the  unique  integral  solution  of  (E^).  A  solution  of  (FDE)  is  by 
definition  a  function  u  e  C([0,T]j  D( A)  )  such  that  u  »  H(G(u)).  By  analogy  with 
Definition  2.1,  we  say  that  u  is  a  strong  solution  of  (FDE)  on  [0,T]  if  u(0)  «  x, 
u  e  W1'1(0,T»X)  nc([0,T);  D(aT)  and  if  u*(t)  +  Au( t)  i  G(u)(t)  a.e.  on  [O.TJ. 

Let  b  e  R),  F  e  l'(0,T;x).  We  shall  say  that  u  is  a  strong  solution  of  the 

Volterra  equation  (V)  on  [0,T]  _if  u  e  I^fO.TrX)  ant1,  if  there  exists  w  e  L’dJ.TjX) 
such  that  w(t)  e  An ( t)  and  u(t)  +  b*w(t)  -  F(t)  a.e.  on  (0,T).  One  establishes  the 
following  equivalence  between  strong  solutions  of  (FDE-,  with  a  particular  G  and  strong 
solutions  of  (V): 

Proposition  3.1.  Let  be  AC([0,TJ;R),  b*  e  BV([0,TJ»R),  Ft  w1'1(0,T>X)  and  b(0)  -  1. 
Let  u  be  a  strong  solution  of  (V)  on  [0,T].  Then  u  is  a  strong  solution  of  (FDE)  on 
[0,T]  with  the  identifications; 

t 

(i)  G(u)(t)  -  f’(t)  -  r*f’(t)  -  a(0)u(t)  -  r(t)x  +  /  u(t  -  s)dr(s) 

0 

(ii)  x  «=  f(0) 

■ 

(iii)  a  -  b> 

( iv)  r  e  L1  (O.TsR)  is  defined  bj;  r  +  a*i  «  a  . 


(3.1) 


Conversely,  let  r  <  BV((0,T]>R),  f*  e  l’cO.TiX),  x  t  57 A)  and  G  be  given  by  (3.1) 
(i).  Let  u  be  a  strong  solution  of  ( FDE )  o£  [0,T].  Then  u  Is  a  strong  solution  of 
(V)  _on  [0,T],  where 

t 

f(t)  -  x  +  /  f ' ( s)ds 
0 

t 

b(t)  =  1  +  /  a(s)ds 
0 

a  -  a*r  -  r  . 

F  e  W1,1(0<T»X),  then  the  Volterra  equation  (V)  is 
equivalent  to  the  evolution  equation  (E^)  where  g  *■  f*  and  where  the  initial  value 
x  -  f(0). 

The  proof  of  Proposition  1  is  straightforward.  The  assumptions  on  b  and  F  permit 
differentiation  a.e.  on  (0,T)  of  a  strong  solution  a  of  (V).  The  differentiated 
equation  is  then  "solved'1  for  Au  by  means  of  the  resolvent  kernel  r  associated  with 
a  «  b* ,  see  (3.1)  (iv),  and  the  variation  of  constants  formula  for  Volterra  equations 
(631.  a  known  result  (1X1  yields  that  a  *  BV( (0,T) }%)  implies  that 

r  e  BV([0,T);R),  a  fact  which  is  used  in  arriving  at  the  formula  (3.1)  (i)  for  G(u).  The 
converse  is  proved  by  reversing  the  steps.  A  part  of  Proposition  3.1  which  motivates  our 
approach  is  contained  in  Mac  Camy  (581  who,  however,  then  studied  (FDE)  by  an  entirely 
different  approach. 

We  remark  that  here  we  have  chosen  to  define  the  resolvent  kernel  by  (3.1)  (iv), 
rather  than  by  r  +  a*r  “  -a  as  was  done  in  [26J.  This  is  more  convenient  for  the  theory 

in  Chapter  4,  and  only  causes  a  change  of  signs  in  the  formula  (3.1)  (i)  of  some  of  the 

terms  in  G(u)(t).  Recall  that  if  r  is  defined  by  (3.1)  (iv),  then  the  solution  of  the 
linear  Volterra  equation  w  +  a*w  -  v  is  given  by  w  -  v  -  r*v,  while  with  the  alternate 

definition  of  r,  w  would  be  given  by  w  -  v  +  r*v,  as  was  used  in  (261 . 


(3.2) 


(i) 

(ii) 

(iii) 


We  remark  that  if  b(t)  6  1  and 
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We  next  use  Bcnilan's  theorem  about  solutions  of  (E^)  to  obtain  some  general  results 
concerning  existence,  uniqueness,  dependence  on  data,  and  regularity  of  solutions  of  (FDE) 
of  independent  interest  and  use  them  to  deduce  corresponding  results  about  solutions  of 
(V). 

Theorem  3.2.  Assume  that  A  is  m-accretive.  x  e  D(A),  and  let 
G  t  C([0,T)>  0(A))  ♦  l'cO.T.-X)  satisfy 


lG(u)  -  G(v)l 


L  (0,t;X)  0 


<  /  Y(s)lu  -  vl  . 


L  (0 ,s;X) 


for  some  Y  «  L  (0,TjR  },  0  <  t  <  T,  and  u,v  «  C(  tO.Tl  ;D(A))  . 


Then  (FDE)  has  a  unique  solution  u  t  C([0,T]j  D(A))  _on  [0,T]. 

We  remark  that  assumption  (3.3)  implies  that  the  value  of  G(u)  at  t  *  [0,TJ 
depends  only  on  the  restriction  of  u  to  [0,t] .  The  idea  of  the  proof  is  very  simple. 
Let  v  ■  B(g)  denote  the  unique  integral  solution  of  (E^)  on  [0,T],  g  *  L*(0,TjX).  We 
seek  a  fixed  point  of  the  map  K  i  C(  [0,T]  ;D( A) )  ♦  C((0,T];  D( A) )  defined  by 
X(u)  ■  H(G(u) ) .  By  property  (2.4)  of  integral  solutions 


lK(u)(t)  -  X(v)  (t)  I  <  /  !G(u)(s)  -  G(v)  (s)  Ids 

0 


(0  <  t  <  T)  , 


for  u.v  *  C((0,T]i  0(A)),  u(0)  -  v(0)  •*  x.  Applying  assumption  (3.3)  it  is  now  an  easy 
matter  to  show  that  is  a  strict  contraction  on  C( I0,T] t  D(A) )  for  j  sufficiently 

large,  so  that  the  map  K  has  a  unique  fixed  point.  For  details  see  (26]. 

Onder  further  assumptions  one  can  apply  the  second  part  of  Beni lan'  »  theorem  to  obtain 
greater  regularity  of  solutions  of  (FDE). 

Theorem  3.3.  In  addition  to  the  assumptions  of  Theortm  1  assume  that  there  is  a  function 
k  i  (0.-)  ♦  (0,-)  such  that 


0.4) 


and  IG(u)(0*)l  <  k(R> 


(0  <  t  <  T) 


whenever  u  «  C((0,T);  D(A))  Is  of  bounded  variation  and  lul  m  <  R.  _If 

I*  CO.TtX) 

X  «  OCA)#  then  the  solution  u  _o£  (FOE)  is  tlpschitc  continuous  on  [0,T].  If  X  is 
also  reflexive,  then  the  solution  u  is  a  strong  solution  of  (FOE)  on  (0,T). 

For  the  proof  of  Theorem  3.3  one  defines  u^  :  (0,T)  ♦  X  by  Ug(t)x  and 
■*♦1  *  *<un>  *  H(G(un)),  n  ■  0,1,...  .  These  iterates  converge  uniformly  and  are  uniformly 
bounded  on  (0,T).  By  Benilan's  theorem  and  assumption  (3.4)  one  shows  that  there  exists  a 
constant  c  >  0  such  that 

t 

varfu  :  (0,tl)  <  c{1  ♦  /  var(u  :  [0,s])ds) 
n*i  •  n 

0 

for  0  <  t  <  T#  so  that  var(u  . .  t  (0#t))  <  c  exp(ct).  Thus  {var(u  s  [0,T])}  and  by 

n**1 1  ft 

(3.4)  (varCGCu  ))  :  [0,T]}  are  both  bounded,  and  {u  },  and  hence  also  u  -  unif  lim 

n  n 

Ujj#  Is  Lipschitz  continuous  on  [0,T].  For  more  details  see  [26]. 

Finally#  the  solution  u  of  (FDE)  depends  on  the  data  A,  G,  x  In  the  following 
sense: 

Theorem  3.4.  let  the  assumptions  of  Theorem  1  be  satisfied.  Let  m-accretlve  operators 

A^  in  X,  mappings  Gn  :  C([0,T)jX)  ♦  I,1  (0,T»X) ,  and  x^eofA)  be  given  for 

n  “  1,2, ...  •  Assume  that  the  Inequality  (3.3)  holds  for  G  replaced  by  Gn,  with  the 

same  Y#  for  n  ■  1,2,...,  and  u,v  «  C([0,T];  D(  A) ) .  For  u  e  C([0,T];  D(  A) )  assume 

that  lim  G  (u)  -  G(u)  in  L^CO.TjX),  lim  x  -  x  e  D(A)  #  and 
°  n««  n 

(3.5)  11m  (X  ♦  U  )“1*  -  (X  ♦  lA)""1*  (z  (  X,  X  >  0)  . 

bet  un  e  C(  [0,T] ;  dTT7>  be  solutions  of  (FDE)  on  [0,T]  with  A  replaced  by  A„,  G 

replaced  by  Cn,  x  replaced  bv  xn#  and  let  u  «  C([0,T);  D(  A) )  be  the  solution  of  (FDE) 

on  [0,T] •  Then  lim  u  -  u  in  C((0#T];X). 
n-*»  n 

The  proof  of  Theorem  3.4  follows  from  the  observation  that  under  our  assumptions  the 
mapping  K(A,x,G)(u)  -  H(A,x#C(u})  of  Theorem  3.2  has  the  property  that  In  the  iterate 
X^#  which  is  a  strict  contraction  for  some  j,  both  j  and  the  contraction  constant 
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only  on  the  function  1  of  (3.3),  and  the  latter  is  assumed  to  be  uniform  in  n;  for 
details  see  [26] . 

We  shall  next  apply  Theorems  3.2,  3.3,  3.4  to  study  the  nonlinear  Volterra  equation 

(V).  if  b  and  F  in  (V)  satisfy  the  assumptions  of  Proposition  1,  it  follows  from  the 

definition  of  G  in  (3.1)  (i)  that 

iG(uMt)  -  G(v)(t)l  <  { |r(0+ )  |  +  var(  r  :  tO,t]))lu  -  vl  _ 

L**(0,t»X) 

where  r  is  the  resolvent  kernel  corresponding  to  b'  «  a  (recall  that 
a  e  BV([0,T];R)  *=■==>  r  t  BV(  (0  ,T]  jR) ) .  Thus  assumption  (3.3)  of  Theorem  3.2  is  satisfied 
with 

Y(s)  «  |r(0)|  +  var(r  :  [0,s])  . 

Moreover,  if  f  e  8V([0,T);X),  (3.1)  (i),  (ii),  imply 

var(G(u)  :  [0,t])  <  C(1  +  var(u  :  t0,t]))  (0  4  t  <  T)  , 

and  IG(u) (0+ ) I  <  C,  where  C  is  a  constant  depending  on  f(0),  f,(0+),  var(f*  :  (0,T]), 
r(0+),  and  var(r  :  to ,‘T] ) »  thus  assumption  (3.4)  of  Theorem  3.3  is  satisfied. 

Let  X  >  0  and  def.'.ne  the  Yosida  approximation  of  the  nr-accretive  operator  A 

on  X  by 

Ax  -  ~  (I  -  Jx),  Jx  e  (I  ♦  XA)’1  . 

*x  s  X  *  X  is  Lipschitz  continuous  with  Lipschitz  constant  j-,  so  a  simple  contraction 
argument  shows  that  the  approximating  problem 

(V  UX  +  b*Vx  ■  f 

has  a  unique  strong  solution  ux  on  (0,T),  under  the  assumptions:  b  e  L^(0,T»R) ,  and 
f  t  L^ ( 0,T;X) .  By  Proposition  3.1  u^  is  a  strong  solution  of 

duX 

(FDEX)  —■  *  AXUX  -  G(UX),  Ux(0)  -  F(0)  . 

One  also  has  lim  (1  +  pA. )  'z  »  (I  ♦  yA)  ^z,  for  y  >  0,  z  «  X. 

XiO 

These  considerations  lead  to  the  following  result  about  solutions  of  (V). 


(iv)  Our  method  can  be  used  to  study  the  nonconvolution  Volterra  equation 


u( t)  +  /  b(t,s)Au( s)ds  »  f(t) 
0 


(0  <  t  <  T)  , 


where  A  and  f  are  as  in  Theorem  3.5,  provided  the  kernel  b,  which  is  defined  on  the 
region  {(t,s)  :  0  <  s  <  t  <  T},  is  sufficiently  smooth  and  b(t,t)  >  0.  The  technique 
for  doing  this  is  outlined  in  (26],  and  is  carried  out  in  detail  by  C.  Rennolet  [72].  For 
different  nonconvolution  equation  results  see  Gripenberg  [37] ,  [38] . 

(v)  In  Theorem  3.5  above  and  3.7  below  it  is  important  to  note  that  the  generalized 
(or  strong)  solution  u  e  C([0,T]:X),  and  therefore  u  f  tP(0,T jX)  for  p  >  1. 

(vi)  If  the  assumptions  of  Theorems  3.5  above  and  3.7  below  are  satisfied  for  every 
T  >  0,  then  the  conclusions  hold  on  [0,**). 

(vii)  The  relation  of  Theorem  3.5  to  other  literature  is  explained  in  [26]. 

(viii)  An  interesting  situation  not  covered  by  the  theory  discussed  above  arises  if 
the  kernel  b  in  (V)  has  the  property 


In  the  case  A  *  3 where  <f  is  convex,  l.s.c.  and  proper,  G.  Gripenberg  [36]  extends 
the  theory  by  replacing  the  assumption  b'  e  BV[0,T]  in  Theorem  3.5  above  by: 

there  exist  TQ  >  0,  cQ  >  0  such  that  if  0  <  t  <  Tfl, 

Var(b'»(t,T0])  <  cQ  log  t"1  . 

Recently  G.  Gripenberg  (34,  Theorem  2]  obtained  the  following  important  generalization 
of  Theorem  3.5  which  will  be  used  in  Chapter  4.  Such  a  result  was  established  by  clement 
and  Nohel  [IQ]  for  the  much  simpler  case  b^  i  0  below  and  A  a  linear  operator. 

Theorem  3.7.  Let  the  assumptions  concerning  f,  A  in  Theorem  3.5  be  satisfied.  Let  the 


kernel  b 


in  (V),  where  b.  satisfies  the  assumptions  of  b  in  Theorem  3.5 


and  where  bj  *  L  (0,T)  and  bj  is  positive,  nonlncreaslnq,  and  log  b2  is  convex  on 
(Q,T).  Then  the  Volterra  equation  (V)  has  a  generalized  solution  u  e  C([0,T];X). 


-2' 


Finally,  we  close  this  chapter  with  two  important  special  cases.  The  first  deals  with 
the  case  of  the  operator'  A  being  maximal  monotone  on  a  real  Hilbert  space  H  and  is  a 
direct  consequence  of  Theorem  3.5.  Recall  that  if  X  •  H,  A  is  m-aceretive  iff  A  is 
maximal  monotone. 

Theorem  3.8.  Let  b  satisfy  the  assumptions  of  Theorem  3.5  for  every  T  >  0.  Let  A  be 

maximal  monotone  on  H.  Let  f  c  w’^(0,“;H),  f<0)  e  D( A) .  Then  (V)  has  a  unique 

generalized  solution  u  f  C(  [£!,“};  D(A)).  If,  in  addition  f'  *  BV^  {  [0 ,»)  jH) , 

f (0 )  e  D(A),  then  u  e  w',2{0.»*H)  and  u  is  a  strong  solution  of  (V)  on  [<>,«■). 

-  ^oc  — —  - 

Remark  3.9.  By  another  theorem  of  Gripenberg  [34,  Theorem  1) ,  this  result  also  applies  to 
(V)  with  kernels  b  =  b^  +  bj  which  satisfy  the  assumptions  of  Theorem  3.7. 

The  second  special  case  of  Theorem  3.5  deals  with  X  •  H  a  real  Hilbert  space  and  the 
operator  a  =  3^  (the  subdifferential  of  <f),  where  <p  :  H  ♦  (-*»,♦*•]  is  a  proper,  convex, 
l.s.c.  function.  The  proof  of  the  next  result  follows  by  combining  Theorem  3.5  with  known 
results  for  evolution  equations  (for  details  see  [26,  Section  4]). 

Theorem  3.10.  Let  the  kernel  b  satisfy  the  assumption  of  Theorem  3,5  for  every  T  >  0. 

Let  A  -  3 w,  where  <t  t  H  *  (-•,-*—]  is  proper,  l.s.c.,  and  convex.  Let 

f  t  w1,2(0,»;H).  If  f (0 )  *  6 (V ) ,  then  (V)  has  a  unique  strong  solution  u  on  [0,») 
loc 

2  2 
such  that  /t  u'  «  Lloc<0,-»H);  if  f(0)  e  D(y>),  then  u'  e  Lloc(0,»;H). 

Remark  3.11.  In  a  different  direction  Kiffe  and  Stecher  [48]  study  existence  and 

uniqueness  of  L2(0,T)  solutions  of  (V)  in  a  Hilbert  space  setting.  They  assume  that 

f  *  L*(0,T»H)  and  they  ise  techniques  of  Barbu  [5]  and  Londen  [54J  to  obtain  their  results 

without  any  differentiability  assumptions  on  the  forcing  term  f,  but  at  the  expense  of 

drastically  restricting  the  growth  of  the  maximum  monotone  operator  A  in  (V).  In  fact, 

this  restriction  rules  the  important  possibility  that  A  is  a  nonlinear  differential 

operator  in  the  spatial  variables,  and  therefore,  their  results  cannot  be  applied  to  the 

physical  problem  in  Chapter  t. 

An  interesting  and  different  variant  of  (V)  was  recently  studied  by  Kiffe  [49].  He 
obtains  existence  of  global  solutions  of  the  equation 
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where  A  can  be  a  nonlinear  differential  operator  in  the  apace  variables  and  where  the 
perturbation  g  is  a  discontinuous  real  function  which  is  not  necessarily  monotone  and 
satisfies  certain  growth  conditions.  The  kernel  b  satisfies  assumptions  similar  to  those 
in  this  chapteri  while  b'  e  BV[0,T]  is  not  assumed,  certain  monotonicity  is  required  of 
b.  The  forcing  term  f  e  W  '*»  the  operator  A  -  where  the  function  is  as  in 
Theorem  3.10.  The  function  i  satisfies  a  compactness  assumption,  and  f(0)  e  ow. 

Remark  3.12.  Other  interesting  variants  of  the  abstract  equation  (V),  motivated  by  the 
heat  flow  problem  formulated  in  Chapter  1,  Section  2,  have  been  studied  by  V.  Barbu  (6] , 
(10],  and  by  H.  Attouch  and  A.  Damlamian  [4].  In  (8,10]  Barbu  generalises  the  dependence 
of  the  internal  energy  e  on  the  temperature;  this  leads  him  to  study  the  equation 

Bu  ♦  b*Au  »  f  (0  <  t  <  •)  , 

where  B  is  a  strictly  monotone  operator.  In  (4]  the  domain  of  the  operator  A  in  (V)  is 
allowed  to  depend  on  time  ti  for  the  heat  flow  problem  formulated  in  Chapter  1,  Section 
2,  the  temperature  u  is  prescribed  at  each  time  t  outside  a  body  fl(t)  in  x  space. 
Assuming  that  fl(t)  depends  smoothly  on  t,  the  temperature  inside  (l(t)  is  determined. 
Remark  3.13.  While  the  results  are  still  rather  incomplete,  an  interesting  study  of 
numerical  approximations  of  solutions  of  (V)  has  been  initiated  by  Mac  Camy  and  Weiss  [61] 
where  other  references  to  numerical  literature  may  be  found. 


3.1.  Introduction.  The  purpose  of  this  chapter  is  to  discuss  the  boundedness  and 
asymptotic  behaviour  as  t  ♦  •  of  solutions  of  the  nonlinear  Volterra  equation 


(V)  u( t )  *  (b*Au)(t)  »  f(t)  (0  <  t  <  •)  • 

The  setting  for  (V)  is  b  s  (0,“)  ♦  R  is  a  given  kernel,  A  is  a  (possibly  multivalued) 
siaximal  monotone  operator  on  a  real  Hilbert  space  H,  and  f  :  [0,~)  ♦  H  is  a  given 
function.  The  exposition  is  largely  based  on  a  forthcoming  paper  by  Clement,  Mac  Camy,  and 
Hohel  (20]. 

The  following  general  assumptions  will  be  assumed  throughout: 

(Hb)  b(  t)  m  B(t) ,  b(0)  >  0,  b^  >  0,  B,B*  e  I.1  ( 0 ,—)  : 

(Hm)  A  maximal  monotone  on  H  > 

(Hf)  f(t)  «  f„  +  F(  t) ,  F  t  W^((0,»»H)),  F*  t  L2(0,-»H),  f„  t  H  I 

here  '  -  d/dt,  H  is  a  real  Hilbert  space  with  scalar  product  (•,•)  and  norm  1  *  |  , 

and  W  denotes  the  usual  Sobolev  space.  The  special  case  of  (Hm): 


<v 


A  »  3  <p,  where  the  function  <f  :  H  ♦  (-",*•)  is  convex, 
lover  semicontinuous,  and  proper 


will  also  play  an  important  role  in  the  theory.  For  definitions  and  standard  results 
concerning  maximal  monotone  operators  and  the  special  case  of  a  subdifferential  the  reader 
is  referred  to  Brezis  [14]. 

He  remark  that  if  one  adds  the  assumption  B'  t  BV^^tO,*)  to  assumptions  (Hb) ,  (Hm) , 

and  (Hf),  then  by  Theorem  3.8  of  Chapter  2.  (V)  has  a  unique  generalized  solution 

v  e  C((0,»);  0(A)),  provided  f(0)  e  D( A) ;  if  also  F’  «  BV,  (0,»;H)  and  f(0)  5  D(  A) , 

IOC 

then  u  is  a  strong  solution  of  (V)  on  (0,»).  If  in  place  of  (Hm)  assumption  (H^)  is 
satisfied,  and  if  f  ( 0 )  <  D{*> ) ,  then  by  Theorem  3.10  of  Chapter  2,  (V)  lias  a  unique  strong 


solution  u  on  |0,“)  such  that  tu*  e  if  f(0)  e  D(*),  then 

»'  «  1*^(0, -»H). 

The  results  on  boundedness  and  asymptotic  behaviour  of  solutions  of  (V)  will  be 

derived  from  a  priori  estimates  obtained  from  the  equivalent  differentiated  form  of  (V): 

*V>  ~  ♦  b(0)Au  ♦  B**Au  t  T*  (0  <  t  <  •).  u(0 )  -  f(0)  . 

tie  shall  distinguish  two  cases:  (1)  A  satisfies  assumption  (Hn)  and  (ii)  A  »  3*  with 

d  satisfying  (H  ) .  Case  (i)  is  developed  in  Section  2,  while  Case  (ii)  is  treated  in 
v 

Section  3;  in  each  case  a  different  energy  method  is  used  to  deduce  suitable  a  priori 
estimates  under  appropriate  additional  assumptions  on  the  kernel  b  end  the  forcing  term 
t»  A  number  of  examples  illustrating  each  situation  is  presented,  in  particular,  the 
theory  developed  in  Section  3  is  used  in  Section  4  to  analyse  the  boundedness  and 
asymptotic  behaviour  of  solutions  of  the  heat  flow  problem  (1.4),  Chapter  1,  under  the 
physically  reasonable  assumptions  motivated  in  Chapter  1,  See.  2. 

The  reader  should  note  that  if  Assumption  (H^)  is  satisfied,  and  if  u  is  the  strong 
solution  of  (V)  of  Theorem  3.10,  Chapter  2,  then  (see  Brezis  [14,  leant  3.3]  p(u(t))  is 

absolutely  continuous  and  one  has  the  "chain  rule”: 

|^*(u(t))  -  (w,  ||)  (w  (  3«Ku){t))  . 

Thus  a  plausible  energy  method  for  the  case  A  •  consists  of  taking  the  scalar  product 
of  (V*)  by  w  «  3v(u(t] )  for  any  solution  u  and  integrating  over  an  arbitrary  interval 
(Q,T)«  Indeed,  this  method  is  used  in  Section  3.  Unfortunately,  there  is  no  analogue  for 
the  chain  rule  when  A  is  a  maximal  monotone  operator,  but  A  *  3v  for  some  proper, 
convex,  l.s.c.  function  g.  For  this  reason  the  development  of  the  theory  in  Section  2  is 
leas  direct  in  that  the  a  priori  estimates  are  derived  from  an  equivalent  equation  to  (▼') 
resulting  essentially  from  applying  Proposition  3.1  and  Remark  3.6  (ill).  Chapter  2,  to 
(V’),  and  then  using  a  different  energy  method  to  obtain  the  a  priori  estimates. 

Zt  should  be  noted  that  Corollaries  2.4,  2.5,  and  Theorem  2.6  of  Section  2  may  be 
viewed  as  natural  generalisations  to  Hilbert  space  of  earlier  results  of  levin  [51)  and 
lond an  [S3]  which  describe  the  limiting  behaviour  as  t  ♦  •  of  solutions  of  (V)  in  the 
scalar  ease  in  which  the  operator  A  is  a  real  function. 


2.  Boundedness  and  Asymptotic  Properties  when  A  Is  Maximal  Monotone.  Throughout  this 
Motion  we  assume  that  assumptions  (Hb),  (HR),  and  (Hf)  are  satisfied  and  that  u  is  a 
strong  or  generalized  solution  of  the  Volterra  equation  (V)  on  [0  ,•*) .  As  explained  in 
Chapter  2  (V)  is  equivalent  to  the  Cauchy  problem 

(V*)  ^4  b(0)Au  4  B*  *  Au  »  F*  (0  <  t  <  -),  U(0)  -  f(0)  . 

let  k  be  the  resolvent  kernel  associated  with  B',  defined  to  be  the  unique  solution  of 

the  linear  Volterra  equation 

si  rfi 

(k)  b(0)k(t)  +  (B**k)(t)  -  -  a.e.  for  0  <  t  <  -  » 

by  standard  results.  Miller  (63],  assumption  (Hb)  implies  that  k  (  L^^fO,*"). 

We  now  use  the  method  of  Proposition  3.1  and  Remark  3.6  (iii)  of  Chapter  2  to 
transform  (V*).  Regarding  (V')  as  a  "linear"  equation  for  Au,  the  variation  of  constants 
formula  for  Volterra  equations  (63]  and  an  integration  by  parts  show  that  (V)  (and  hence 
also  (V))  is  equivalent  to  the  Cauchy  problem 

(2.1)  b(0T  dt  +  dt  <k*Ul  4  *“  »  ft  (0  <  t  <  -),  nCO)  -  f(0)  , 

where  f  1  i  (0,*)  ♦  H  is  the  function  given  by  either 

«.2>  f,(t)  -^-P'ft)  4  f(0)k(t)  4  (k*F*)(t)  (0  <  t  <  •> 

or 

(2-3)  ^(t)  -  bjoTF’<t)  4  k(0)f (t)  4  (k'*f)tt)  (0  <  t  <  •)  . 

We  shall  use  an  energy  method  based  on  taking  the  scalar  product  of  (2.1)  by  u,  and 
also  by  /t  u,  and  we  obtain  a  priori  estimates  by  integrating  over  an  arbitrary 
Interval  [0,T].  We  will  first  state  the  general  result  for  (2.1)  and  then  interpret  it 
for  (V). 

Theorem  2.1.  let  u  be  a  strong  or  generalized  solution  of  the  Cauchy  problem  (2.1)  on 


12.6)  C  ♦  «  >  0  . 

(a)  £f  f1  €  L2(0,«jH)»  1  then  u*  l"(0,*»iH)  ^  L1(0,»}H)| 

(b)  If  also  /t  k*  e  l'jO,-)  and  /t  f^  «  L2{0,«}H),  then  /t  u  e  L*(0,«iH)  H  l2(0,«*jB) . 
Consequently.  |u(t)|  “  o{^)  **.  *  *  m  and  u(t)  ♦  0  strongly  as  t  ♦ 

We  remark  that  no  claim  is  made  that  the  rate  |u(t)|  "  o(^)  *8  *  *  "  18  optimal. 
The  coercivity  assumption  (2.4)  concerning  the  maximal  monotone  operator  A  is 
natural  for  the  problem  in  light  of  comparable  assumptions  in  evolution  equations. 
Assumption  (2.S)  and  the  hypotheses  concerning  k.  k*  will  be  justified  in  Lemmas  2.2, 

2.3,  below.  Two  different  classes  of  kernels  b  in  (V)  are  considered,  each  of  which  lead 
to  the  energy  inequality  (2.S),  the  first  with  n  *  0,  the  second  with  «l  >  0,  and  for 
each  of  which  /t  k*  «  L1  (0 ,*•) .  These  technical  lemmas,  together  with  appropriate 
assumptions  on  the  forcing  function  f  in  (V) ,  permit  an  easy  interpretation  of  Theorem 
2.1  for  solutions  of  (V).  This  will  be  done  in  Corollaries  2.4  and  2.S  below.  The  proof 
of  Lemma  2.2  appears  in  Appendix  1.  Lemma  2.3  is  an  extension  of  a  result  of  Mac  Caiay  [59] 
which  in  its  present  form  was  recently  established  by  M.  Tangredi  [82] . 

Lena  2.2.  (a)  Let  b  iiatisfy  assumption  (Hb)  with  b_  >  0,  and  let  b  satisfy  the 

frequency  domain  condition 

4b 

<F|  there  exists  4  >  0  such  that  b—  ♦  Inf  [rh  Xm  B(in)]  >  6,  where 
•  •  h«*' 

B(iti)  “  /  exp(-int)B{t)dt.  Then  the  resolvent  kernel  k  of  B*  satisfies 

0 

It  *  l/(0,«). 

(b)  If  also  8'  €  L2(0,-:,  then  k  e  L2(0,-)j  if  also  B"  «  L1  (0,—) ,  then 

k*  «  L^O,-). 

(c)  If  the  assumptions  of  (a)  are  satisfied,  B"  «  L*{0,*»),  and  B  is  a  kernel  of 
positive  type  on  (0,«),  then  for  every  T  >  0  and  for  every  w  «  L2(0,T) 

*  A 

/  «(t)  (k*w) (t)dt  >  0  . 

•  dt 


33 


|d)  If  the  assumptions  of  (a)  and  (b)  are  satisfied,  and  /t  B*  *  l'(0,«*)  ^  L2(0,*»), 
/i  B-  *  L^O.-),  then  /tki  L^O,-)  n  l2(0,«),  and  /t  k*  «  L^O,®). 

Lemma  2.3.  Let  b  satisfy  assumption  (H^)  with  ba  "  0,  and  let 

(1)  tVB)  t  L1  (04**)  (j  -  0,1,2}  m  *=  0,1, 2, 3).  t3B  «  L^O.-)  , 

(H)  B  be  strongly  positive  on  (0,**)  . 

Let  k  be  the  resolvent  kernel  of  B * .  Then; 

(a)  k  t  c’lO,®)  ; 

m 

(b)  k(t)  »  k„  +  K(t),  ka  •>  (/  B(t)dt)'1  >  0  ,  K(B)  «  l\o  ,•)  (»  -  0,1,2)  » 

0 

(c)  If  also  B,B* ,  /t  B,  /t  B*  «  L  (0,«)  one  has  K,  /t  K  <  L  (0,«)  » 

(d)  •  for  every  T  >  0  and  for  every  w  r  L*(0,T)  there  exists  n  >  0  such  that 


as  . 

-nioB(in)  «  n  /  sinnt  B(t)dt  >  o  (n  «  R)  • 

0 

Thus  if  ba  >  0  la  any  constant,  b(t)  «  ba  ♦  B(t)  satisfies  the  frequency  domain 
condition  (F)  with  5  ■  ba,  and  (see  Lemma  2.2(a))  k  e  L*(0,~).  If,  in  .addition,  B 
satisfies  the  remaining  smoothness  and  integrabllity  assumptions  of  Lemma  2.2,  all 
conclusions  of  Lemma  2.2,  and  assumptions  (2.5)  with  n  ■  0,  and  /t  k*  *  l'(0,«)  of 


Theorem  2.1  are  satisfied 


Consider  again  B  in  (2.7).  In  addition,  assume  that 


(2.8)  the  measure  dB'  has  a  nonzero  absolutely  continuous  pa rtf 
then  (see  (67,  Corollary  2.2]),  B  is  strongly  positive  on  (<>,“)  (for  example. 

It.  C(0,»),  (-1)kB(k)(t>  >0,  0  <  t  <  -,  k  -  0,1,2,  B*  (t)  JTO).  Thus  if  B  satisfies 
(2.7),  (2.8),  and  the  integrability  and  smoothness  assumptions  of  Lemma  2.3,  and  if 
b(t)  5  B(t)  (ba  “  0),  then  all  conclusions  of  Lemma  2.3,  assumptions  (2.5)  with  n  >  0, 
and  /t  k*  «  L1(0,»»)  of  Theorem  2.1  are  satisfied. 

Mext,  consider 

■  -X  t 

(2.9)  B( t)  "  J  B  e  3  cos  w  t  (B  >  0,  X  >  0,  w  e  R) 

j«1  3  3  3  3  3 


with  strict  inequalities  holding  for  at  least  one  j  (if  »j  “  0,  j  ■>  1,...,m,  B 
satisfies  both  (2.7),  (2.8)).  This  function  B  is  strongly  positive  on  (0 ,“)  (see 
(67] ) ,  since  by  direct  calculation 


R.  B(in)  -  j  l  B  X  (~ - 

•  i«1  33  X2  + 


(n  -  v. ) 


x;  «•  (n  ♦  wa)‘ 


r)  (n  »  R) 


Moreover,  B  satisfies  all  other  assumptions  of  Lemma  2.3.  Thus  if  b(t)  -  B(t) 
(b—  ■  0),  all  conclusions  of  Lemma  2.3,  assumptions  (2.5)  with  n  >  0  and 
/t  k*  e  L3(0,«)  of  Theorem  2.1  are  satisfied. 

For  the  kernel  B  in  (2.9)  one  has 


-i»Xm  B(in) 


2,  2  .  ,2  2, 

?  n  "  (n  ~  V 

i-1  3  (X2  ♦  -  n2)2  ♦  4X2n2 


(n  «  R)  . 


Thus  b(t)  ”  bm  *  B(t),  where  ba  >  0  is  any  constant,  satisfies  the  frequency  domain 
condition  (F)  of  Lemma  2.2  if  X^  >  w^  (J  ■  1,. ..,■).  Evidently,  b  is  a  kernel  of 
positive  type  on  (0,*).  Therefore,  if  b(t)  ■  b—  ♦  B(t),  b—  >  0,  B  defined  by  (2.9)  with 


-35' 


Hj  >  W j  ( j  -  1 . . 

(2.5)  with  n  ”  0  and 
Incidentally,  if 


all  conclusions  of  Lemma  2.2  (but  not  of  Lemma  2.3) 
/t  k*  t  L1  (<>,••)  of  Theorem  2.t  are  satisfied. 

>  0  is  any  constant,  and  if 


assumptions 


(2.10) 


B(t) 


■  -X  .t 
£  Be  3  sin  w  t 
i*1  3  3 


(Xi  >  0,  >  0) 


with  strict  inequalities  holding  for  at  least  one  j,  then  the  frequency  domain  condition 
(H  of  Lemma  2.2  is  satisfied  with  6  =  b^.  However,  such  a  kernel  b  is  not  of  positive 
type. 

Lemmas  2.2  combined  with  appropriate  assumptions  on  A  and  f  yield  the  following 
easy  interpretation  of  Theorem  2.1  for  solutions  of  (V). 

Corollary  2.4.  Let  assumptions  (H^)  with  b^  >  0,  (Hjj,)  and  (Hj)  with  f—  arbitrary  be 
satisfied.  In  addition,  assume  that  b  satisfies  the  hypotheses  of  Lenna  2.2,  and  that 
/t  F'  «  L2  (0,*;H) .  Let  u  be  a  strong  or  generalized  solution  of  (V)  _on  [0,*»).  If  the 

coerclvity  assumption  (2.4'  holds  with  C  >  0,  then  u  and  /t  u  e  L  (0,";H)  L  (0,*»;H) 
and  u(t)  ♦  0  strongly  as  t  *  •• 

Indeed,  define  f.,  by  (2.2).  By  (B{)  and  Lemma  2.2  (k  e  L*(0,«»  )  D  L2(0,“), 
k*  C  L*(0,"))  one  trivially  has  f^  e  L2(0,*;H).  By  Lemma  2.2  one  also  has 
/t  k  «  L1  (<),«•)  n  L2  (0 ,»)  and  /t  k*  «  L*  ( 0 ,") .  These  together  with  the  assumption 
/t  r*  «  L2(0,**;H)  used  iu  (2.2)  show  that  /t  f1  «  L2(0,*»H);  the  fact  that 
/t  (k*F*)  e  L2(0,»jH)  in  (2.2)  follows  from  the  straightforward  estimate 
t  t 

/  t|/  k(t  -  s)F'(s)ds|2dt  <  2lkl%  l/tF'l2 
0  0  L*(0,*)  L‘C(0,-,H) 

4  2l/t  kl2  IF' I2,  (VT  >  0)  . 

1(0,-)  L  (0,->H) 

By  Lecxaa  2.2  again,  (2.S)  holds  with  n  -  0.  Thus  if  c  >  0  in  (2.4),  the  result  of 
Corollary  2.4  follows  by  applying  Theorem  2.1. 

Lemma  2.3  combined  with  appropriate  assumptions  on  A  and  f  yield  a  different 
interpretation  of  Theorem  2.1  for  solutions  of  (V). 
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Corollary  2.5.  Let  assumptions  (H^)  with  “  0,  (HR)  and  (Hj)  with  f—  ■  0  be 
satisfied.  In  addition,  assume  that  b(t>  -  B(t)  satisfies  the  assumptions  of  Lemma  2.3, 
and  that  f  5  F  also  satisfies  f.  /t  f,  /t  f*  «  L2(0,“;H).  Let  u  be  a  strong  or 
generalized  solution  of  (V)  on  (0,«).  if  the  operator  A  satisfies  the  coerclvlty 
assumption  (2.4)  with  c  >  0  (or  even  c  >  -n,  where  n  >  0  is  the  constant  in  Lemma 
2.3d),  then  u  and  /t  u  t  L*(0,»»H)  ft  L2(0,-jH)  ,  and  u(t)  ♦  0  strongly  as  t  ♦ 

The  proof  of  Corollary  2.5  is  similar  to  that  of  Corollary  2.4,  except  that  fj  must 
now  be  defined  by  (2.3),  and  Lemma  2.3  is  used  in  place  of  Lemma  2.2.  Note  also  that  the 
additional  assumptions  concerning  f,  /t  f  are  essential. 

The  important  case  b^  -  0  in  (H^) ,  b  =  B  satisfying  the  assumptions  of  Lemna  2.3, 
and  tm  *  0  in  (Hj)  is  not  covered  by  Corollary  2.5.  In  this  situation  Theorem  2.1  must 
be  modified  in  the  following  manner. 

Theorem  2.6.  Let  the  assumptions  (Hb)  (b_  -  0),  (HB) ,  (Hf)  with  f_  arbitrary,  and  the 
assumptions  of  Lemma  2.3  be  satisfied.  In  addition,  assume  that  F,  /t  P, 

/t  f'  *  L2(0,«jH).  Let  u  be  a  strong  or  generalized  solution  of  (V)  on  (0,“),  let 
b*  the  unique  solution  of  the  limit  equation  corresponding  to  (V )  i 

(V^)  .  nm  *  {/  B(t)dt) Aua  a  f,  . 


Let  the  operator  A  satisfy  the  coerclvity  condition: 


(2.11) 


if  ▼  Au  and  v—  e  Au^  and  T  >  0,  then 
T  T  , 

/  (v(t)  -  vm,  u( t)  -  u  )dt  >  c  /  |u(t)  -  uj  dt 

0  0 

for  some  c  >  0  (c  >  -q  Is  sufficient!  see  Lssm  2.3d)  . 


Then  u  -  u—  and  /t  (u  -  u— )  *  L  n  L2(0,*iH);  consequently  u(t)  ♦  strongly 

S£  t  ♦  •  and  |u(t)  -  u^l  ■  o(^)  as  t  ♦  •. 

Bemark  2.7.  Since  b  2  B  satisfies  the  hypothesis  of  Lemma  2.3,  B  is  strongly  positive  on 

m 

(0,*) ,  and  therefore  /  B(t)dt  >  0.  since  the  operator  A  is  maximum  monotone  on  H  the 
0 


-37- 
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limit  equation  (VL)  has  a  unique  solution  for  any  tm  «  H»  In  particular,  if  f—  “  0, 
u—  ■  0  and  in  this  case  Theorem  2.6  reduces  to  Corollary  2.S. 

Corollaries  2.4.  2.S  and  Theorem  2.6  together  form  the  natural  generalisation  to 
Hilbert  space  of  corresponding  scalar  results  for  (V)  due  to  Levin  (SI)  and  Londen  (S3). 

Sketch  of  Proof  of  Theorem  2.1.  (a)  Take  the  scalar  product  of  (2.1)  with  u  and 

Integrate  from  0  to  T.  Using  (2.4),  (2.S)  we  obtain 


2b(0) 


|u(t ) |  +  (n 


♦«)/ 

o 


(u(t>|  dt  < 


2b(0) 


If CO  1 1 


+  / 
0 


(ft(t),u(t))dt 


Since  n  ♦  c  >  0,  and  f,  t  L^(0,<*;K),  the  assertion  (a)  follows  by  standard  estimates 

i  i  * 

t 


|  (b)  Next  take  the  scalar  product  of  (2.1)  with  tu  and  integrate  from  0  to  T.  An 

! 

integration  by  parts  yields 


(2.12) 


T 

2b(0) 


iu(t»r 


t(u(t),  (k*u) (t) )dt  +  c 

<  2b(0)  f  ,u(t)l  dt  +  / 


/  t|u(t)|2dt 
0 

(tft(t),u(t))dt  . 


A  straightforward  calculation  shows  that  (see  (20)) 

T  . 

/  t(u(t) ,  (k*u)(t))dt  -  I  +  J  , 


where 


T  T 

X  •  /  </t  u(t),  f-  (k*/t  u) (t) )dt  >  n  / 

0  Z  0 


t|u(t)|  dt 


9  -  /  (/t  u(t),  /  k»(t  -  t)(/t  -  /r)u(T)dt)dt  . 


and  where 


Using  ft  -  ft  <  ft  -  1,  0  <  t  <  t,  u  «  L  (0,-j  H),  and  /t  k‘  c  L  (0,»)  yields 


|J|  <  (/  t|u(t)l2dt)1/2  lul 


/  /t|k* (t) |dt  . 


L  (0,«jH)  0 


Substitution  of  the  estimates  for  I ,  J  in  (2.12)  gives  the  final  inequality 


5^07  *'«**'/  .»■= 

T 

♦  lul  l/t  k* I  (/  t|u(t)|2dt)1/2+  l/t  f  I  lul 

L2(0,-|H>  L  (0,«)  0  1  L  (0,-»H)  L2(0,-;H) 

The  conclusion  -ft  u  «  l"*(0,*jH)  n  l2(0,»jH)  follows  by  standard  estimates  using 
C  ♦  n  >  0,  u  *  Ii2(0,«jH)  by  (a),  and  the  assumptions  /t  k*  e  L^(0,")  and 
/tfj«  L2(0,»jH).  This  completes  the  sketch  of  the  proof. 

Proof  of  Theorem  2.6.  The  proof  will  be  reduced  to  that  of  Theorem  2.1  by  the 

m 

following  steps.  First  by  lemma  2.3  /  B(t)dt  *  k”1  >  0.  Therefore,  the  limit  equation 

0 

(?tl  can  be  written  in  the  form 


*•".  +  to-  »  k-f- 


which  is  the  same  as 


(2.13) 


i?o7  k u- ♦  h  (k#v  ♦  to- .» k-f- ♦  {k(t>  -  *->«- 


Next,  subtracting  (2.13)  from  (2.1)  gives 


(2.14) 


b(0)  dt 


jjT  (u  -  Mm)  *  k*(u  -  uM)  ♦  Au  -  Au^  a  T.  (t)  (0  <  t  <  •) 


where  by  an  elementary  calculation 


(2.15)  Pt(t)  F' (t)  k(0)r(t)  «t)fB  ♦  (k'*F)(t>  -  uJC(t)  . 

Leona  2.3  and  the  assumptions  concerning  F  clearly  imply  that  Fj  satisfies  the  same 
assumptions  as  fj  in  Theorem  2.1.  The  method  of  proof  of  Theorem  2.1  applied  to  (2.14) 

(2.15) ,  where  the  coercirity  assumption  (2.11)  is  used  in  place  of  (2.4),  now  yields  the 
needed  a  priori  estimates  for  u  -  u^  and  St  (n  -  us) ,  and  completes  the  proof. 

Example  2.8.  We  give  an  example  of  a  maximtsa  monotone  operator  A  in  (V)  which  is  not  a 
subdifferential,  and  for  which  the  theory  developed  in  this  section  is  applicable.  Let 


OCR"  be  a  bounded  open  set  with  smooth  boundary  30.  Let  K  be  the  Hilbert  space 
L2(ft).  Let  6  be  a  maximum  monotone  graph  with  0  t  0(0)  and  with  primitive  j  (i.e. 
0  “  Jj).  Let  Aj  be  the  operator  defined  by 

D(Aj)  «  (u  t  u«  H®(0)  n  h2(0),  0(u>  «  L2(Q)J  . 


AjU  “  -Au  +  0(u) 


(U  *  D( Aj ) )  . 


It  is  clear  [15]  that  Aj  is  maximum  monotone  on  H  since  A^  “  3<^ ,  where 
t  ■  ♦  (**•«•]  is  the  proper,  convex,  l.s.c.  function  given  by 


*j(u)  « 


4-  /  |Vu|2dx  +  /  j(u)dx  if  u  c  H^Jl)  and  j(u)  *  l\o> 


0  0 
+•»  otherwise  . 

3u 


r>  3U  1 

Define  L(u)  •  £  b  -r —  (b  *  R,  u  e  H  (0).  By  the  divergence  theorem  L(u)  is  sx>notone 

i-1  1  ”i  1  0 

and  (u,L(u))  ”  /  uL(u)dx  -  0.  Finally,  following  Pazy  (70,  Ex.  3.5]  define 
0 

A'AjU, 

By  a  perturbation  theorem  of  Crandall  and  Pazy  [27] ,  A  is  maximum  monotone  on 

H  “  L2(0),  and  by  an  easy  calculation  using  Green's  theorem  and  the  Poincare  inequality 

there  exists  a  constant  e  >  0  such  that 

(An.u)  ■  (A.u,u)  *  -  /  uAudx  ♦  /  u0(u)dx  >  /  |Vu|2dx  >  elul2  . 

0  0  0  L(Q) 


Thus  A  satisfies  the  coercivity  assumption  (2.4)  for  every  T  >  0. 

Bemark  2.9.  The  concept  of  strong  positivity  of  a  kernel  plays  an  important  role  in 
stability  theory  for  Volterra  equations  (see  Halanay  [40],  Kohel  and  Shea  (67J,  Staffans 
P4],  [75],  [76],  [77],  [78],  Crandall,  Londen  and  Nohel  [28]).  As  we  have  seen  kernels 
B  satisfying  (2.7),  (2.8),  as  well  as  oscillatory  kernels  B  of  the  form  (2.9)  are 
strongly  positive  on  (0,»).  If  b(t)  -  +  B(t)  where  tom  >  0  and  B  is  strongly 

positive  on  (0,*)  one  can  extract  some  information  about  the  behaviour  of  solutions  of 
OT)  as  t  ♦  •.  In  particular,  one  can  establish  ue  L2(0,»jH)  (Theorem  2.1  part  (a))  by 
another  energy  method  directly  from  (v)s 


Proposition  2.10.  Let  the  general  assumptions  (Hjj) ,  (Hm) ,  (Hj)  be  satisfied.  Let  u  be  a 
strong  or  generalized  solution  of  ( V >  en  [0,“).  I £  F  r  L2(0,";H),  if_  B  la  strongly 
positive  on  (0,~),  if  the  coerclvltv  assumption  (2.4)  is  satisfied  with  e  >  0/  and  if 
. f  »  0  whenever  •  0 ,  then  u  t  L2 ( 0 , H ) . 

Proof.  Let  v  e  Au  and  let  0  <  T  <  •  be  given.  Take  the  scalar  product  of  (V)  by  v 
and  integrate  from  0  to  T  obtaining 

T  b-  T  2  T 

(2.16)  /  (v(t),u(t))dt  +  m-  !/  v(t)dt|  +  Q.tvrTJ  <  /  ( v(t) ,f (t) )dt  , 

0  2  0  B  0 

where 


Q-IviT]  -  /  (v(t),(B*v)(t))dt  . 

B  0 

Since  B  e  1^(0, m)  is  strongly  positive  and  F,F*  «  L'tO.^fH),  a  result  of  Staffans 
[70,  Proposition  4.1)  shows  that  there  exists  a  constant  Y  >  0  such  that 


1/  (v(t),F(t))dt|  <  Y(QB[ViT)} 

0 


1/2 


Using  this,  the  coercivity  assumption  (2.4)  (e>0),  and  the  obvious  estimate 


b_  T 


1/  (v(t),fjdt|  <  If  v(t)dt|2  +  IfJ*  (b.  >0) 


in  (2.16)  yields  the  inequalities 


b_  T 


iV2 


l*J‘ 


(2.17)  e  /  |u(t>|  dt  +  —  1/  v(t)dt|  +  Q  [vit]  <  YiCjlvjT)}  ♦  -5 -  (b.  >  0) 

0  0  • 


and 
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(2.18)  e  /  |u(t)|2dt  +  O  [v,T]  <  Y£ehtv;T)}1/2  (b,  -  0.  f„  -  0)  . 

0 

The  result  of  Proposition  2.9  now  follows  from  (2.17),  (2.18)  by  completing  the  square  on 
the  CB  terms. 

Incidentally,  we  have  also  shown  that 

(2.19)  sup  G_[v;T]  <  •  (v  e  Au)  , 

T>0  ^ 


and  if  b  >  0,  also  that 


j 


T 

sup  if  v(t)dt|  <  “  . 

T>0  0 

It  follows  from  (2.19)  using  another  result  of  Staffans  (see  Crandall,  Nohel ,  Londen  [28, 

Lemma  3.1]  in  a  Hilbert  space  setting)  that 

sup  | (B*v) (T) |  <  •  (vs  Au)  . 

T>0 

Unfortunately,  there  appears  to  be  no  direct  way  to  establish  also  that  u(t)  is  uniformly 

continuous  on  [(),«•)  if  assumption  (Hn)  is  satisfied  and  A  *  3 s?  (if  A  «  3*5  this  can  be 

2 

done  as  in  Theorem  3.1  below).  The  uniform  continuity  together  with  u  e  L  (0  would 

imply  that  u  ♦  0  as  t  ■»  •  strongly  in  H.  This  provides  at  least  one  motivation  for 
the  indirect  method  of  Theorem  2.1. 

3.  Boundedness  and  Asymptotic  Properties  When  A  »  3*.  Let  the  general  assumptions  (Hb) , 
(H^),  (Hf)  be  satisfied  and  let  u  be  a  strong  or  generalized  solution  of  (V)  or  (0,«>). 

In  this  section  we  shall  obtain  different  boundedness  and  asymptotic  results  for  the  case 
A  “  i<e,  and  when  bm  >  0  in  (H^) .  These  results  which  are  motivated  by  the  physical 
problem  discussed  in  Chapter  1,  Section  2  are  deduced  from  a  priori  estimates  which  ure 
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1 

i 


obtained  directly  from  the  equivalent  Cauchy  problem: 

(V*)  ♦  b(0 )Au  ♦  B*  *  *u  »  P*  (0  <  t  <  •),  u(0 )  *  f(0)  . 

at 

Theorem  3.1.  Let  the  general  assumptions  <H^)  with  b—  >  0,  (H^)»  (H^)  be  satisfied  and 

let  u  be  a  strong  or  oenerallted  solution  of  (V)  on  10,*).  If  the  kernel  b  satisfies 
the  frequency  domain  condition  (F)  of  Lemma  2.2  and  It 

(3.1)  Inf  rtz)  >  —  , 

z«H 

then 

(3.2)  Sup  v»(u(t))  <  »  i 

0<t<- 


lf  V  e  3^(u),  then 

(3.3)  V  e  L2(0,»jH)  , 

(3.4)  |j:  «  » 

and 

(3.5)  u  la  strongly  uniformly  continuous  on  to,*)  . 

If  also  lira  tf(u)  ■  +•,  then 
!«!♦* 

(3.6)  Sup  |u(t)|  <  •  , 

0<t<- 

and 

(3.7)  lira  ?(u(t))  »  <pm  ■  Inf  v(r)  exists  . 

t***  *«H 

Moreover,  If  the  Inclusion  3c(v)  »  0  Implies  w  -  0,  then 

(3.8)  u(t)  — x  0  .  (weakly)  ai  t  *  •  . 

The  frequency  domain  condition  (F)  la  satisfied  by  several  classes  of  kernels  b  with 
bm  >  0  as  was  seen  in  Section  2  (see  examples  of  b  -  b—  ♦  B  with  B  given  by  (2.7), 

(2.9) ,  (2.10)).  Thus  Theorem  3.1  generalizes  a  recent  result  of  S.  O.  London  [54, 

Corollary  2]  and  a  result  of  V.  Barbu  (6,  Theorem  2). 

The  assumptions  concsmlng  V  in  Theorem  3.1  are  not  sufficient  to  obtain  strong 
convergence  of  u(t)  to  zero  as  t  ♦  •.  For  this  result  one  needs  the  coercivlty 
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condition  (2.4)  with  C  >  0.  If  (2.4)  is  satisfied  with  v  e  3y(u)  it  is  a  standard 
result  (see  Brezis  (14])  that  the  inclusion  3*?(w)  >  0  has  w  -  0  as  the  only  solution, 
and  that  0  e  DO*?).  Then  the  definition  of  the  subdifferential  (14]  implies  that 

(u)  >  *(0)  (u  *  H)  , 

and  therefore  assumption  (3.1)  of  Theorem  3.1  holds.  This  motivates  the  following  results 
which  complement  Corollary  2.4  for  the  case  A  "  3^.  Kote  that  in  Theorem  3.2  below  only 
the  frequency  domain  condition  (F),  but  not  the  assumption  that  B  is  a  kernel  of  positive 
type  (see  Lemma  2.2),  is  needed.  Also  note  that  here  the  assumption  on  F  is  less 
restrictive. 

Theorem  3.2.  Let  the  general  assumptions  (H^)  with  b^  >  0,  (H^),  (Hj)  be  satisfied,  and 

let  u  be  a  strong  or  generalized  solution  of  (V)  on  (0,“).  Let  b  satisfy  the 

frequency  domain  condition  { F) ,  and  for  v  e  3<p(u)  let  the  coerclvltv  condition  (2.4)  with 

C  >  0  be  satisfied.  Then  conclusions  (3.2 )— (3.5)  of  Theorem  3.1  hold,  and 
s 

u  (  L  (0 ,**>H) ,  which  implies  that  u(t)  ♦  0  strongly  as  t  ♦ 

Remark  3.3.  Zf  b(t)  5^)0  in  (V),  a  case  not  excluded  in  Theorems  3.1  and  3.2,  the 

above  theorem  and  its  proof  yield  a  simple  boundedness  and  asymptotic  behaviour  result  for 
the  evolution  equation 

~  ♦  b^S^fu)  a  g,  u(0)  »  ufl  , 

where  g  ■  F'i  compare  Brezis  (14,  Theorem  3.11]  where  g  e  L1(0,»>H). 

Remark  3.4.  If  the  coercivity  condition  (2.4)  with  c  >  0  and  A  »  is  replaced  by  the 

■ore  general  condition:  for  every  T  >  0  there  exists  c  >  0  such  that 

Y  Y 

(2.4*)  if  ▼  «  o(u),  then  /  (y(t),u(t)  -  *)dt  >  /  |u(t>  -  z|2dt 

0  0 

for  some  s  *  H,  then  it  is  easy  to  show  that  the  inclusion  9^(v)  a  0  has  w  »  s  as  the 
only  solution,  and  that 

y(n)  >  y(z)  (u  «  H)  . 


Then  the  method  of  proof  of  Theorem  3.2  easily  yields  that  u(t)  ♦  s  strongly  es  t  ♦ 


tewrlt  3.5.  In  Theorem  3.1  and  3.2  the  assumption  bm  >  0  in  (Hfa)  is  crucial)  If  b—  ■  0 
the  frequency  domain  condition  (F)  cannot  be  satisfied  (see  examples  (2.7),  (2.9),  (2.10)) 
for  any  5  >  0.  On  the  other  hand,  in  these  theorems  tm  in  (Hf )  is  arbitrary  and  the 
case  f#  ■  0  is  not  ruled  out,  provided  ba  >  0.  If  b^  •  0  in  (Hfa) ,  one  can,  of  course, 
still  apply  Corollary  2.5  if  fm  -  0,  and  Theorem  2.6  if  fm  *  0,  with  A  -  3*. 

Proof  of  Theorem  3.1.  (a)  Let  0  <  T  <  •  be  arbitrary;  take  the  scalar  product  of 
(V*)  with  v  «  3*(u)  and  integrate  over  [0,T).  Using  ~  V(u(t))  «  (v(t),  j  (t))  a.e. 
(see  Brezis  (14])  one  obtains 

T  ,  T 

(3.9)  p(u(T))  +  b(0)  /  |v(t)|  dt  ♦  Q_,  (vjT]  ■  /  (P'(t),v(t))dt  ♦  o(f(0))  , 

0  ^0 

where 

T 

Q_, tw;T)  -  /  (v(t) ,B,#v(t) )dt  . 

0 

We  next  apply  a  frequency  domain  method  (see  Nohel  and  Shea  [67])  to  Qg. .  Define 
v^(t)  -  w(t)x(0,T]  and  its  Fourier  transform 

m  m 

▼T(n)  -  /  s"1,'tvT(t)dt  . 

Xxtend  B*  evenly  to  (— ,0)  by  B'(-t)  •  B'(t)  (0  <  t  <  •).  In  the  following 
calculation  use  is  made  of  the  hypothesis  B'  <  L1{0,«),  and  the  Parseval  and  convolution 
theorems: 

T  T  T 

Q_,(V»T]  -  /  (v(t),B'»v(t))dt  -  ^  /  (v(t),  /  B'(t  -  T)v(T)dT)dt 

0  2  0  0 

•  •  •  •  • 

•  j  /  <pTtt),  /  b* (t  -  T>vt(t>dT)dt  ■  I2  »'Cn)dn  . 

•  00  4» 

tines  S'  is  even,  B'(q)  *  2 Re  B*(in)  (n  «  K),  where  •  denotes  the  Laplace 
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transform.  The  assumptions  B.B*  c  L*  ((>,••) 

*  A  * 

B'(ln)  -  in  B(in)  -  B(0)  yield  Re  B'(in) 


and  the  familiar  formula 
* 

-  -nl»  B(in)  -  B(0) .  Therefore, 


Qjji  IvjT] 


—  J  |vT(n)l2t-nim  B(in) 


B(0)}dn  . 


*  2  1  “  2 

Substituting  this  result  into  (3.9)  and  using  /  |v(t)|*dt  «  —  /  |v  (n)i  dn,  as  well  as 

0  21  T 

b(0)  -  B(0)  “  bm,  the  frequency  domain  condition  (F),  and  Parseval's  theorem  again,  yields 


(3.10)  *(u(T))  +  «  /  |v(t)|2dt  <  |v>(f(0))|  +  /  |  (F*  (t)  ,v(t) )  |dt  <0  <  T  <  ->  . 

0  0 

2 

The  assumption  F*  *  l  (0,*),  Cauchy-schwarz  and  an  elementary  inequality  give  the 
estimate 

«  T  2  *  •  , 

(3.11)  F(u(T»  *  -  f  |v(t)  |  dt  <  |v>(f(0))|  ♦  ^  /  |F*  (t)  1  at  <  -  (0  <  T  <  •)  . 

0  0 

Assumption  (3.1)  used  in  (3.11)  yields  conclusions  (3.2),  (3.3)  and  (3.6). 

Beturning  to  (V*)  and  using  B’  e  L^CO,**),  v  e  L2(0,»:H),  P*  e  L2(0,*»jH)  gives 

conclusion  (3.4).  Combining  (3.3),  (3.4)  with  ^  *(u(t))  -  (v(t),  (t))  yields 

flt  ot 

4r  P(u(t))  «  (0,**) ,  and  this  together  with  assumption  (3.1)  implies  that  lim?(u(t)) 

exists.  To  establish  all  of  (3.7)  we  use  the  definition  of  subdifferentials  for  every 
T  i  3*(u)  and  for  every  w  e  H  *(u(t))  «  #(w)  +  (v(t),u(t)  -  w) ,  0  <  t  <  -.  Since 
»*  l  (0,**jH)  and  v  e  L2(0,«»h)  there  exists  a  sequence  {t  }  ♦  •  as  n  ♦  •  such  that 

II 

*T(t|i),u(t#)  -  w)  ♦  0  as  n  •*  "j  this  proves  (3.7),  and  from  it  easily  (3.8).  To  prove 
(3.5)  take  T  <  t  and  use  (3.4)  and  Cauchy-Schvarz  obtaining: 

l«(t)  -  «(T)|  <  /  |f|  (s)  Ids  <  /t  -  T  {/  (s)|2ds}1/2 

T  0 

«  x/t  -  t  (0  <  T  <  t  <  •)  . 
this  completes  the  proof  of  Theorem  3.1. 


-46' 


T 

>  h 


Proof  of  Theorem  3.2.  As  remarked  in  the  paragraph  preceding  Theorem  3.2  the 

coercivity  condition  (2.4)  implies  that 

Inf  #(*>  >  p(0)  >  —  , 

X«H 

so  that  assumption  (3.1)  is  satisfied.  Thus  conclusions  (3.2)~(3.5)  follow  immediately 
from  Theorem  3.1.  Zn  view  of  (3.5)  the  conclusion  u(t)  ♦  0  strongly  as  t  ♦  •  follows 
once  it  is  shown  that  u  c  L2(0,*;H).  But  using  assumption  (2.4)  with  c  >  0  and 
VC  l2(0,“;H)  for  v  c  3#(u)  (proved  in  (3.3))  one  has 

T  T  T  T 

c  /  |u(t)|2dt  <  /  (v(t)  ,u(t))dt  <  f  /  |u(t)I2dt  +  ~~  /  |v(t)|2dt  . 

0  0  *  0  '•  *  0 


T  m 

I  /  |u(t)|2dt  <  ~  /  |v(t)|2dt  :  -  (0  <  T  <  -) 


Since  t  >  0,  this  completes  the  proof  of  Theorem  3.2. 


4.-  Application  to  a  Problem  of  Heat  Plow  in  a  Material  with  Memory..  Zn  this  section  we 
study  the  heat  flow  problem  (1.4)  in  one  space  dimension  formulated  in  Chapter  1.  Section 
2.  We  use  the  existence  and  uniqueness  theory  of  Chapter  2  and  the  boundedness  and 
asymptotic  theory  of  Chapter  3.  Section  3.  to  deduce  the  principal  result.  Theorem  4.5 


below. 


For  clarity  of  exposition  we  restate  the  heat  flow  problem  (1.4),  Chapter  1,  Section 


(fr  (b.u  ♦  p*ul  -  en«<u  L,  “  Y*®(«  >_  ♦  h  (0  <  t  <  0  <  x  <  1) 

K  0  °  *  *  x  * 

u(0,x)  -  ^(x)  (0  <  x  <  1),  u(t.O)  -  u(t,1)  SO  (t  >  u)  , 

where  subscripts  denote  differentiation  with  respect  to  x.  We  assume  that  the  conditions 
(which  were  motivated  in  Chapter  1)t 


bQ  +  Re  S(ln)  >0  (!)  <  R)  , 


Cft  -  /  Y(t)dT  >  0  , 


where  6(ln)  ■  /  B(t)exp(-int)dt,  ere  satisfied.  We  also  assume  that  the  function 

0  , 

C  t  R  ♦  R  satisfies  assumptions  (o)  of  Unna  1.3  (Chapter  1),  uQ  c  Ho(0,1),  and  that 

12  2  2 

the  external  heat  supply  h<  L  (0,<*jL  (0.1))  n  l  (0,»it  (0,1)).  Under  these  assumptions 
ve  have  seen  in  Chapter  1  that  the  initial-boundary  value  problem  (4.1)  is  equivalent  to 


the  abstract  Volterra  equation 


S  +  b*  Au  t  f 


(0  <  t  <  -)  t 


in  the  present  application 


b(t)  -  *  (p*C)(t) 

b0 


f(t)  -  g-~~  +  (P*G)  ( t,  • ) 

bo 


(0  <  t  <  -)  , 


(0  <  t  <  «)  , 


where  p  is  the  unique  solution  of  the  resolvent  equation 


bflp(t)  ♦  (B*p)  (t)  -  - 
0  b0 


(0  <  t  <  -)  , 


C(t>  -  C  -  /  Y(T)dT  (0  <  t  <  -)  , 

0  0 


C(t,x)  “  b  x  (x)  +  /  h(t,x)dx  (0  <  t  <  •»  0  <  x  <  1)  , 

0  0 


end  the  nonlinear  operator  A  ■  3*  satisfies  assumption  (H^)  with  <P  given  by  the  proper, 
convex,  l.s.c.  function  defined  in  and  satisfying  the  properties  of  teams  1.3  in  Chapter 
1.  We  recall  that  in  two  or  three  space  dimensions  the  corresponding  heat  flow  problem  in 
a  bounded  domain  0  with  smooth  boundary  T  also  satisfies  the  Vblterra  equation  (V) 
with  b  and  f  as  above,  but  with  L3(0,1)  replaced  by  L*(fl) ,  and  with  A  •  i<f 
where  <9  is  the  proper  convex,  l.s.e.  function  defined  in  and  satisfying  the  properties  of 
Remark  1.4  in  Chapter  1.  We  recall  also  that  the  key  properties  of  the  kernel  b  and 


i 


of  the  forcing  tern  f  for  the  heat  flow  problem  are  stated  in  Lemmas  1.1  and  1.2  of 
Chapter  1.  With  these  properties  in  mind  all  that  is  needed  in  order  to  apply  Theorems  3.1 
and  3.2  of  Chapter  3  to  the  problem  under  study  is  to  show  that  the  frequency  domain 
assumption  (F)  of  Lemma  2.2  can  be  satisfied  for  physically  reasonable  classes  of 
relaxation  functions  6,  Y«  In  this  direction  we  have: 

Lemma  4.1.  Let  bg,  Cq,  8,  Y  satisfy  the  assumptions  of  Lemma  1.1,  Chapter  1.  Define 
the  kernel  b  _in  (V)  by  equation  (b) .  Then  the  frecruency  domain  assumption  (F)  of  Lemma 
2.2  is  equivalent  to  the  condition:  there  exists  6  >  0  such  that 

(c„  -  Rey(in))(b  ♦  ReB(in))  -  ImY(in)lmfl( in) 

(4.2)  Inf  — - - - : - - -  >  «  . 

<n«R)  |b0  ♦  B(in)l 

Proof  of  Lemma  4.1.  Define  the  constant  b^  >  0  and  the  function  B  as  in  Lemma 
1.1,  Chapter  1.  Taking  the  Laplace  transform  of  B  one  computes 

-  i  ,e0  "  VCin)  . 

B(in)  -  -  i  (— - s - bj  (n  «  R)  . 

n  bQ  +  8(in) 

Thus 

.  c  -  y(in) 

b  -  n  la  B(in)  •  Re( - s - )  (n  e  R)  » 

b0  ♦  (S(in) 

from  which  the  condition  (4.2)  is  an  immediate  consequence. 

Using  Lemma  4.1  one  can  construct  a  large  number  of  examples  of  functions  8  and  Y 

such  that  assumption  (F)  is  satisfied.  In  particular  one  has  the  following  physically 

Important  special  cases.  Rote  that  in  Corollaries  4.2  and  4.3  below  the  physical 
t  t 

conditions  b  ♦  /  6<T)di  >0  (0  <  t  <•)»«-/  y(T)dT  >0  (0  <  t  <  •)  are  both 

0  0  0 

satisfied  (although  they  are  not  explicitly  needed  in  the  theory),  because  the  functions 
$  and  Y  are  positive,  and  assumption  (y)  is  assumed  to  hold.  For  a  different  example 
In  which  (y)  is  satisfied  but  the  above  physical  conditions  need  not  hold  see  Remark  4.8 
below. 
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Corollary  4.2.  Let  b^  >  0,  >  0  and  8,  Y»  tB,  tY  r  L^(0,**).  Also  assume  that  B  and 

T  are  positive,  nonlncreasing  and  convex  on  (0,“),  and  that  the  assumption  (Y)  Is 
satisfied.  Then  assumption  (F )  is  satisfied  if  either  for  a  fixed  >  0  the  constant 
«g  >  0  Is  chosen  sufficiently  larcre,  or.  If  for  a  fixed  >  0  the  constant  b^  >  0  is 
sufficiently  large. 

Remark  4,3.  (i)  If  B  “  Y  =  0  (the  standard  heat  flow  problem)  (  F)  is  satisfied  for  any 

b 

choice  of  b  >  0 ,  e  >  0  with  5  »  — . 

0  0  c. 

0 

(ii)  If  8  =  0  and  y  satisfies  the  assumptions  of  Corollary  4.2,  (F)  is  satisfied 


c0  "  /  T(t)dt 

for  any  choice  of  b  >  0,  c.  >  0  with  S  «  - - - . 

0  0  b. 

o 

(ill)  If  Y  =  0  and  8  satisfies  the  assumptions  of  Corollary  4.2,  (F)  is  satisfied 

for  any  choice  of  b  >  0,  c_  >  0. 

0  0 

Sketch  of  Proof  of  Corollary  4.2.  The  proof  will  make  use  of  Lemma  4.1;  we  establish 

1  » 

(4.2).  Since  B.Y  (  L  (0,*)  and  are  positive,  nonincreasing  and  convex.  Re  B(ih)  and 

Re  Y(ln)  are  nonnegative.  The  function 


*  a  m 

ImY(in)Im$(in)  -  /  Y(t)sin  ntdt  /  B(t)sin  ntdt  (n  e  R) 

0  0 

is  even,  continuous,  zero  when  n  "  0,  nonnegative,  and  has  limit  zero  as  n  ♦  • 
(Riemann-Lebesgue  lemma).  The  denominator  in  (4.2)  satisfies 


o  <  bg  4  |b  +  B(in)|2  <  2b2  ♦  3 </  B(t)dt)2  (n  e  R)  . 

v  0 

Moreover, 

m 

♦  ReB(in)  >  bg  >  0  (n  e  R)  , 

(so  that  (PR)  is  satisfied),  and 


ReY(ih)  >  cn 


/  Y(t)dt  >  0 

0 


(•)  «'  R)  . 
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Therefore,  the  existence  of  6  >  0  such  that  (4.2)  holds  Is  established  for  choices  of 
bQ  and  Cq  as  asserted.  This  completes  the  proof. 

Another  physically  Important  case  for  the  heat  flow  problem  Is  the  following  special 
case  of  Lemma  4.1  and  Corollary  4.2. 

Corollary  4.4.  Let 


(4.3) 


n 

S(t)  -  £  b  , 
k-1 


(0  <  t  <  •)  , 

(0  <  t  <  -) 


with  bfc  >  0,  8^  >  0,  ck  *  0'  \  *  0  and  strict  inequalities  hold  for  at  least  one  pair 

a  c. 


b|^,  6^  and  one  pair  ck,  Let  >  0,  eQ  >  0,  and  cQ  -  J  ^  >  0. 


Then  the 


k-1  'k 


frequency  domain  condition  (F)  13  satisfied  if 


■  cl  _  n  b  *c. 

■>,('.  -  I  f)  >  1 1  £  5*)!  I  A 

le-1  Tk  k-1  pk  k-1  'k 


The  proof  of  Corollary  4.4  Is  a  consequence  of  showing  that  there  exists  a  4  >  0 
such  that  (4.2)  holds.  The  inequality  (4.4)  follows  by  using  elementary  calculus  to  find 
the  infimuK  over  n  «  R  of  the  expression  in  (4.2)s 


(s-  I 


o'.  n  1 
I  7 


K  .2 


,  s  cl  n  b, 

t)  -  "*(  Z  y-SH  Z  y— ) 

!\  k-1  ♦  n  k-1  e*  ♦  «r 


B  <  n  on 

l‘„  4  z  T-S  4  <  z  — '—J ) 

v  i_<  a-  a  i_<  a*  «* 


k-1  6‘  ♦  n 


k-i  ♦  n 


Mo  claim  is  made  that  the  constant  j  in  (4.14)  is  optimal. 

Me  next  combine  the  properties  stated  in  Lemmas  1.1,  1.2  and  1.3  of  Chapter  1  and 
Lemma  4.1,  and  Corollaries  4.2  and  4.3  above  with  the  abstract  theory  to  establish  the 
following  result  for  the  physical  heat  flow  problem  (4.1)  in  a  one-dimensional  material 


•SI. 


with  memory.  To  see  that  a  more  general  result  (not  necessarily  physical)  with  0  and  y 
oscillatory  can  hold  we  refer  to  Remark  4.8  below. 

Theorem  4.5.  Let  b^  >  0,  cfl  >  0,  let  0,  y,  t8,  ty  e  l/(0,— )  and  let  0  e  L2(0,— ). 
Assume  that  0  and  y  are  positive,  nondccrcasing,  convex,  and  that 

<y)  eo  “  /  >  0  • 

0 

Assume  that  0  i  R  ♦  R  satisfies  assumptions  (a)  of.  Chapter  1,  that  the  initial 

temperature  ufl  e  hJ(0,1),  and  that  the  external  heat  supply  h  «  L^O," ;H)  Ft  l2(0," jH) , 

2 

where  H  «  L  ( 0 , 1 ) .  Then  the  heat  flow  problem  (4.1)  has  a  unique  strong  solution  u  on 

d  sy  2 

[0.-)  x  (0,1)  such  that  -r—  t  L.  (0,-;H).  Moreover,  if  either  for  a  fixed  b_  >  0.  the 

"  IOC  1  '  ■  "  -----  -  Q  — - 

constant  >  0  is  sufficiently  large,  or  for  a  fixed  c^  >  0,  the  constant  3  0  is 
sufficiently  large,  then  the  solution  u  has  the  properties; 

U  t  l"(0,-;H)  H  l2(0,«;H),  f  L2<o,-;H)  , 

and  lim  u(t)  ”  0  strongly  in  H. 

t-*" 

Remark  4.6.  For  heat  flow  in  more  than  one  space  dimension  let  Q  be  a  bounded  body  in 

R?  or  R^  with  smooth  boundary  T.  Then  the  temperature  u  satisfies  an  equation  of  the 

form  (4.1)  with  the  operator  -a(u^)^  replaced  by  -V^fXIVul )Vu)j  the  boundary  condition 

is  u(t,x)  »  0  (0  «  t  <  ",  x  «  T),  ant  the  initial  condition  is  u(0,x)  -  uQ(x) 

(x  e  ft).  If  H  is  the  Hilbert  space  L2(8),  if  the  function  X  «  R  ♦  R*  satisfies 

1 

assumptions  (X)  of.  Chapter  1,  if  uQ(x)  «  HQ(0),  and  h  e  L  (0,-;H)  Ft  L  (0,«jH) ,  then 
the  results  of  Theorem  4.S  holds,  provided  the  constants  bQ  >  0,  cQ  >  0  and  the 
relaxation  functions  0  and  y  satisfy  the  assumptions  stated  in  Theorem  4.5. 

Remark  4.7.  let 


a 

-0kt 

0Ct)  -  I 

V 

k-1 

VI 

“V 

T(t)  -  [ 

V 

k-1 
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(0  <  t  <  -) 


where  b^  >  0,  >  0,  >  0,  >  0  end  strict  Inequalities  hold  for  at  least  one  pair 

M  Cv 


b.,  0.  and  one  pair  c  ,  y.  .  Let  b-  >  0,  cft  >  0,  and  c. 


S  '■fc 

~  I  r=  >  0. 

k-1  Tk 


Let  o,  uQ,  h 


satisfy  the  assumptions  of  Theorem  4.S.  Then  by  Corollary  4*4  all  conclusions  of  Theorem 
4. 5. hold  if  the  inequality  (4.4)  is  satisfied. 

Proof  of  Theorem  4.5.  Under  the  assumptions  of  the  theorem  the  heat  flow  problem 
(4.1)  is  equivalent  to  the  abstract  Volterra  equation  (V)  with  the  kernel  b  given  by 
equation  (b),  the  forcing  term  given  by  equation  (f),  and  the  operator  A  «  3 <f  where 
<f  t  H  ♦  is  the  proper,  convex,  l.s.c.  function  defined  in  Letitna  1.3,  Chapter  1  (or 

Re  mark  1.4,  Chapter  1  in  more  than  one  space  dimension).  To  establish  the  existence  and 
uniqueness  of  a  strong  solution  of  (V)  (equivalent  to  (4.3)),  we  apply  Theorem  3.10, 

Chapter  2.  Lemma  4.2,  Chapter  1,  shows  that  the  assumptions  of  Theorem  3.10,  Chapter  2 
concerning  f  are  satisfied  with  f(0,x)  "  uQ(x)  «  H^(0,1)  “  D(o)  .  Moreover,  (H^)  is 
satisfied.  Lemma  1.1,  Chapter  1,  shows  that  assumptions  (Hb)  are  satisfied.  Thus  to  apply 
Theorem  3.10,  Chapter  2,  we  must  still  verify  that  B*  e  [0 ,") .  From  the  expression 

for  B  in  Lemma  1.1,  Chapter  1,  we  compute 

(4.5)  B* (t)  -  -  *  <Y*PHt>  (0  <  t  <  •)  . 

J»0  0 

Since  0  is  monotone  by  hypothesis,  the  resolvent  equation  (p)  and  a  standard  argument 
(see  e.g.  Bellman  and  Cooke  (12])  show  that  p  is  monotone.  Finally,  since  y  is 
monotone,  it  follows  that  B*  e  BV[0,“).  Thus  Theorem  3.10,  Chapter  2,  yields  the 
existence  and  uniqueness  of  a  strong  solution  u  of  (V)  on  (0,*)  such  that 
u*  *  L^tO.-jH). 

Ve  shall  next  apply  Theorem  3.1  of  this  chapter.  Concerning  the  kernel  b  Lemma  1.1, 
Chapter  1,  shows  that  assumptions  (Hfe)  are  satisfied  with  b„  >  0.  Moreover,  Corollary 
4.2  shows  that  b  satisfies  the  frequency  domain  condition  (F)  if  bg  and  eg  are  chosen 
as  in  the  statement  of  Theorem  4.S. 

Lemma  1.3,  Chapter  1  (or  Remark  1.4,  Chapter  1  in  the  case  of  more  than  one  dimension) 

shows  that  *(y)  >0  (yen),  lim  *(y)  •  *•,  and  that  the  inclusion  3«(w)  »  0  has 

hrl*" 

w  •  0  as  the  only  solution.  Lamas  1.2,  Chapter  1,  shows  that  assumptions  (Hf)  are 


satisfied.  Therefore,  by  Theorem  3.1  the  solution  u  has  the  properties* 


•up  *(u(tn  <  ",  sup  |u(t)|  <  •»,  —  «  L  (0,*»).  u(t)  Is  uniformly 
0<t<**  0<t<- 

continuous  on  (0,*»). 

lim  *(u(t))  -  0,  and  u(t)  - k0  as  t  ♦  m  « 

t-*- 

Lemma  1.3,  Chapter  1  also  shows  that  under  assumption  (o)  the  coercivity  assumption 
(2.4)  is  satisfied  for  every  T  >  0  with  e  -  pQ»2  >  0  (or  another  positive  constant  in 
the  case  of  more  than  one  space  dimension  -  see  Remark  1.4,  Chapter  1).  Therefore,  by 
Theorem  3.2  one  also  has  u*  L2(0,~»H)  and  u(t)  ♦  0  as  t  ♦  •  strongly  in  B.  This 
completes  the  proof. 

Remark  4.8.  Suppose 

”®1* 

0(t)  “  b^e  coslt  (^,(1  >  0,  0  <  t  <  «■) 

"V 

T(t)  -  e^e  cosut  **1*^1  *  °'  0  *  t  C  ' 

and  assume  that  bg  >  0,  cQ  >  0.  Also  suppose  that  a  and  h  satisfy  the  assumptions  of 
Theorem  4.S.  Although  the  assumptions  concerning  f  in  Theorem  4.5  are  not  satisfied, 
one  still  has  by  Lemma  1.1,  Chapter  1,  that  assumptions  (Hb)  hold  with  b—  >  0  provided 

Vi 

<tl  2>°- 

Moreover,  B*  <  SV  (0,*)  from  (4.5),  and  the  existence  and  uniqueness  of  a  strong 

IOC 

solution  of  (V)  (equivalent  to  (4.1))  such  that  u*  (  L|oe(0,*;H)  follows  from  Theorem 
3.10,  Chapter  2.  Thus  to  obtain  all  of  the  conclusions  of  Theorem  4.5  we  need  only  verify 
that  the  frequency  domain  assumption  (F)  holds  in  order  to  apply  Theorems  3.1  and  3.2.  By 
Leans  4.1  it  suffices  to  find  sufficient  conditions  on  the  constants  b„,  eg,  b,,  0^,  X, 

“  so  that  (4.2)  holds.  An  elementary,  but  tedious  calculation  shows  that  A  >  0 


•ueh  that  (4.2)  holds  exists  in  the  esse  y*  ><*)*,  0*  >  X* ,  provided  assumption  (y)  above 
holds,  and  provided  >  0  is  chosen  sufficiently  large. 

While  no  claim  is  made  here  that  the  above  functions  0  and  Y  represent  physically 
plausible  relaxation  functions,  it  is  of  some  interest  that  the  theory  can  still  be 
applied.  In  this  connection  it  may  also  be  noted  that  here  the  function 


*  Vi  -V  v  -v 

C( t)  •  c.  *  /  Y(t)dx  -  c.  -  — - -  (1  -  e  co siot)  -  — - r  e  sinut  . 

0  Y1  +  “  Y,  -►  « 


In  a  genuinely  physical  problem  as  motivated  above  one  would  need  to  require 
C(t)  >0  (0  <  t  <  •),  as  veil  as  assumption  (Y)  •  However,  in  the  application  of  the 
theory  the  physical  requirement  C<t)  >  0  (0  <  t  <  •)  is  not  used  and  indeed,  for 

axa mple. 


Vi 


T1  4 


2  C1 


could  be  negative,  even  though  C(")  •  cg 


Vi 


>  0  holds 


Chapter  4 

Existence  and  Asymptotic  Behaviour  of  Positive  Solutions 
of  Monllnear  Velterra  Equation  for  Heat  Flow 


1.  Introduction.  In  thia  chapter  we  discuss  the  positivity  of  solutions  and  their 
asymptotic  behaviour  as  t  ♦  •,  of  the  nonlinear  Volterra  equation 
(V)  u(t)  ♦  (b*Au)(t)  »  f(t)  (0  <  t  <  -) 

in  the  general  setting:  b  :  [0,*»)  ♦  R  is  a  given  kernel,  A  is  a  nonlinear  (possibly 
multivalued)  m-accretive  operator  defined  on  a  real  Banach  space  X,  f  s  [<>,•)  *  X  is  a 
given  function;  the  integral  in  (V)  is  understood  in  the  sense  of  Bochner.  The  assumptions 
which  are  imposed  on  b,  A,  t  are  motivated  by  the  problem  of  nonlinear  heat  flow  in  a 
material  with  memory  formulated  in  Chapter  1  to  which  the  general  positivity  and  asymptotic 
theory  developed  in  Section  3  will  be  applied  in  Section  4.  A  different  application  of 
9*B*r*l  theory  is  given  in  Example  3.4,  Section  3,  to  a  nonlinear  conservation  law  with 
memory. 

This  chapter  which  generalises  and  complements  earlier  work  of  Clement  and  Nbhel  [18) 
on  positivity  and  of  Clement  [17]  on  limiting  behaviour  of  positive  solutions  of  (V)  is 
primarily  based  on  a  forthcoming  paper  by  Clement  and  Nohel  [19].  The  generalization 
enables  us  to  apply  the  theory  to  the  physical  problem  in  Section  4.  General  existence, 
uniqueness  and  continuous  dependence  results  for  solutions  of  (V)  which  need  not  be 
positive  have  been  established  by  Crandall  and  Nohel  [26]  and  by  Gripenberg  [34]  (see  also 
Theorems  3.5  and  3.7,  Chapter  3);  these  will  be  referred  to  as  needed. 

*•  will  motivate  the  assumptions  on  the  kernel  b  which  are  needed  for  positivity  of 
solutions  of  (V)  and  which  will  be  needed  throughout  the  analysis  by  means  of  a  simple 
linear  problem  at  the  end  of  this  section.  These  considerations  suggest  the  concept  of 
complete  positivity  of  the  kernel  b  (Definition  1.1  below)  which  plays  an  important  role 
in  the  analysis.  Some  properties  and  a  useful  characterization  of  completely  positive 
karnels  are  obtained  in  Section  2. 
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We  shall  consider  equation  (V)  in  the  slightly  less  general  form 
<Vg)  u(t)  ( b*  Au )  ( t )  5  uQ  +  (b*g)(t)  (0  <  t  <  ®)  . 

We  assume  throughout  the  following  minimal  assumptions: 

\  b f  l’  (0,«), 

loc 

( H ^ \  A  m -accretive  in  a  real  Banach  space  X  , 

[  u.  e  D(A) ,  and  g  e  (0,«°;X)  . 

u  loc 

The  motivation  for  taking  f  =  Ug  +  b*g  in  (V)  is  given  in  Section  3  (see  argument  at  the 
beginning  of  Section  3  following  (V  ) ) .  The  main  results  of  this  chapter,  described  in 
Section  3,  give  a  rather  complete  description  of  the  asymptotic  behaviour  of  the  positive 
solutions  of  the  abstract  equation  (V)  as  t  ♦  including  a  priori  estimates  for  their 
rates  of  decay.  The  results  are  then  applied  to  the  physical  problem  in  Section  4. 

The  additional  assumption  we  shall  make  on  the  kernel  b  in  order  to  insure 
positivity  of  solutions  was  first  introduced  in  [18],  it  is  motivated  by  the  following 
remark.  If  b  3  1  then  (Vg)  reduces  to  the  evolution  equation 


(du  ,  , 

—  +  Au  3  g 

dt 

u(0)  -  uQ 

It  is  well-known  ['13]  that  if  the  resolvent  *  (I  +  AA)  1  of  A  maps  a  closed  convex 

cone  P  of  X  into  itself  for  every  A  >  0,  then  u(t)  e  P  for  all  t  >  0,  provided 

that  Ug  «  P  and  g(t)  «  P  a.e.  on  to,*).  Let  us  take  for  instance 

X  -  {u  r  c[a,b] |u(a)  =  u(b)  «  0} 

2 

equipped  with  the  supremun  norm;  D(A)  *  {u  e  X  |  u  e  C  [a,b]  and  uxx  e  x}  and 
Au(x)  -  for  u  *  D( A) .  It  is  standard  that  A  is  m-accretive  in  X.  Moreover, 

if  P  «  {u  f  X|u(x)  >  0  x  t  [a,b]},  then  J^P  C  P  for  every  A  >  0;  thus,  as  is 
classical,  the  solution  o':  the  heat  equation  is  nonnegative  provided  that  the  initial 
value  Ug  and  the  forcing  form  g  are  nonnegative. 

We  want  to  consider  a  class  of  kernels  b  under  which  the  solution  (V),  resp.  (V^) , 
preserves  this  positivity  property.  This  requirement  is  useful  and  natural  in  the 


application  to  tho  mod*?!  of  boat  flow  in  a  material  with  memory  discussed  in  Section  4,  an 
in  Example  3.4  of  Section  3. 

Consider  (Vg)  with  Au  “  -uxx  with  D(A)  as  in  the  above  example.  It  is  easy  to 

give  necessary  conditions  to  be  imposed  on  b  in  order  that  positivity  is  preserved  by 

(V  )  whenever  uA  and  er  are  positive.  Let  X  denote  the  principal  eigenvalue  and  u 
Q  U 

the  corresponding  principal  eigenfunction  of  A,  normalized  by  max  u(x)  =  1.  Clearl 

•n  2  it  _  xc[a,b] 

X  -  (- - )  and  u(x)  =  sin(—  )(x-a).  If  u_  ~  au,  g(t)  =  A0(t)u  with  a  >  0,  and 

' b-a  b~a  0 

6(t)  >  0  where  £  €  L?  (0,°°),  then,  as  can  be  verified  directly,  the  strong  solution  o 
loc 

^g)  iS 

(1.1)  u(t)  =  (os(Xb)(t)  +  ( S*r ( Xb) )  ( t )  ]  u  (0  <  t  <  “)  , 

where  the  functions  s(b)  and  r(b):[0,“)  ♦  R  are  respectively  solutions  of  the  linear 
Volterra  equations 

(s( b) )  s(b) (t)  +  ( b* s ( b ) ) ( t )  =  1  (0  <  t  <  ») 

( r( b) )  r(b) ( t)  +  (b*r(b))(t)  =  b(t)  (0  <  t  <  »). 

Recall  the  standard  fact  (see  e.g.  R.  K.  Miller  [63])  that  if  be  L^oc(0,“),  t*’° 

functions  s(b),  r(b)  are  uniquely  defined  and  s(b),  r(b)  e  Lioc  Moreover,  if 

F  e  hi  (0,«)  the  unique  solution  of  the  linear  Volterra  equation 
loc 

(1.2)  u(t)  +  (b*u)  ( t)  =  F(t)  (0  <  t  <  -) 
is  given  by 

(1.3)  u(t)  =  F(t)  -  ( r(  b)  *F)  ( t)  (0  <  t  <  »)  . 

In  particular,  taking  F  =  1  in  (1.2),  one  has 

t 

(1.4)  s(b) ( t)  »  1  -  /  r ( b) ( T ) dt  (0  <  t  <  »)  , 

0 

so  that  s(b)  is  absolutely  continuous  on  [0,«)  whenever  b  f  ^loc^'"*'  The  ^u,‘ct^on 
sib)  is  called  the  fundamental  solution  of  (1.2),  while  the  function  r(b)  is  called  the 


resolvent  kernel  associated  with  b 


Returning  to  the  solution  (1.1)  of  (V^)  with  Au  »  ,  X  >0*  u,  u^,  g  defined 

above,  we  note  that  u(x)  >  0  for  x  t  (a,b).  Thus  the  solution  u(t)  will  be 
nonnegative  for  every  a  >  0  and  for  every  $  t  ,*•) ,  $  >  0.  only  if  the  functions 

r(b),  s(b)  are  nonnegative  on  (0,-).  Moreover,  if  one  Imposes  the  requirement  that  the 
solution  (1.1)  of  (Vg)  should  be  nonnegative  and  independent  of  the  length  of  the 
interval  (a,b),  it  is  clear  that  both  of  the  functions  r(Xb)  and  s(Xb)  must  be 
nonnegative  for  every  X  >  0.  We  remark  that  these  latter  necessary  conditions  imposed  on 
the  kernel  b  have  been  shown  to  be  sufficient  to  guarantee  the  preservation  of  positivity 
by  the  solution  operator  of  the  nonlinear  equation  (V)  in  the  general  case  of  A  m- 
accretive  on  X  (see  (18,  Theorem  4.5]). 

The  above  considerations  suggest  the  following  concept  of  complete  positivity  of  the 
kernel  b: 

Definition  1.1  We  shall  say  that  the  kernel  b  is  completely  positive  on  (0,T] 
it  b  t  L^  (0,T)  and  if  the  functions  r(Xb)  and  s(Xb)  -  1  -  1*r(Xb)  are  nonneqatlve 
on  [0,T]  for  every  X  >  0. 

Some  known  sufficient  conditions  which  insure  the  complete  positivity  of  the  kernel 
b  on  (0,T]  ares 

(I)  b  t  l'(0,T)  is  nonnegative,  nonincreasing,  and  log  b  is  convex  (see  Miller 
[62],  Levin  (52],  Clement  and  Nohel  (18]). 

(II)  (special  case  of  (i))  b  »  l'(0,T)  and  b  is  completely  mono  tonic  on  (0,T) 
(see  Miller  (62)). 

2.  Completely  Positive  Kernels.  In  this  section  we  give  an  alternate  and  useful 
characterisation  of  completely  positive  kernels  (Theorem  2.2)  which  will  be  needed  for  the 
development  of  the  asymptotic  properties  of  positive  solutions  of  the  abstract  Volterra 
equation  (Vfl) .  For  this  purpose  we  consider  the  linear  scalar  Volterra  equation  (1.2)  in 
the  form 

•  ♦  2*1  ■  ^  ♦  b*g 
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(2.1) 


vhere  b  t  t'lO.T),  Ug  «  R,  g  t  l'(0,T),  and  T  >  0.  Zts  unique  solution  (see  (1.1), 
(1.4))  Is  given  by 

(2.2)  u(t)  -  UgS(b)  (t)  +  (r(b)*g)  (t)  (0  <  t  <  T)  . 

In  the  following  proposition  we  list  some  elementary  properties  of  completely  positive 
kernels  which  are  needed  in  the  sequel. 

Proposition  2.1.  Assume  that  b  is  completely  positive  on  [0,T]  for  some  T  >  0.  Thant 

1)  b  is  nonnegative  on  (0,T)  and  for  every  p  >  0,  s(pb)  is  nonneoatlve  and 
nonincreasing  on  [0,T]. 

2)  For  every  p  >  0,  r(pb)  is  itself  completely  positive  on  (0,T). 

Mext,  assume  b  is  completely  positive  on  [0,Tl  for  every  T  >  0.  Then; 

3)  lit  b  e  L*(0,«),  then  for  every  p  >  0 

m 

lia  s(pb)ft)  -  (1  +  p  /  b(t)dT)"1  , 

0 


lr(pb)l  -  (p  /  b(T)dr) (1  +  /  b(T)dT)”1  . 

L  (0,»)  0  0 

4)  Xf_  b  /  t'(fl,«|,  then  for  every  p  >  0 

11a  s(pb)(t)  »  0  and  lr(pb)l  »  1. 

t*»  b(0,«) 

5)  2£.  b  /  and  b  e  AC(0,»),  then  for  every  p  >  0,  r(pb)  e  C[0,»)  and 

11a  r(pb)(t)  -  0  . 

The  proof  of  Proposition  2.1  is  elementary  and  is  omitted;  for  details  see  (19). 

In  the  next  result  we  give  an  alternate  and  useful  characterization  of  completely 
positive  kernels  b.  Some  arguments  used  are  similar  to  those  of  (33]. 

Theorem  2.2.  Let  T  >  0,  b  «  I»'(0,T),  b  $  0.  Then  b  is  completely  positive  on  [0,T] 
if  and  only  if  there  exists  o  >  0  and  k  e  L*(0,T)  nonneoatlve  and  non  increasing 
satisfying 
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(2.3) 


cb(t)  ♦  k*b(t)  -  1 


t  «  10, T] 


Remarks;  (i)  It  follows  from  (2.3)  that  a  >  0  If  and  only  If  b  e  L**(0,T).  If  this  is 
the  ease  b  «  a  1  s(a  'k)  and  thus  b  t  AC[0,T).  Conversely  if  b  e  AC(0,T],  then 

•  b(0)  '  >  0.  Moreover,  observe  that  if  a  >  0,  then  k  e  BV(0,T)  (equivalently 
k(0*)  <  *)  if  and  only  if  b'  <  BV[0,T)). 

The  importance  of  the  remark  a  >  0,  (k  e  BV(0,T))  is  that  for  kernels  b  satisfying 
the  assumption: 

(■)  b  t  AC [0 ,T]  ,  b(0)  >  0  ,  b*  «  BV [0 ,T] 

the  existence  and  uniqueness  of  a  generalized  solution  u  e  C((0,T)»  D( A) )  of  the  abstract 
▼olterra  equation  (V)  has  been  established  by  Crandall  and  Mohel  (Chapter  2,  Theorem  3.5), 
whenever  the  operator  A  is  m-accretive,  and  f(0)  e  D( A) ,  f  e  W1,1(0,T;X).  For  the 

special  case  X  -  H  a  real  Hilbert  space  and  A  -  3<fi  we  refer  to  Chapter  2,  Theorem  3.10. 

Recently  Gripenberg  (34,  Theorem  2]  has  extended  this  result  to  the  case  of  kernels 
b  ■  bj  ♦  b^,  where  b,  satisfies  the  above  regularity  assumption  and  where  bj  «  L^O.T), 
bj  la  positive,  nonincreasing,  and  log  bj  is  convex  on  (0,T)  with  A  and  f  as  above 
(see  Chapter  2,  Theorem  3.7).  This  result  with  bf  :  0  and  A  linear  was  already 
•stabliahed  by  Clement  and  Nohel  (18).  These  more  general  completely  positive  kernels  b 
correspond  to  the  case  a  «  0.  The  problem  of  existence  of  generalized  solutions  of  (V) 
with  only  the  assumption  that  b  is  completely  positive  is  under  study  end  will  be  treated 
elsewhere. 

(ii)  It  follows  from  Theorem  2.2  and  Remark  (i)  that  if  b  is  completely  positive,  then 
b  need  not  be  noninereasing;  it  also  need  not  be  convex  and  a  fortiori  log  convex.  Choose 

•  -  1  and  k(t)  -  1  for  t  «  (0,1J  and  k(t)  -  0  for  t  >  1}  then  b  -  s(k)  is 
completely  positive.  But  as  shown  in  Levin  [52],  b*  ■  -r(k)  is  negative  on  some  interval 

[0,o)  with  a  «  (1,2)  and  positive  in  (a, 2).  Thus  b  is  not  nonincreasing  on 

[0,2].  Moreover,  assume  b  to  be  convex  on  [0,*).  Then  b  is  strictly  increasing  for 

t  >  a,  and  moreover,  lim  b(t)  -  «.  But  this  is  impossible,  since  b(t)  <1  as  seen 
from  (2.3)  and  the  fact  that  k,b  are  nonnagative  and  a  •  1.  Thus  b  is  not  convex. 
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(Ill)  if  b  1*  completely  positive  and  absolutely  continuous  on  [0,T],  then  It  follows 
from  (2.3)  that  b(t)  <  b(0)  for  t  t  [0,T]. 

(Iv)  It  follows  from  Theorem  2.2  that  if  b  t  L?  (0,«*),  b  Is  positive,  decreasing,  log 

loc 

b  Is  convex  and  b(0+)  ■  "  ,  then  the  linear  Vol terra  equation  of  the  first  kind 
(2.4)  k*b(t)  -  1  .  t  >  0 

possesses  a  unique  solution  k  t  lJo^( 0 ,“)  which  is  nonnegative  and  nonincreasing.  However, 
given  k  t  L*(0,T],  k  nonnegative  and  nonincreasing,  equation  (2.4)  may  not  have  a 
solution  in  L^O.T).  (Take  k(t)  i  1).  Thus  when  a  -  0,  equation  (2.3)  does  not 
provide  a  way  to  generate  completely  positive  kernels  which  are  not  absolutely  continuous 
on  (0,T). 

We  omit  the  quite  technical  proof  of  Theorem  2.2.  It  makes  repeated  use  of  the 

following  result  due  to  Levin  (52).  If  u  satisfies  u  +  b*u  «  f  with  b  e  L^(0,T),  b 

nonnegative  and  nonincreasing,  f  e  L1(0,T),  nonnegative,  nondecreasing,  then  u  is 
nonnegative  on  [0,T].  The  proof  of  Theorem  2.2  may  be  found  in  [19]. 

3.  Qualitative  properties  of  abstract  Volterra  equations  with  completely  positive 
kernels.  In  this  section  we  study  some  properties  of  generalized  solutions,  including 
positivity  and  the  asymptotic  behaviour  of  positive  solutions  as  t  ♦  ",  of  the  nonlinear 
abstract  Volterra  equation 

(Vj)  a  ♦  b* Au  »  «0  ♦  b*g  (t  >  0)  . 

Although  our  results  are  stated  for  generalized  solutions,  it  is  obvious  that  the  results 
bold  for  strong  solutions,  whenever  strong  solutions  are  shown  to  be  generalized  solutions 
(see  Remark  (i)  following  Theorem  3.1  below). 

The  Justification  for  taking  f  ■  Ug  ♦  b*g  in  (V)  is  as  follows.  Zf  b  satisfies 
assumption  (H)  (Section  2,  Remark  (i)  following  Theorem  2.2),  if  f  e  W1,1(0,T»X),  and  if 
f(0)  «  OUT.  then  there  exists  a  unique  u^  «  DU)  and  a  unique  g  t  L1(0,TjX)  such 
that 


0.1) 


fit)  -  ♦  (b*g)(t) 


(0  <  t  <  T) 


Indeed.  Ug  “  f(0)  (  D( A)  and  9  is  the  unique  solution  of  the  linear  equation 
(3.2)  b(0 )g( t)  +  <b’»g)(t>  -  f*(t)  (0  <  t  <  T)  . 

Conversely,  if  b  satisfies  assumption  (H)  and  u^f  0(A)  ,  g  t  l/(0,T»X),  then  f  given 
by  (3.1)  satisfies  f(0) (  0(A),  f  t  W1#1(0,T»X).  We  shall  make  the  following  general 
assumptions: 


MW 

(H) 


A  is  m-aecretive  in  X 


u0 1  d(A) 


9  *  Lloc 

b  is  completely  positive  on  [0,“)  . 


The  basic  preliminary  result  assuming  the  global  existence  of  solutions  of  (Vg)  under 
assumption  (H)  is  known: 

Theorem  3.1.  If  A,  Ug,  g  and  b  satisfy  assumption  (H)  then: 

(1)  ^f  Uj  and  u2  are  the  generalized  solutions  of  (Vg)  corresponding  to  the  data 
Ug  9^,  1  “  1.2.  then  the  following  estimate  holds: 

(3.3)  |Ul(t)  -  u2(t)l  <  lu0>1  -  u0  2l  +  ( b* I g1  -  92D(t)  t  >  0  a.e. 

(2)  _if  p  is  a  closed  convex  cone  in  X,  _if  J^(P)  C  P  for  every  X  >  0,  and  if 
Ug  «  P  and  g(t)  e  P  a.e.  on  [0,«),  then  u(t)  *  P  a.e.  on  t0,«)»  moreover,  if 
▼  -  u  «  P  implies  J^v  -  J^u  e  P  for  every  X  >  0,  u,v  e  X,  and  if  “gf2  ~  uo,1  *  p* 
9j(t)  -  gj(t)  €  P  a.e.  on  (0,«),  then  Uj,  i  “  1,2,  the  corresponding  generalized 
solutions  of  (Vg>  satisfy  u2(t)  -  Uj(t)  *•  P,  a.e.  on  (0,»). 

Remarks:  (i)  The  existence  of  a  generalised  solution  in  the  linear  case  under  the 
assumption  b  completely  positive  was  proved  in  [18].  In  the  nonlinear  case,  when 
b  *  AC [0,T] ,  b(0)  >  0  and  b  «  BV(0,T]  ,  or  when  b  <  t1(0,T),  b  is  positive, 
nonincreasing  and  log  b  is  convex  on  (0,T),  the  existence  of  generalised  solutions  of 
<v  follows  from  results  Crandall  and  Nohel  [26]  and  Gripenberg  [34] ,  already  discussed  in 
Chapter  2  and  in  Remark  (i)  following  Theorem  2.2.  Moreover,  if  more  regularity  is  assumed 
on  b  and  f,  then  (see  (26),  (34])  the  generalised  solution  is  also  a  strong  solution  of 
(V,). 
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(II)  Estimate  (3.3)  was  proved  in  HI). 

(III)  We  sketch  the  proof  of  the  positivity  result  asserted  in  (2)i  the  details  appear  in 
Clement  and  Kohel  (18];  the  last  assertion  in  (2)  can  be  established  in  the  same  way. 


Consider  the  approximating  equation  of  (V^)  resulting  from  replacing  the  operator  k 
by  its  Yosida  approximation 

Ax  -  -1  <1  -  Jx>.  Jx  -  (I  -  XA)”1! 

CVX)  Ux  ♦  b*\«x  -  uQ  ♦  b«g  (0  <  t  <  -)  . 

By  the  definition  of  Ax  equation  (V^)  is  the  same  as 

°A  4  I  b*uX  -  “o  4  b*9  4  X  b*JXuX  * 

Using  (2.2)  with  g  replaced  by  g  +  ■j  Jxux»  one  easily  checks  that  (Vx>  is  equivalent  to 
the  integral  equation 


‘V 

where 


and 


With 


«X  “  f  A  4  W  ' 

*x  *  b)«0  ♦  Xr(i  b)*g 

W‘>  ■  "i  »*W  • 


'i  •  *  h  ■ 


the  Integral  equation  (Vx)  is  equivalent  to  the  nonlinear  equation 


<v 


\  m  Vvx  4  V  • 


Let  0  <  T  <  •  be  arbitrary.  By  the  complete  positivity  of  b  and  by  the  hypothesis  of 
Theorem  3.1,  part  (2),  *x(t)  tv  a.e.  on  [0,T].  Noting  that  Wx  maps  t1(0,T;X)  into 
Itself  and  recalling  that  the  operator  Jx  is  a  contraction  (since  A  is  m-accretive) , 
one  prjves  (for  details  see  (18]  and  (26])  that  some  iterate  of  Wx  is  a  strict 

contraction  on  L^IO.TiX)  for  n  sufficiently  large.  Thus  the  Integral  equation  (vx)  has 
•  unique  solution  vx  c  L^(0,TrX)  given  in 

"  11a  for  any  vfl  e  L1(0,T;X)  . 
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Therefore,  the  Approximating  equation  (V^)  has  the  unique  solution  u^  ■  ♦  v^.  But 

f^  «  P  a.e.  on  10,TJ,  and  if  v0<t)  <  P  a.e.  on  (0,Tl,  the  complete  positivity  of 

b  end  the  assumption  J^(P)  C  **  for  every  X  >  0  Insure  that  W^(v0>(t)  *  P  a.e.  on 

(0 ,T]  and  the  same  holds  for  W^(v0)(t)  for  every  n.  Consequently  u^C t )  f  P  a.e.  on 

(0,T] ,  and  if  u  is  a  solution  of  (V  )  on  [0,Tl  such  that  u  “  weak  lira  u.  in 

9  X*0  A 

L^(0,T>X) ,  then  u(t)  t  P  a.e.  on  fO,Tj.  Since  T  >  0  is  arbitrary  this  completes  a 
sketch  of  the  proof. 

We  next  obtain  some  results  concerning  t\e  asymptotic  behaviour  of  solutions  of  (V^) 
as  t  ♦  ••  We  first  consider  the  case  b  *  t'(0,»). 

+*  r 

Theorem  3.2.  Let  A,  u„,  g,  b,  satisfy  the  general  assumptions  (H)  with  b  f  0  and 
b  «  L'(0,-). 

m  •* 

(1.)  Let  g  t  L  (0,*»>X)  and  assume  there  exists  g  e  X  such  that 

lira  lg(t)  -  g  I  *  0.  Let  u  be  the  generalized  solution  of  (V  )  and  define 
t«*  _  _  «• 

U  -  J  (u  +  bg  J,  where  b  ■  f  b(t)dt  >0.  Then  the  following  estimate,  which 
-0  0  „  ‘ 

implies  strong  cor.veroence  of  u(t)  to  u  as,  t  ♦  •,  holds: 


/  b(Tjdl 

lu(t)-u"l  <  - -  lu  -u™ '  +  (b*lg  -  g“l)(t)  (0  <  t  <  -)  . 

b  0 


(2.)  In  addition,  let  b  «  L  (0,»)  and  lira  b(t)  ■  0.  Let  g  «  g1  +  g2,  where 

t4*  1  P 

g^  satisfies  the  assumptions  of  g  in  part  1 ,  and  where  g^  e  L  (0,»>X)  +  L  (0,*»fX), 

P  «  (1,«).  Let  u  be  the  generalized  solution  of  <V  )  and  let  u  «  J_  (u  ♦  bg  ) . 

9  h 

the  following  estimate,  which  implies  strong  convereence  of  u(t)  to  u *  as  t  ♦  •» 


(3.5) 


/  b(s)ds 

lu( t)  *  u"l  <  - -  lu-  -  u"l  ♦  <b*lg  -  g“l)(t> 

b 

♦  (b*l92D(t)  (0  <  t  <  •)  , 


where  g  ■  lim  g. (t). 

t*«  „ 

Remark.  Part  2  Is  proved  below.  Part  1  of  Theorem  3.2  was  proved  in  Clement  til)  end  the 

proof  will  be  omitted;  it  uses  ideas  similar  to  the  proof  of  part  2. 

Next  we  consider  the  case  where  b  /  L*(0,“)  which  is  needed  for  the  application  In 

m 

Section  4.  In  order  to  establish  the  strong  convergence  of  u  to  u  as  t  ♦  •  we  shall 
require  that  the  nonlinear  operator  A  in  (V)  satisfy  a  rather  strong  coercivity 
condition. 

MS 

Theorem  3.3.  Let  A,  u0,  g,  and  b  satisfy  the  general  assumptions  (H)  with 

b  f  L*(0,«*)  and  suppose  A  is  coercive  in  the  sense  that  there  exists  u  >  0  for  which 

A  -  MX  is  accretive  in  X. 

1.  bet  g  be  in  l“{0,«;X)  and  let  g“  e  X  such  that  lim  lg(t)  -  g"l  «  0.  tat 
u  be  the  generalized  solution  of  (Vq>  and  let  u  be  the  unique  element  in  X  satisfying 

Aa~  a  g".  Then  the  following  estimate  which  Implies  strong  convergence  of  u(t)  to  u ™ 
as  t  ♦  holds; 


(3.6)  l«(t)  -  «*l  4  l  r(Wb)(T,dT  ,w0  - 

♦  •*1(r(wb)*lg  -  g" l)(t)  (0  <  t  <  •)  . 

2.  in  addition,  let  b  be  *C [0 ,«•]  and  g  ■  *  9j  *here  g^ ,  g2  satisfy  the 

assumptions  of  Theorem  3.2,  part  (2),  with  g.  -  lim  g  (t).  Let  a  be  the  generalised 

m  mm 

solution  of  (Vg),  and  let  u  be  the  unlgue  element  in  X  satisfying  Au  t  g^.  Then  the 

following  estimate,  which  implies  strong  convergence  of  u(t)  to  a  as  t  ♦  holds; 


Iu(t>  -  u" I  <  /  r(ub)(x)dx  lu.  -  u**l 

t 

♦  M_1Cr{wb)*lgt  -  g"l)(t) 

♦  M_1(r(wb)*lg2l)(t)  . 


Remarks .  (i)  Since  b  Is  completely  positive  and  b  f  L*(0,»),  it  follows  (see 

Proposition  2.1)  that  r(ub)  <  l'(0,“),  and  therefore,  if  the  assumptions  of  part  1  hold, 

(3.6)  implies  lim  lu(t)  -  u  I  ■  0. 
t—» 

When  b  also  satisfies  b  *  AC[0,*]  it  follows  (see  Proposition  2.1]  that 

r(«b)  «  L^O,-)  n  C(0,«),  and  lim  r(wb)(t)  -  0.  Therefore  (3.7)  implies 

lim  lu(t)  -  u"l  -  o,  if  the  assumptions  of  part  2  hold, 
t— 

(11)  As  is  clear  from  the  proofs,  the  assumption  g<  L  (0,“;X)  and  there  exists 

m  m 

g  such  that  lim  lg(t)  •  g  I  •  0  in  part  (1)  of  Theorem  3.2  can  be  weakened  to 
1  *— 

9 1  l  (0,»;X)  and  there  exists  9  *  X  such  that  lim  (b  •  Ig  -  9  l)(t)  •  0.  Similar 

loc  t—» 

generalizations  can  be  made  in  Theorem  3.2,  part  (2),  and  in  Theorem  3.3. 

Proof  of  Theorems  3.2,  part  2.  As  in  the  proof  of  Theorem  3.2,  part  (1)  in  [17],  we  first 

prove  the  result  with  A  replaced  by  A^,  X  >  0,  and  then  we  pass  to  the  limit  as 

14  0.  For  X  >  0,  let  u^  be  the  strong  solution  of  the  approximating  equation 

u^  ♦  b*A^u^  ■  ♦  b*g  (t  e  (0,«)>  . 

Using  the  definition  of  A^  and  applying  (2.2)  we  see  that  u^  satisfies  the  equation 
(3.U)  u^  -  r(X  'b)*J^u^  ♦  s(X  'b)u#  +  Xr(X  'b)*g 

for  t  *  tO,-J .  Since  A^  is  also  m-accretive,  there  is  a  unique  u^  satisfying  the 
limiting  equation 

(3.0)  ♦  bA.u?  ■  u.  +  bg™  . 

X  XX  0  1 

Using  the  fact  that  b  e  l\o,«)  and  9  ■  9j  ♦  g2  we  can  rewrite  (3.9)  in  the 
equivalent  form 

(3.10)  u"  ♦  b*A^*  "  uo  ♦  b*g  ♦  b*(g"  -  g^)  -  b*g2  -  {w~  , 
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where 


<(t)  •  /  b(s)ds  and  w  A  ^  g 

Let  n  <  (0,**)  ♦  R  be  the  unique  solution  of  the  linear  equation 

!>  ♦  X^b^n  “  t  l 

then  obviously 

(3.11)  qw^  ♦  X  'b*!^  -  £w^  (0  <  t  <  •)  • 

Using  (3.10),  (3.11),  (2.2)  and  the  definition  of  A^  we  obtain 

(3  12)  -  rU^M-J^u"  ♦  sd1  b)uQ  +  Xr(X_1b)*g 

♦  Xr(X“1b)»(g"  -  g,}  +  Xr(X”1b)*g2  -  qw“  . 

Subtracting  (3.12)  from  (3.S)  we  obtain 

(3.13)  lu][(t)  -  u^l  <  (r(l-1b)*lu“  -  uxD(t)  +  X(r(X_1b)*lg"  -  9,0(1:) 

♦  X(r(X”1b)«lg2l(t)  ♦  |q|(t)  lw“l  . 

It  la  shown  in  Clement  (11)  (see  argument  following  (3.18)  in  [17])  that  n  >  0.  Thus  by 
•sing  the  same  argument  ax  in  (17]  one  gets  (take  convolution  of  (3.13)  with  X  *b)x 
13.14)  lux(t)  -  u"l  <  C(t)lw"l  ♦  (b*lg,  -  g"l)(t) 

♦  (b*lgaO(t)  (0  <  t  <  •)  . 
the  conclusion  (3.S)  folic ws  by  using  (3.9)  and  rewriting 


/  b(s)ds 

C(t)lw  "l  »  - -  lu  -  n*l  , 

/  b(s)ds 

0 


and  then  letting  X  ♦  0.  Mote  that 


Proof  of  Theorem  3.3.  We  first  establish  the  results  with  A  replaced  by 
■  tfl  ♦  B^,  X  >  0,  where  B^  Is  the  Yosida  approximation  of  B,  defined  by 
B  ■  A  -  ml.  Note  that  B  is  s^accretive  in  X.  Let  u^  be  the  strong  solution  of  the 
approximating  equation  to  (V^)  written  in  the  form: 

0*15)  g  +  wb*U,  +  Ub*w”'B,Ut  »  U.  ♦  (ilb*W  %  . 

X  X  X  X  0  ’ 

Since  the  kernel  b  /  L*(0,“),  we  transform  this  equation  into  a  form  which  has  the 
property  that  its  new  kernel  will  be  in  L^(0,**)  and  completely  positive.  Indeed,  if  we 
take  the  convolution  of  (3. 15)  by  r(ub),  subtract  the  result  from  (3.15),  and  use  the 
definition  of  r(ub)  we  get  the  approximating  equation  equivalent  to  (3.15): 

(3.16)  u^  ♦  r(wb)*w  'b^u^  ■  u^  +  r(ub)*(u  1g  -  uQ)  . 

From  Proposition  2.1,  r(ub)  is  completely  positive  and  r(ub)  e  L*(0,**),  with 

as 

/  r(ub)(5)di  -  X. 

0 

To  prove  Theorem  3.3,  part  1,  tie  wish  to  apply  Theorem  3.2,  part  1,  to  (3.16).  If 

g  satisfies  the  assumptions  of  Theorem  3.3,  part  1,  so  does  u  *g  -  u^.  Thus  all 

assumptions  of  Theorem  3.2,  part  1  are  satisfied  with  b  replaced  by  r(ub),  A  replaced 
-1  -1 

by  «  8^,  g  replaced  by  <■>  g  -  uQ,  and  u  by  u^.  We  obtain  (by  (3.4)): 


(3.17)  lu^(t)  ~  0*1  <  /  r(«b)(T)dx  lu^  -  u™l  ♦  m'"1(r(wb)*ig  -  g" l)(t)  (0  <  t  <  •)  , 

m 

where  u^  is  the  unique  solution  of  the  limiting  equation 

(3.18)  ®X  +  "  BXUX  ”  “  ' 

which  exists  because  B^  is  m-accretive  and  u  >  0.  Note  that  (3.17)  is  the  estimate 
(3.6)  with  u  replaced  by  u^. 

If  g  satisfies  the  assumptions  of  Theorem  3.3,  part  2,  and  b  <  AC(0,“),  if  follows 
from  Proposition  2.1.  that  r(ub)  <  L*(0,")  n  C[0,«*].  Thus  with  9  ■  g,  ♦  g2  we  can  apply 
Theorem  3.2,  part  2,  to  (3.16)  and  we  obtain  (from  (3.5))  the  estimate: 

lu.(t)  -  uTl  <  /  r(wb)(t)dT  lua  -  u7l 

(3.19)  *  *  t  0  a 

♦  (m'‘1r(«b)«l91  -  g*l)(t)  ♦  («“1r(wb)«lgJ[l)(t)  (0  <  t  <  -)  , 


-69- 


where  ie  the  unique  solution  of  the  Uniting  equation 

•  wl  m  -1  • 

0.20)  u^  ♦  w  B^u^  ■  «  91  . 

Mote  that  (3.19)  is  the  estimate  (3.7)  with  u  replaced  by  u^. 

Since  8  is  m-accretive  lim  u.  ■  u  ,  where  in  the  case  of  (3.18)  u  satisfies  the 

X40  * 

Uniting  equation 

•  m  • 

0.21)  u  +  u  Bu  a  w  g  « 

or  equivalently  u  satisfies  the  limiting  equation  ' 

(3.22)  Bu  a  g  . 

Similary.  in  the  case  of  (3.20)  we  find  that  u*  satisfies  the  limiting  equation 

(3.23)  Bu  i  g1  • 

Zt  remains  to  prove  that  lim  u.  -  u  in  L^(0,T;X)  for  every  T  >  0,  where  u  is 

X40 

the  generalized  solution  of  (V^) .  Having  done  so  it  is  immediate  that  the  estimates 

0.17).  (3.19)  hold  with  u.  replaced  by  u  thus  obtaining  (3.6)  and  (3.7).  He  know  that 

T  -  i 

lim  /  lu,  (t)  -  u( t) Idt  ”  0,  where  u.  t  L,  satisfies 

X40  0  X  *  lOC 

0.24)  •  u#  ♦  b#g  . 

Introduce  the  notation  6  «  X(1  +  X»)"\  Since  A  •  ml  +  8,  one  easily  checks  that 

A^  -  w(1  ♦  «X)_1I  ♦  Bd  . 

Me 

Thus  the  solution  u^  of  0.24)  satisfies  the  equation 

0.25)  ♦  b*Bj  u^  -  uQ  +  b*(g  -  »(1  ♦  uX)_,u^)  , 

or  by  (3.15)] 

0.26)  Hg  ♦  b*Bjug  •  «q  ♦  b*(g  -  «Uj)  . 

Vo  compare  solutions  of  (3.25)  and  (3.26)  we  apply  the  inequality  (3.3)  of  Theorem  3.1.  and 
we  obtain  (note  that  u01  ■  u02  ■  Ug) 

lnA  -  «jl  <  wb*{1  ♦  Xw)'1  iZk  -  u4»  ♦  XwJ(1  ♦  X»)"tb*luxl  , 
and  hence  also  the  estimate 

0.27)  l«x  -  «4I  «  wb* lux  -  Uj!  ♦  <w3b*luxl  . 

Xt  follows  from  0.27)  that  for  every  T  >  0 


(3.28)  /  lu  -  u,  I ( t)dt  <  «w2lbl  •  (l  ♦  Ir(-wb)  I  )  /  lu.(t)ldt  . 

0  I*  (0,T)  fc(0,T)  0 

T 

Since  u.  converges  to  u  In  L  (0,T)X),  we  finally  obtain  lim  /  lu.  -  ulftldt  -  0, 

*  X*0  0 

which  proves  that  the  u^  converges  in  L1  (0,T;X]  to  the  generalized  solution  u  of 

(Vg)  for  every  T  >  0.  This  completes  the  proof  of  Theorem  3.3. 

Example  3.4.  A  Conservation  Law  with  Memory.  As  an  illustration  of  Theorems  3.1  and  3.2 
we  consider  the  existence  and  qualitative  properties  of  positive  solution  of  the  problem 

t 

(e)  u(t,x)  ♦  /  b(t-s)$(u(s,x) )  ds  »  u  (x)  (t  >  0,  x  t  R)  . 

0  xu 

We  assume  that  $  *  C*  (R)  is  a  given  function.  If  b  =  1  problem  (c)  is  equivalent  to 
the  nonlinear  conservation  law  in  one  space  dimension: 

«t  ♦  8(u)x  “  0  *  u(0,x)  ■  uQ(x)  (x  e  R)  . 

Although  no  particular  physical  significance  is  claimed  for  (e),  it  evidently  contains  the 
usual  conservation  law  as  a  special  case.  The  latter  has  been  studied  extensively  from 
special  points  of  view.  Crandall  (24]  has  shown  that  if  $:R  ♦  R  is  a  given  smooth 
strictly  Increasing  function  (actually  $  continuous  is  sufficient)  such  that  8(0)  «  0, 
then  the  operator  A  defined  by  Au  ”  8(*Ox  on  the  Banach  space  X  “  L1 (R) ,  with 
0(A)  ■  (u  e  t1(R)s8(ux)  «  I^lR)}  (see  (24,  Definition  1.1  and  Theorem  1.1]),  is  m- 
accretive  on  X,  and  D( A)  ■  X.  Moreover,,  one  has:  J^(0)  "  0  and  J^u  4  j^v  (X  >  0), 
whenever  u  <  v,  u,  v  t  L1  (R) . 

Za  (c)  assume  that  b  r  b  •  0  completely  positive  on  (0,**),  and 

Up  *  D(A)i  to  be  specific  take  b  nonnegative,  nonincreasing  and  log  b  convex  on 
(0,m),  b  *  0.  Then  by  Gripenberg's  result  (34),  see  Theorem  3.7,  Chapter  2,  and  Remark 
(1)  following  Theorem  2.2,  and  by  Theorem  3.1,  problem  (e)  has  a  unique  generalized 
ablution  ui  u  is  nonnegative  whenever  Ug  is  nonnegative,  and  u^  >  u^  whenever 
Ugj  >  «pj.  If,  in  addition,  b  «  L1  «>,•),  then  this  generalised  solution  u  converges 


strongly  In  L*(R)  u  t  ♦  •  to  the  element  *  D( A)  which  is  the  unique  solution  of 
the  limit  equation 


«„(x)  ♦  (/  b(t)dt)  $ (u^( x) )  -  u  (x)  (x  «  R)  t 

0  x  o 


exists  and  is  uniquely  defined  since  /  b(t)dt  >  0  and  A  is  m-accretive. 

0 


4.  Application  to  Konllnear  Heat  Flow  in  a  Material  with  Memory.  In  this  section  we  apply 
the  theory  developed  in  Sections  2  and  3  of  this  chapter  to  discuss  the  global  existence , 
uniqueness,  positivity,  and  decay  of  positive  solutions  of  the  nonlinear  heat  flow  problem 
formulated  in  Chapter  1,  Section  2.  This  problem  was  discussed  from  a  different  point  of 
view  in  Chapter  3,  Section  4.  In  Chapter  3,  Section  4,  boundedness  and  decay  results  were 
obtained  without  consideration  of  the  positivity  of  the  data.  Here  the  principal  concern 
Js  the  existence  and  decay  of  positive  solutions  when  the  data  are  positive. 

Referring  to  Chapter  1,  Section  2,  we  again  study  the  heat  flow  initial-boundary  value 
problem  (see  (1.4)  Chapter  1,  also  (4.1)  Chapter  3)i 

Ij:  H>0  «(t,x)  ♦  (S*u) (t,x)]  -  cfl  o(ux(t,x)>x  -  (y*o(ux>x)(t,x) 

♦  h(t,x)  (0  <  t  <  •,  0  <  x  <  1) 


(4.1) 


«(t,0)  -  u(t,l)  SO  (t  >  0) 


l  u(0,x)  "  u#(x)  (0  <  x  <  1) 

As  before  we  assume  that  the  constants  b„,  Cq  are  positive,  that  0,  y  *  L1  (0 ,-) ,  and 
that 


(0) 


b#  ♦  /  0(t)dt  >  0 


(Y)  e  -  /  Y(t)dt  >  0  . 

0 

We  else  assume  that  the  function  o  s  R  ♦  R  satisfies  assumptions  (o)  of  Lemma  1.3, 

Chapter  1,  and  that  the  external  heat  supply  satisfies  (at  least)  the  assumption 
1  2 

h  e  L,  (0,“jL  (0,1)).  Defining  the  functions 
loc 

t 

(C)  C(t)  -  c  -  /  Y(T)dT  (0  <  t  <  •)  , 

0  0 

t 

(G)  G(t,x)  -  b  u  (x)  ♦  /  h(T,x)dt  (0  <  t  <  -,  0  <  x  <  1)  , 

00  0 

the  heat  flow  problem  (4.1)  is  equivalent  to  (see  Chapter  1,  Section  2)  the  abstract 
Volterra  equation 

(»t)  u  ♦  0*u  +  c*Au  -  G  (0  <  t  <  -,  0  <  x  <  1)  > 

here  we  have  taken  the  constant  b.  as  1  without  loss  of  generality.  The  nonlinear 
operator  A  -  3 *>,  defined  in  Lemma  1.3,  Chapter  1,  satisfies  the  properties  established  in 
lemma  1.3,  Chapter  1.  provided  o  satisfies  assumptions  (o),  and  if  u^e  H^(0,1),  then 
G(0,x)  »  uQ(x)  e  hJ'O,  1 )  ”  D (<P).  We  note  also  that  in  the  case  of  more  than  one  space 
dimension,  the  relevant  heat  flow  problem  formulated  in  Chapter  1,  Section  2  is  equivalent 
to  the  Volterra  equation  (V^)  with  the  nonlinear  operator  A  defined  in  Remark  1.4, 

Chapter  1;  the  theory  developed  below  applies  equally  well  to  this  case. 

The  Volterra  equation  (Vj)  may  be  written  in  the  standard  form  (V^)  of  Section  3  above 
by  defining  the  resolvent  kernel  r{0)  of  0  to  be  the  unique  solution  of  the  linear 
equation: 

(r(0))  r(0)  ♦  0*r(0)  -  0  (0  <  t  <  ->  » 

clearly,  if  0  «  L1(0,«),  then  r(0)  c  L  (0,«»)  (at  least).  Next,  define 

IOC 

b  t  IC,«)  by 

Cb)  b  -  C  -  r(0)*C  , 

where  C  is  the  function  defined  In  (C).  Then  the  variation  of  constants  formula  shows 


that  (V|)  is  equivalent  to  the  Volterra  equation 

u  ♦  b* Au  “  G  -  r(0)*G 


(0  <  t  <  «)  l 


taking  bQ  -  1  in  (G)  one  sees  that  (v1),  as  veil  as  the  heat  flow  problem  (4.1),  is 
equivalent  to  the  Volterra  equation 

(4.2)  u  ♦  b*Au  -  uQ  ♦  l*(h  -  r(0)*h  -  uor(0))  (0  <  t  <  -)  . 

The  result  of  applying  the  theory  of  Section  3  on  positivity  and  decay  of  solutions  of 
(Vg)  to  the  heat  flow  problem  (4.1)  is: 

Theorem  4.1.  Let  B  be  bounded,  nonneoatlve,  nonincreasinc  and  convex  on  [0,*»).  Let 
T  be  nonneoative,  nonlncreasing,  log  convex  and  bounded  on  [0»“).  Let  assumption  (y) 
hold,  and  let 

(4.3)  0* (t)  +  0(t)  <  0  a.e.  for  t  e  [0,-)  . 

C° 

Let  the  assumption  (a)  be  satisfied,  and  let  A  ■  3^  where  >P  is  defined  in  Lemma  1.3, 
Chapter  1. 

1.  If  uQ  «  L2(b,1)  and  if  the  forcing  function  he  L^^no,**)  x  [0,1)),  then  the 

nonlinear  Volterra  equation  (4.2)  (equivalent  to  the  heat  flow  problem  (4.1))  possesses  a 

2  2 

unique  strona  solution  u  on  (0,»),  such  that  t  u*  e  L.  (0.«:L  (0.1)):  if 
~  loc 

Vq  e  nj(0,1),  then  u*  e  L2oc(0,-;L2(0,1>.* 

2.  If  the  data  uQ  and  h  satisfy  uQ  ^(x)  <  u0  2(x)  a.e.  on  [0,1]  and  If 
h\t.x>  <  h2(t,x)  a.e.  on  [0,«)  x  [0,1]:  then  the  corresponding  strong  solutions  u^ 

(1  ■  1,2)  satisfy  u^t.x)  <  u2(t,x)  a.e.  [0,**)  x  [0,1]:  in  particular.  If  ufl(x)  >  0 
and  h(t,x)  >  0  a.e.  on  [0,1]  and  [0,*»)  x  [0,1)  respectively,  then  u(t,x)  >  0  a.e. 
on  [0,-)  x  [0,1]. 

3.  If,  In  addition.  0  e  L1(0,-),  and  if  h  -  h,  ♦  hj,  where  h.j  e  l"(0,»:L2(0,1>) 

•  2  m 

and  there  exists  h.  «  L  (0,1]  such  that  lin  h.(t)  -  h.  _  »  0  ,  and  where 

p  2  L2(0,1) 

h2  «  I.  (0,«;L  (0,1))  for  some  p  >  1,  then  the  strong  solution  u  of_  [4.2]  converges 

strongly  in  L2(0,1)  as  t  ♦  •  to  the  element  u  e  L2(0,1):  u  Is  the  unique  solution 


•  h1  v  *  “ 

9,  -  “  (l  ♦  Y  "  /  Y(t)dt,  C(-)  -  c0  -  /  Y(t)dt  >  0  . 

In  particular.  If  h^  «  0,  then  u  ■  0. 

Before  giving  the  proof  of  Theorem  4.1  we  state  a  lemma  which  establishes  some 
properties  of  the  kernel  b  defined  by  (b). 

Lemma  4.2.  Let  0,  '1,  C  satisfy  the  assumptions  of  Theorem  4.1.  Then  b  defined  by  (b) 
is  completely  positive  on  {0,«*),  b  satisfies  the  assumption  (H)  (Remark  (i)  following 
Theorem  2.2).  and  o,k  associated  with  b  in  Theorem  2.2  satisfy  a  “  c^*  >  0  and 
k  €  L1  (<>,•)  with 


(4.4) 


/  k(T)dt  -  [-£—  Y(1  + 

0  c0  1  1 


where  0  ■  /  0(t)dt,  Y  “  /  Y(t)dt.  Moreover ,  b  (  and  b*  t  L1  ((),«•). 

0.0 

Proof.  Since  the  functions  0  and  C  «  AC,  [0,»)  it  follows  that  the  functions  r(0) 

xoc 

and  b  t  ACloc(0,«)  (see  definitions  (r(0))  and  (b)  respectively).  Note  that 
b(0)  -  >  0.  Define  a  -  b(0)  1  and  let  k  be  the  solution  of  the  linear  Vclterra 


equation 

(k)  b(0)y  +  b'*y  “  (0  <  t  <  -)  . 

Since  b’  e  lJ  [(>,«•),  k  e  L?  (0,"),  and  since 
IOC  10c 


~  (b»k)(t)  -  b{0)k(t)  +  (b'«k)(t)  -  -  , 

one  has  by  integration  that  k  satisfies  the  linear  Volterra  equation 
(4.5)  eb(t)  ♦  (k*b)(t)  -  1  (0  <  t  <  •)  . 

Xn  order  to  show  that  b  is  completely  positive  it  suffices,  by  Theorem  2.2,  to  show 
that  k  Is  nonnegative,  nonincreasing  and  bounded  on  (0,*).  '  We  first  observe  that  the 
assumptions  made  on  Y  imply  that  C,  *C'  are  conve.c  and  log  (*C* )  is  convex  on 
(0,«) .  This  in  turn  Implies  that  log  C  is  convex  on  (0,*),  see  C.  Cripenbcrg  (33). 


■75 


Since  C  in  nonncgative,  nonincreasing  and  belongs  to  C  is  completely 

positive  on  |0,a).  Moreover  C  also  satisfies  assumption  (H)  (see  Restart  (i)  following 

theorem  2.2).  It  follows  from  Theorem  2.2  that  there  exists  a  >0.  end 

c 

k  <  (0,“),  nonnegative,  nonincreasing  and  bounded  satisfying 

c  loc 

(4.6)  a  C(t)  ♦  (k  ♦C)(t)  -  1  <0  <  t  <  •)  • 

c  c 

Note  that  a  »  c  '  -  b'1  ( 0 )  ■  a  .  From  the  definitions  of  b  in  (b),  of  r(0),  end  from 

C  0 

(1.2)  and  (1.3)  (see  Section  1,  Chapter  4)  it  follows  that 

(4.7)  C(t)  -  b(t)  «■  (B*b)(t)  (0  <  t  <  •)  . 

Substituting  (4.7)  into  (4.6)  yields 

ob  ♦  (k  ♦  aS  ♦  k  *8)*b  -  1  . 
c  c 

end  thus  (4.5)  implies  that 

(4.8)  k(t)  -  k  (t)  ♦  «8(t)  ♦  (k  *8){t)  (0  <  t  <  «)  . 

c  c 

Since  k  *  BV(0,»)  we  have 
c 

.  t 

5^  (08  ♦  k/8](t>  -  «8*(t)  +  ke(0)6(t)  +  /  8(t  -  t)dkc(T)  a.e. 

on  t0,«).  Hypothesis  (4.3)  and  the  identity  (4.6)  imply  that 

S0'(t)  +  k  (0)8(t)  -  —  I8'(t)  ♦  2121  0{t)j  (  o.  Moreover,  since  k  is  nonincreasir.g 
c  c _  c„ 

0  0 

sad  8  is  nonnegative 

t 

/  8(t  -  T)dk  (T)  <  0  (0  <  t  <  •)  . 

0  C 

Thus  k  is  nonnegative  and  nonincreasing  on  (0,») .  Therefore,  one  also  has  k  *  BV(0,“) 
if  k  «  X,  (0,1].  But  ke  and  8  era  bounded  and  8  «  lVo.I)  imply  k  «  lT[0,1]  (note 
that  here  the  assumption  8  *  L*(0,«)  is  not  needed).  Using  Theorem  2.2  again  it  follows 
that  b  is  completely  positive  and  satisfies  assumption  (K). 


1  *1 

Ve  next  establish  that  b  f  L  ((>,••)*  Since  C  *  -Y  *  L  and 

11*  C{t)  -  C(-)  -  c  -  /  Y(a)ds  >  0,  it  follows  that  C  f  L^O,**).  If  b  e  1/(0, «>,  it 
would  follow  from  (4.7)  and  the  assumption  8*  L  (0,«)  that  C  (  L  (0,«*),  a 
contradiction.  Thus  b  /  L^tO,*1). 

We  next  prove  that  b'  e  L1{0,“1.  Indeed,  from  (b)  it  follows  that 
b*(t)  -  C'(t)  -  C(0)r(B)(t)  -  (r(B)*C')(t)  . 

But  C*  »  -Y  «  L^O,**)?  moreover,  r(B)  «  L^O,**),  since  B  is  nonnegative, 
nonincreasing,  convex  and  B  *  l'(0,«)  (use  the  Paley-Wiener  theorem  and  the  fact  that  B 
is  positive  definite). 

Finally,  we  show  that  k  *  L^tO,-).  From  (4.8)  and  the  fact  that  B  «  L1  «>,•),  it  is 

sufficient  to  prove  k  <  L1(0,-).  From  (4.6)  and  the  fact  that  C  is  positive, 
c 

nonincreasing,  C(”)  >  0  and  k  is  nonnegative,  we  have 

t 

C(-)  /  k  (T)dT  <  C*k  <  1  (0  <  t  <  -)  , 

0  C  C 


which  proves  that  k  «  l'(0,*»).  Formula  (4.4)  follows  easily  from  (4.8)  and 
c 

differentiated  form  of  (4.6).  This  completes  the  proof  of  Lemma  4.2. 

n  -8.  t 

Bemark.  In  Lemma  4.2,  if  8(t)  ■  J  b^e  with  b^  >  0  and  0  <  Bj  <  Bj 

lt"*1  Y(0) 

then  condition  (4.3)  is  satisfied  if  B.  >  -  holds.  Indeed,  since  log  8 


and  nonincreasing,  it  suffices  to  require 


11* 


B’(t) 
B(  t) 


Y(0) 

c_ 


the 


is 


«  V 

convex 


Proof  of  Theorem  4.1.  We  begin  with  the  proof  of  the  existence  and  uniqueness  of 

2  2 

strong  solutions  of  the  Volterra  equation  (4.2).  Defining  f:  (0,*>)  *  L  (0,1)  *  L  (0,1) 

fcy 

(4.9)  f  -  u.  ♦  1*(h  -  r(8)*h  -  uft*(8)) 

w  u 

we  have 

(4.10)  f*  -  h  -  r(6)*h  -  uQr(8)  (0  «  t  <  •,  0  <  *  <  l)  . 


It  follows  from  Lemma  4.2  that  the  kernel  b  satisfies  the  assumption 


CH)  b(0)  >  0.  b  c  tc.  b*  «  BV_  t0,->  , 

loc  lo c 

12  2  5 

and  that  f  t  W^tO.-.-L  (c,1))  whenever  h  e  L^tfO,-)  *  tO.l)).  Since  under 
assumptions  to).  A  -  3 •?,  where  <t  is  defined  in  Lemma  1.3,  Chapter  1,  the  existence  and 
uniqueness  of  strong  solutions  u  with  the  properties  asserted  in  conclusion  1  of  Theorem 
4.1  -follows  from  the  result  of  Crandall  and  Nohel  (26,  Theorem  4  and  Theorems  3.5  and  3.10, 
Chapter  2] . 

We  next  establish  the  asymptotic  results  asserted  in  Theorem  4.1.  since 
r(6)  t  L*  (0,— )  (proved  in  the  demonstration  of  Lenina  4.2),  and  since  h  -  h^  +  h2  one 
has  from  (4.10)  that 

f*  -  (ht  -  r(0)*h1)  +  (h2  -  r($)*h2)  -  r(B)u0  (0  <  t  <  0,  0  <  x  <  1)  , 
where  h1  -  r(6)*h1  e  l"(0,-;L2(0,1)),  and 

«• 

11*  (h  -  r(S)*h  ){t)  -  (1  -  /  r(6)(T)dT)h*  -  s(BH-)h~  «  (1  +  B)-1h"  » 
t-*-  o  11  1 

moreover, 

(hj  -  r(B)*h2  -  r(6)uQ)  «  L1(0,-;L2(0,D)  +  LP(0,-;L2(0,1 ))  , 
with  1  <  p  <  «. 

By  using  (4.10),  the  fact  that  the  kernel  b  defined  by  (b)  is  by  Lemma  4.2 
completely  positive,  and  Theorem  2.2,  we  can  write  the  Vol terra  equation  (4.2)  in  the 
equivalent  form 

U.11)  u  ♦  b*Au  -  uQ  ♦  b*(uf*  +  k»f)  (0  <  t  <  -)  . 

To  arrive  at  (4.11)  we  use  the  relation  ob  ♦  k  •  b  «  1  in  the  right  hand  side  of  (4.9) 
and  recombine  terms  making  use  of  (4.10).  Thus  (4.11)  is  in  the  basic  form  (Vg)  of  Section 
3  with 


*4.12) 

q  -  af  ♦  k*f ' 

*0  4  t  <  -)  . 

From  Lea 

u  4.2  k  e  L^O,-). 

ind  g  -  gt  +  g2. 

where  (using  h  -  hj  ♦  hj  In  (4.10)) 

*4.13) 

»1  “  a(h1  - 

*  k*(h1  - 

r(6l*hj )  (0  <  t  <  -J  , 

(4.14) 

h  "  »2,1  ♦  ®2,2 

*0  4  t  <  -) 

vtth 
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93  1  ■  -au0r($)  -  ar(8)*hj  ♦  k»(h2  -  r(8)*h2  -  u#r(8).  gJ2  -  ah2  . 
i  •  2  d  2 

Clearly  g^  t  L  ( 0 , " ; L  (0,1))  and  g22  e  Lp(0  ,“,-L  (0,1)).  From  Lemma  4.2  one  has  that 
b  Is  completely  positive  on  (0,*),  b  f  L1(0,*>),  and  b'  e  L*  ((>,•).  Thus  all  assumptions 
of  Theorem  3.3,  part  2,  are  satisfied.  We  conclude  that  estimate  (3.7)  of  Theorem  3.3, 

•  m  m  m 

Section  3#  holds,  and  therefore  lim  lu(t)  -  u  I  *  0,  where  u  ■  A  gt  with  g  given 

t+*»  *  ’ 

m 

in  the  statement  of  Theorem  4.1)  note  that  to  evaluate  g^  use  is  made  of  (4.4)  and  of 
Proposition  2.1. 

Finally,  ve  establish  the  "comparison*'  result  asserted  in  Theorem  4.1,  part  2.  Let 

2  2 
P  ■  {u  t  L  (0,1)  :  u  >  0})  P  is  a  closed  convex  cone  in  Lz{0,1)  and  v  -  u  e  P  if  and 

only  if  u  <  v.  Moreover,  it  is  standard  that  if  u  <  v  then  J^u  <  J^v  for  every 

X  >  0,  where  ”  (I  +  XA)  *.  We  shall  prove  the  result  for  solutions  of  the  Volterra 

equation  (V^)  which  is  equivalent  to  (4.2).  As  usual  we  shall  prove  the  result  for 

solutions  of  the  approximating  equation  (v  )  of  (V,)  in  which  A  is  replaced  by  the 

Tosida  approximation  A^,  X  >  0,  and  then  obtain  the  result  by  letting  X  ♦  0. 

Let  u0 e  L2(0,1),  h ie  (0, 1 ) ,  i  -  1,2,  satisfy  uQ1  <  u02  and  h^t)  <  h2(t) 

a.e.  on  [0,T] ;  let  u^  ^  be  the  strong  solutions  of  the  approximating  equation 

,  ♦  X  ’c*u.  »  X  .  ♦  u.  ,  ♦  1*h. 

X,i  X,i  X,i  X  X,i  0,i  i 

(i  -  1,2)  X  >  0)  0  <  t  <  T)  . 


"lA1 


It  follows  by  an  elementary  calculation  (which  uses  the  definitions  of  b,  r(X-1b),  and 
the  relation  C  =  b  +  8*b)  that 

(4.15)  nX  i  -  r(X*'1b)*J^u^  ±  ♦  f^  i  (i  -  1,2;  X  >  0)  , 
where  f.  ,  are  solutions  of  the  linear  Volterra  equation 

(4.16)  *X#1  ♦  ***X>1  ♦  l  -  u01  ♦  1*hA  (i  -  1,2)  . 

Bence  by  a  familiar  calculation  one  has 

fX,i  “  u0,i"{e  *  X"1c>  *  hi*s(®  ♦  x”1<:>  U  -  1.2)  , 
and  free  (4.15),  (4.16)  the  difference  ux  2  ~  UX  1  Mtlafla*  the  Volterra  equation 
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(4.17J 


MX,2  “  UX,1  “  r(X  b,*<JXUX.2  "  JXMX,1) 


♦  C«0#j  -  ♦  X"’c>  +  (h2  -  hj)*s(0  ♦  X^C)  . 


Since  0  ♦  X  'c  Is  positive,  nonincreasing,  it  follows  from  Levin's  result  (see  Section  2) 
that  s(0  ♦  X^’c) (t)  >  0.  Thus 


*x  *  (uQ  j  •  u0  ()  s{0  +  X  ’c)  +  (hj  -  h1)*s(0  ♦  X  ’c) 


satisfies 

*^(t)  >  0  a.e.  on  (0,«). 

Next, 

define  vx  -  ux#2  -  «x>,  - 

using  (4.17)  we  have  that  vx  satisfies 

(4.18) 

vx  -  r(X"1b)*(Jx(vx 

+  UX,1  +  *X)  “  'VUX,1)>  * 

As  in  (18]  (see  also  Remark  (lii)  following  the  statement  of  Theorem  3.1,  Section  3, 
Chapter  4)  one  shows  that  -  11a  v  ^ ,  where 


(4.19) 


"l.nfl  *  r(X’1b,*<JX{vX.n  +  UX,1  +  V  *  JX(«X.in 


where  v^  ^  «  L1  (0,T;L2 (0, 1 ) )  is  arbitrary.  Choosing  ^(t)  (  P  a.e.  on  (0,T], 

T  >  0  arbitrary,  one  shows  easily  that  v.  (t)  e  P  a.e.  on  (0,T)  for  all  positive 

AyA 

integers  n.  This  also  uses  the  fact  that  r(X_1b)  >  0,  that  is  an  increasing  sup 

with  respect  to  the  ordering  <,  and  that  *^(t)  c  P  a.e.  on  (0,Tl.  Thus  v^(t)  *  P 
a.e.  on  (0,T],  and  for  X  >  Ot 

(4.20)  «x  2(t)  -  «x  ^(t)  “  *^(t)  ♦  v^(t)  >  0  a.e.  on  (0,T]  . 

Since  T  >  0  is  arbitrary,  (4.20)  holds  on  (0,*),  and  the  conclusion  follows  by  letting 
140.  This  completes  the  proof  of  Theorem  4.1. 
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Appendix  1 


B' 


Proof  cf  Lemma  2.2;  fa)  Consider  the  resolvent  equation  (k)  of  B*.  Since 
«  L*  (<>,••),  the  Paley-Wiener  theorem  1691  yields  that  k  t  I.1  (<),•*)  if  end  only  if 
P(z)  “  b{0)  ♦  B*  (z)  ■  bw  ♦  zB(z)  *  0  (Re  z  >  0)  . 


With  z  ■  x  ♦  iy 

*  '  * 

Re  P(z)  ■  bm  +  x  Re  B(z)  -  y  Im  B(z)  (x  >  0) 

*k  * 

2m  J>(z)  ■  x  1m  B(z)  +  y  Re  B(z)  (x  >  0)  . 

Since  P(z)  is  analytic  in  Re  z  >  0  and  continuous  in  Re  z  >  0,  Re  P(z)  and  Im 
F(s)  are  harmonic  for  x  >  0.  Hence  by  the  maximum  principle  for  harmonic  functions, 

4k  A 

P(z)  *  0  for  x  >  0  if  either  Re  P(iy)  “  b_  -  y  Im  B(iy),  or  Im  P(iy)  *  y  Re  B(iy) 
are  different  from  zero  for  -“  <  y  <  But  by  the  frequency  domain  condition  (F) 

Re  P(iy)  >  0  for  «*»  <  y  <  •,  and  thus  k  e  L^tO,*). 

(b)  Since  B’  «  l/fO,-)  ^L2(0,»)  and  k  e  L1  «>,••),  one  has 
2  2 

B*  *  k  ■  k  •  B*  «  t  (0,*) ,  and  the  result  k  e  !•  (0,**)  follows  by  Inspection  of  equation 
(k).  It  also  B"  *  (0,**) ,  then  B*  *  C[0,«)  (so  that  |B*  (0)|  <  •)  and  we  may 

differentiate  (k)  to  obtain 


b(0)k'(t)  «•  B*  (0)k(t)  ♦  (B"*k)  (t)  -  -  (0  <  t  <  -)  , 

and  clearly  k*  «  l'oj,**). 

(c)  If,  as  is  the  case  here  k*  «  L*  (<>,«■),  the  energy  inequality  in  (c)  is  derived 
by  the  following  simple  argument  (see  the  method  of  [67,  Theorem  1]).  Extend  k'  evenly 
for  t  <  0,  and  let 


.<t> 


|w(t>  if  t  «  [0  ,T) 
0  otherwise  . 


I 

/ 


0 


J  *  «g  • 

w(t)  *r  (k*w)  (Odt  -  k(0)  /  w“(t)dt  /  w(t)(k*«w)(t)dt 

0  0 

T  T  T 

-  k(0)  /  w2(t)dt  +  J  f  w(t)  /  k*  (t  -  T)v(T)dTdt 
0  0  0 


-  kC0>  /  w|(t)dt  ♦  j  /  vT<t>  /  k'Ct  -  T)wT(T>dtdt  . 


Letting  w^tn)  -  /  eiT,tw^(t)dt,  (n  <  R), 


the  Parseval  and  convolution  theorems  give 


/ 

0 


*<t)  ~  (k«v)(t)dt 


/  |wTCn)!2dn  “•  /  lHp(n)l2k*(n)dn  . 


But  k*  (n)  ”  2Re  k'  (in),  where  *  Is  the  Laplace  transform,  and 

as  * 

Be  ka(in)  “  Retink(irj)  -  k(0)l.  Therefore, 


/  v(t)  ^  (k*w)(t)dt  -  ~  /  |wT(r.)|2Re[ink(ln)]dn 

0  -• 


Bow  an  easy  calculation  from  equation  (k)  yields 


Hetink(in)]  -  Re 


Re  B(in> 


b(in)  .  ,  .  b. 

(Re  B(in))  ♦  (la  B(in)  -  — ) 

It 


W  Re  B(in) 


na(Re  B( in))2  ♦  dm  B(in)r>  -  b.,)2 


>  0  (-•<»<•), 


where  the  last  Inequality  follows  from  the  assumption  that  B  is  a  kernel  of  positive  type 

on  10,-)  (which  is  equivalent  to  Re  B(in)  >  0  (67,  Thau  2)r  note  that  it  is  impossible  to 

* 

bound  Retlqk(in)]  away  from  aero,  even  if  B  is  strongly  positive  on  t0,-)). 


-02- 


|d)  Multiply  equation  (k)  by  /t: 

b(0)/t  k(t)  ♦  /t  (B'nitt)  - 

An  elementary  calculation  involving  /t(B**k) 


/ t  B* ( t) 
b(0) 


CO  <  t  <  •)  . 


shows  that  /t  k  satisfies 


b|0)^t  k( t)  ♦  /  B'Ct  -  t)/x  k|T)dx 

0 


~  ~  /  (/t  -  /t)B‘(t  -  T)k(T)dt 


(0  <  t  <  •)  . 


Since  /t  -  /r  <  /t  -  T  for  0  <  t  <  t,  and  since  also  /t  B*  «  L^(0,»)  by  assumption 
and  k  f  L1 (0,«)  by  (a),  the  integral  on  the  right  side  of  the  last  equation  defines  a 
function  in  I,'  (0 ,*) .  Then  /t  k  «  L 1  (0 ,“)  by  the  argument  of  part  (a).  The  additional 
assumptions  and  elementary  estimates  also  yield  /t  k  *  L2(0  ,«•) . 

Finally,  differentiating  the  equation  (k)  as  in  part  (b),  multiplying  the  resulting 
equation  by  /t,  and  using  elementary  estimates  yields  /t  k*  e  L^IO,*).  This  completes 


the  proof  of  Lemma  2.2 
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